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Section 1 - Summary 


MOTIVATION FOR THE VECTOR MAGNETOMETER. STUDY AND ITS 
TASKS 

The study was undertaken to establish the feasibility of developing a 
satellite -borne magnetometer required to make a more accurate map of 
the earth's magnetic field. 

In the late 1970's a Scout-launched Application Explorer Satellite, the 
MAGSAT, will carry advanced vector and scalar magnetometers into earth 
orbit for the purpose of making a high accuracy global survey of the earth's 
magnetic field. 

To obtain high resolution, a low polar orbit will be used. To prevent 
temperature variations and to obtain uniform atmospheric conditions, the 
nominal orbit will be sun- synchronous to the plane of the earth's ter- 
minator — a dawn-to-dusk orbit. 

The data from this mission will provide an updated quantitative global 
map of the earth's vector magnetic field. This data will be used to revise 
the regional magnetic charts used for navigation and other purposes through- 
out the world. 

In addition, the data will indicate the presence of magnetic anomalies in 
the earth's crust. When correlated with geological and geophysical models 
of the crust, this information can indicate potential targets for natural 
resource exploration. 

The overall accuracy requirements desired for the vector magnetom- 
eter instrument are high. Not only must the instrument measure the vector 
field with an absolute error per axis of 5Y and an alignment stability of 
5 seconds of arc, but the orientation of the instrument with respect to the 
earth must also be precisely known. Because an instrument with such accu- 
racy and alignment stability has never been flown, the present study was 
undertaken under a contract fromNASA-GSFC to establish the feasibility 
of such an instrument. The contract required the development of a design 
approach to the problem of precise measurement of vector field aboard a 
low-altitude, earth-orbiting satellite, with analytical verification that the 
design approach developed is capable of meeting the requirements of the 
specification. Consideration of the problem of calibrating the instrument 
was also part of the specified effort. 

The tasks actually performed in the study are indicated in Table 1-2. 

A substantial effort was devoted to the problem of selecting the most prom- 
ising sensor configuration from the standpoint of alignment stability. Once 
this configuration was selected, a specific structural and thermal design 
for the sensor was developed and analyzed. 

The need for high accuracy and stability imposes severe requirements 
on the magnetometer electronics. Most of the study was devoted to the 
problem of accurate measurement and digitization of the signals produced 
by the magnetometer. 

Calibration and alignment of the magnetometer to the required accuracy 
poses a difficult problem. One can establish the location of an optical ref- 
erence with second-of-arc accuracy, but it is not as easy to find the direc- 
tion of a reference magnetic field that accurately. This problem was also 
studied, with emphasis on the actual experimental situation and the method 
of reducing the data to obtain the most reliable scale factor and alignment 
calibration. 


1-2 



TtMt 


zero 


( 


The study produced a feasible design approach, as well as a number of 
efimtc conclusions and some recommendations for the development of the 
flight magnetometer and for extensions of the study to resolve some problem 
areas uncovered. A summary of the design approach is presented next. 

TABLE 1-1. OBJECTIVES OF THE MAGSAT MISSION 


• Provide updated global magnetic field maps. 

• Locate magnetic anomalies in the earth's crust. 


TABLE 1-2. TASKS PERFORMED IN THE VECTOR 
MAGNETOMETER STUDY 


The purpose of the study was to e 
vector magnetometer with an absolute 
stability of 5 seconds of arc. 


stablish the feasibility of a flight 
accuracy of 5Y and an alignment 


Task 

• Sensor Configuration 


• Sensor Design 


i • Elect r onic s 


• Calibration 


Scope 


- Establish the sensor configuration 
least sensitive to misalignment 
errors. 


- Select a specific sensor design 
approach and confirm its perform- 
ance by analysis. 

Select a specific electronics config- 
uration and verify its performance 
by analysis. 

- Develop a plan for calibrating the 
magnetometer to the required 
accuracy. 
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DESCRIPTION OF THE DESIGN APPROACH 

The Vector Magnetometer design studied uses a triaxial fluxgate sensor 
assembly with triaxial vector feedback to eliminate the effects of fluxgate 
misalignment. The fluxgate electronics has been specifically designed for 
superior stability using currently available tec hnology. 

The requirements for the Vector Magnetometer for MAGSAT are 
among the most stringent ever conceived. During the study we selected 
a design and demonstrated by analysis that it exceeds the requirements - 
in most cases — by substantial margins. The principal features of this 
design are listed opposite and summarized below. 

Pr incipal Features of the Design - To eliminate the effect of sensor 
misalignments, three orthogonal ring-core fluxgate sensors are operated 
within a coil system which provides a triaxial vector feedback field that 
cancels the earth's field, so that the sensors operate in a null-sensing 
mode. When this is done, the alignment of the system depends almost 
solely on the alignment of the feedback coils and not on the alignment of 
the sensor elements, which are difficult to make dimensionally stable. 

Also, the sensors contribute only to the offset and noise of the system an 
not to instabilities in the overall gain. This approach was used on the 
Pioneer 10/11 Vector Helium Magnetometer. 

An exploded view of the sensor is shown in the figure opposite. The 
feedback coils are constructed on spherical graphite-epoxy forms, which 
provide high structural rigidity and thermal conductivity. Analysis of this 
construction shows that all alignment and stability requirements can be 
met with substantial margins. 

An external set of coils provides essentially complete cancellation of 
the external field (to 0. 027 at 60 cm from the center of the sensor). 

The sensor housing is made almost entirely of non-conducting materials 
to eliminate thermo-electric currents. The sensor assembly is thermally 
insulated and operated at very low power, to control self-heating effects. ^ 
The multiple -layer thermal shield consists of insulating layers and a shunt 
conductive layer to conduct heat to the temperature-controlled baseplate. 
Active temperature control of the baseplate using an ac-dnven printed- 
circuit heater reduces internal gradients. 

The electronic design follows a conventional arrangement. 1 he llux- 
eate cores are excited by a tuned drive circuit to provide high excitation 
with low drive losses. The second harmonic signal from the sensor is fil- 
tered and amplified, synchronously demodulated, and further amplified to 
provide a dc voltage. This voltage is converted to a current which is fed 
back to the triaxial coils. The circuits are implemented using modern 
operational amplifiers to provide the required stability. The critical cur- 
rent feedback components must be thermally stabilized. 

Conversion of the dc voltage to a digital signal to be telemetered may 
be accomplished by either a single-range 17-Bit ADC or a multiple -range 
12-Bit ADC with an offset field generator. Only the detailed design of this 

latter approach however was investigated. 

Performance - In the course of the study, we performed an analysis 
of this design based on a detailed error allocation. Our analysis shows 
that the MAGSAT requirements can be exceeded by the proposed design, in 
most cases by substantial margins. The discussion of the specification 
compliance, presented next, indicates in detail the level of performance 
that could be achieved. 
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TABLE 1-3. PRINCIPAL DESIGN FEATURES 


Feature 

• Triaxial Vector Feedback 

• Rigid Spherical Feedback 
Coil Forms 

• Multiple Thermal Shield 

• Active Temperature Control 

• Current Source Temperature 
Control 

• Dual- slope Integrating ADC 


Advantage 

- Eliminates effects of sensor 
misalignment. 

- Provides high thermal and struc- 
tural stability. 

- Very low internal temperature 
gradients. 

- Ensures fixed-temperature 
operation. 

- Ensures gain stability and linearity. 

- Eliminates zero offset errors. 



Figure 1-1. Exploded View of the Sensor. Triaxial 
feedback is provided by coils mounted on rigid ther- 
mally stable spherical forms. 
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COMPLIANCE WITH STUDY CONTRACT SPECIFICATION 

The design approach was selected after an optimization analysis. Analysis 
of the performance shows that it will meet or exceed the specified perform- 
ance, within the specified weight, power, and environmental constraints. 

The contract for the present study specified the performance require- 
ments of the Vector Magnetometer, and required that the performance of 
any design approach be confirmed by analysis. 

Performance Requirements - The specified orthogonality and offset 
requirements are all similar to those of magnetometers we and others have 
built and flown before, and we see no difficulty in meeting them. 

The noise requirement can be met, but we found that the peak noise 
requirement in the specification is unnecessarily stringent in view of the 
fact that the specified quantizing error is ten times larger. With the agree- 
ment of NASA technical personnel, we adopted a modified noise specifica- 
tion — the highest noise we thought was acceptable. By doing this, we pro- 
vided higher design margins for meeting the critical requirements for 
absolute accuracy and alignment stability. A compromise design with lower 
noise would make the sensor larger and impact the weight. 

Selecting an Approach to Exceed the Requirements - Our approach to 
the study has been to exceed the MAGSAT requirements rather than merely 
come close to them. We knew when we started the study that the use of tri- 
axial feedback would enable us to have a good chance of meeting the align- 
ment stability requirements, but we did not expect to be able to exceed them 
by a large margin. In fact, in our analysis of the triaxial feedback design, 
we have been unable to identify an error which we cannot control by design 
to well within the specification. The analysis is presented in the body of 
the report. 

The analysis presented does not consist only of a detailed analysis of 
the specific design which was studied. To select the design approach and 
narrow down the range of possible design characteristics, we performed 
tradeoff analyses to show that the approach selected was actually the best 
approach, or very near to it. 

The principal features of the design which meets these specifications 
have already been described. In Table 1-4 we have identified those design 
features which enable us to meet each requirement. 

Analysis of Performance - As explained previously, our ability to 
exceed these requirements is a direct consequence of two design features — 
the use of triaxial feedback and the fabrication of spherical feedback coil 
forms and other structural components from graphite-epoxy materials. 
These materials provide the required high rigidity and high thermal con- 
ductivity, but yet are light in weight and do not appear to exhibit any signi- 
ficant thermoelectric effect. 

The analysis which we have performed has put us in position to make a 
number of specific recommendations, which are treated in the following 
pages. 
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TABLE 1-4 COMPLIANCE WITH STUDY CONTRACT SPECIFICAH ON 


Zero Offset 
Sensor 
Electronics 
RSS Total 


Required 

Performance 


0. 47 


Estimated 

Performance 


0. 37 


0.47 


Sensor 
Electronics 
RSS Total 

Scale Factor 
Sensor 
Electronics 
RSS Total 

Sensor Alignment 

External Field 


1. 07 


4. 17 


0. 57 


0. 35 


0. 357 


0. 06 


1. 17 


2. 5 arc sec 0. 6 arc sec 


0. 01 7 


assembly and on Page 5-13 for the electronic 
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PRINCIPAL CONCLUSIONS OF THE STUDY 

The vector magnetometer should employ triaxial vector feedback to reduce 
misalignment sensitivity. By use of a strong feedback coil structure 
made from isotropic, stable material, and by use of suitable thermal 
shielding, the stability and alignment of the sensor can be made to exceed 
requirements. Either a single-range analog-to-digital converter (ADC) 
system or one with multiple offset ranges will provide the electronics 
stability required. . 

While the noise and offset requirements for the magnetometer have been 
met many times before, the scale factor and alignment stability and cali- 
bration are the most stringent ever required. The principal conclusions 
reached in our study of how these requirements can be attained are listed 
opposite. 

The scale factor stability of the magnetometer is determined by the 
stability of the sensor feedback coils and the feedback electronics. The 
alignment can depend not only on the alignment of the feedback coil, but 
also on the difficult-to-control alignment of the magnetic sensor itself. 

Sensor Design Study: The Need for Triaxial Vector Feedback - The 

recommended design approach employs fluxgate sensors because they are 
simple and meet the noise requirements. To attain the required scale 
factor stability and linearity, these sensors must be enclosed in an overall 
magnetic feedback loop, cancelling the earth's field with a feedback coil 
around the sensor so that the sensor operates in a null mode. Then, only 
the sensor noise and offset matter. In our analysis, we found that the sen- 
sitivity of the overall alignment to misalignment of the sensors themselves — 
in contrast to the feedback coils — could almost be eliminated by nulling all 
three components of the vector field at each sensor. Mathematically, this 
finding relates to the fact that the orientation of a null vector (of zero 
length) need not be specified. Thus, our specific design approach consists 
of a system in which three orthogonal sensors are surrounded by three 
orthogonal feedback coils. 

In investigating this system, we found that the external field can be 
cancelled almost exactly (to <0. 02 Y at 60 cm) by use of external cancella- 
tion coils. By careful selection of the coil geometry and the use of a 
simple, rigid structure using isotropic graphite-epoxy materials, as well 
as by providing thermal shielding, we were able to obtain an estimated 
scale factor and alignment stability of 1. 5Y and <5 seconds of arc, respec- 
tively. The weight and power requirements of the sensor are low: 1. 82 kg 

and 300 mW, the power being almost entirely that of the thermal control 
heater. The feedback coil and sensor drive power is only 35 mW maximum. 

Electronic Design Study - Because low sensor power input is required 
for thermal stability of the sensor housing, a reactive drive scheme is 
needed to eliminate the losses produced by conventional current or voltage 
drive methods. Except for the use of this technique (to be flown on MJS) 
the sensor drive and ac signal handling circuits can be conventional. The 
important problem is the design of the feedback system and the ADC. The 
study found that either a single- range ADC or a limited- range ADC with 
several offset ranges could meet the requirements; the design approach 
studied in detail was based on the multiple-range configuration. 

To achieve the necessary stability of the feedback coil current, it is 
necessary to control the temperature of the feedback current source to 
±5°C. 
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Calibration Study - The study concluded that the absolute determination of 
the scale factor and alignment of the instrument is very difficult. The veri- 
fication of axis and scale factor stability is much less difficult. There 
appears to be no difficulty in principle in determining the orientation of the 
magnetometer housing to the order of one second of arc by optical means. 
However, the real problem is to determine the actual physical orientation of 
the sensor magnetic axes. This cannot be done without knowledge of the abso- 
lute orientation of the calibration field — not of the field-producing coils, but 
the field itself. Only the relative orientation of the coils at the GSFC facility 
is known, and that only to about 10 seconds of arc. If this facility is to be 
used, it will have to be upgraded. 

TABLE 1-5. CONCLUSION OF VECTOR MAGNETOMETER STUDY 


Conclusions of Sensor Design Study 

• Triaxial feedback is needed to reduce the misalignment sensitivity. 

• Large external coils are needed to attain a low external field with 
low power and good field uniformity. 

• Only a small coil structure is needed (less than allowed envelope) 
even if only two coils are used. 

• By using spherical coil forms made with graphite-epoxy material, 
we can get a field stable to 1, 5Y and <5 seconds of arc. 

• Use sensor temperature control to ±1°C. 

• Require electronics temperature control to ±5°C. 

Conclusions of Electronics Design Study 

• Reactive sensor drive required to reduce sensor dissipation. 

• Either a single-range or multiple -offset- range ADC will work. 

• Electronics temperature control required (±5°C). 

Conclusion of Calibration Study 

• Principal problem is absolute calibration. 

• Upgrading of GSFC calibration facility is required. 
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Section 1 - Summary 


RECOMMENDATIONS FOR HARDWARE DEVELOPMENT AND 
FURTHER STUDY 

The flight instrument should definitely employ a triaxial feedback scheme 
using a graphite-epoxy structure for stability. The details of the design 
should be checked by measuring prototype fluxgates and feedback coils. 

The theory of the fluxgate performance should be developed further, and 
the calibration problem needs further work. — — 

From the conclusions reached in the study, we have developed definite 
recommendations for the next steps in the development of the vector mag- 
netometer. Some recommendations relate to the design of the hardware, 
while others relate to the attainment of better understanding of how to 

design and test such a device. , . 

Hardware Recommendations - The principal conclusion of the study is 
that triaxial vector feedback should be used to attain low sensitivity to flux- 
gate sensor misalignment by using a three- axis coil system surrounding 
three orthogonal fluxgate sensors. This follows directly from a mathemat- 
ical analysis of the sensitivity problem. The alternative the use of ree 
separate single- axis sensors mounted at right angles - does not provide 
the required low axis alignment sensitivity because the feedback loop for 
each single-axis sensor is closed in only one direction. , 

The size and detailed design of the triaxial feedback coils depends upon 
an engineering analysis in which the size of the fluxgate sensors, and hence 
the feedback coils, was selected by optimization of the fluxgate and coil 
dimension to minimize the size and weight while meeting a somewhat 
relaxed noise specification. Before committing the design to a hardware 
development phase, the assumptions should be checked by building and test- 
ing sensors and coils with the recommended characteristics to verify the 

It is generally agreed, and confirmed by our work, that the structural 
and thermal stability required can be attained if a graphite-epoxy structure 
is used. This material has the combination of high rigidity, low thermal 
expansion, and high thermal conductivity characteristics necessary to pro- 
vide the required structural and thermal stability. We have recommended 
the use of an isotropic material with random fiber orientation to avoid 
thermal distortion caused by anisotropic expansion. The use of this ma e- 
rial should be investigated further with the fabrication and test of a proto- 

tvpe structure. . . , 

Although the electronics design required is conventional and straight- 
forward, we do recommend that the critical subsystems such as the 17-bit 
ADC for the single-loop system as well as the DAC, OFG, and Feedback 
Current Generator for the two-loop system be studied further to verify the 

flight worthiness of the components used. 

Recom mendations for Further Study - Other than the excellent metal- 
lurgical work reported by NOL, the literature search performed for the 
study did not uncover any published analyses of the problem of obtaining 
optimum performance of the sensor - i. e. , analyses of the signal-to-noise 
ratio as opposed to separate calculations of the signal sensitivity or the 
noise. Consequently, we developed analyses of this area ourselves. 

Our analyses of sensor performance contradict some conventionally 
accepted ideas about fluxgate design (although we understand that results 
in agreement with our analyses have been obtained experimentally). 


1-10 





was already realized- if on^i^to ^ P^° blem em P hasized a difficulty which 
rate facilities are ne^ed to o thlS m ° re accu ' 

s\7w^tThe r ° aCh ^ 

for those ^ses in whi<fht^ t Derf^rnifi? Uir f e ^i^ )er ^j^ manCe iS attained > except 
which the design requirements are impLmented. ^ ^ mannCr in 



• Fabricate and test graphite-epoxy structure. 

# s E y“^ ^ e rs r . thineSS ° f Si ^ le - lo °P and ‘—loop electronics 



# facn^y CaUbration 6tud y to delude plans for upgrading GSFC test 



I 


1-11 





SECTION 2 

MAGNETOMETER DESIGN REQUIREMENTS 

Review of Vector Magnetometer Performance 
Requirements 


Components of the Vector Magnetometer 
Allocation of the Errors to Sources . . . . 



Section 2 - Magnetometer Design Requirements 


REVIEW OF VECTOR MAGNETOMETER PERFORMANCE 
REQUIREMENTS 

The Vector Magnetometer for the Applications Explorer Mission (AEM) 
will be used to obtain a global map of the earth's magnetic field. To pro- 
vide a significant improvement over existing maps, the magnetometer 
must be exceedingly stable, and its orientation precisely known. 

In the following paragraphs we will discuss briefly the requirements 
for the Vector Magnetometer. The remainder of this study addresses the 
problem of how these requirements can be met. The requirements are 
summarized in the Table. 

Need for Accurate Measurements - To make measurements which 
advance the state of the art it is necessary to measure the earth's field 
with an overall accuracy of 5Y (IV = 10“9 Tesla). This requirement is 
difficult because there are a number of error sources which contribute to 
inaccuracy and because the earth's magnetic field can be as large as 
64000Y, so that the required accuracy is better than 1 part in 10 4 . The 
contribution to the inaccuracy from stability in the instrument offset must 
be less than ±0. 4Y over the operating temperature range; the digital reso- 
lution must be ±1Y. Of these requirements, the one for offset stability is 
the most difficult to achieve, since the allowable drift is only six parts in 

io 6 . 

Noise and Bandwidth Considerations - The bandwidth of the instrument 
is to be 25 Hz, with a capability of limiting the bandwidth to 1 Hz. Over 
the 25-Hz bandwidth, the noise is specified to be 0. 1Y peak for a 10-second 
scan with a 25-Hz bandwidth (about ±0.038Yrms, or less). It has been 
agreed that this requirement is unnecessarily stringent in view of the fact 
that the quantizing error is lYpeak. For the study, we have adopted the 
criterion that the peak noise should be at most equal to the peak quantizing 
error. In addition, the rms noise should be at most 0. 38Y. 

Physical Characteristics - The total weight of the magnetometer, 
including the sensor electronics and the signal and power cables in the 
boom, is to be less than 5 kg. The sensor volume is limited to a 30-cm 
diameter sphere, while the electronics volume is limited to 4500 cm as 
shown in the table. 

Importance of Alignment Stability - The overall accuracy of the field 
measurement is affected by the misalignments of the sensor axes. To 
provide the required accuracy the orientation of the axes must be known 
with an angular error of less than 7 8 |jl radians or 16 seconds of arc; this 
angular misalignment produces a signal of 5Y along the axis of the mag- 
netometer normal to the full scale 64000-Y field. To measure the orienta- 
tion of the sensor — mounted on a 3-6 meter flexible boom - an optical 
attitude determination system will be used. A number of error sources 
contribute to the total error in determining the attitude of the instrument. 
The budget for the errors internal to the instrument is shown in the table. 

The orientation of sensor axes is to be repeatable over the temperature 
and environmental range to 5 arc-seconds repeatable. The magnetic axes 
shall be orthogonal to within ±0. 1 degree. This requirement also seems 
unnecessarily stringent. Even with a maximum inaccuracy of 5Y, one 
could correct for lack of orthogonality as large as 1 degree with a worst- 
case error of only ±0. 1Y. 
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The angular alignment between the sensor and the geometric axes of the 
sensor housing, as well as the alignment between the optical reference and 
the geometric axes of the housing, must be measurable to 1 arc-second. 

Since the error can presumably be in either direction, we interpret this 
specification to mean ±1 arc-second for each of the two error sources, for 
a worst-case error of ±2 arc-seconds. 

gxternal Field Cancellation - Because a scalar magnetometer is mounted 
nearby, the external field of the Vector Magnetometer must be small. We 
have adopted a value of 0. 5Y at 60 cm (24 inches) from the center of the 
Vector Magnetometer sensor. 

Environmental Constraints - The performance specifications of the mag- 
netometer are to be met before and after exposure to the shock, vibration, 
and acceleration tests, as well as during exposure to the simulated or actual 
orbital thermal environment. The details are given in the instrument speci- 
fication. In the table we have shown the temperature and the maximum accel- 
erations experienced in test. 

• these difficult requirements necessitates careful systems analy- 

sis of the instrument design. We now proceed to discuss the implications of 
these requirements for the design of the magnetometer. 

TABLE 2-1. PERFORMANCE REQUIREMENTS FOR VECTOR MAGNETIC 
FIELD MEASUREMENTS 


The requirements for accuracy 

and alignment are particularly stringent. 

Measurement 

3 orthogonal components of earth's 
magnetic field 

Orbit Characteristics 

Near-polar, low altitude (325-550 km) 
Insolation: 60 min; eclipse 30 min. 

Dynamic Range per Axis 

±64000Y 

Overall Accuracy 

±sy 

Offset Stability 

*0. 4Y 

Resolution 

±1Y 

Bandwidth 

25 Hz (limiting to 1 Hz) 

Noise 

0. 1Y peak for 10 second scan, with 25-Hz 
bandwidth (original specification); 1Y peak 
or 0. 3Y rms (modified specification) 

Axis Stability 

5 arc-seconds 

Axis Orthogonality 

90 degrees ±0. 1 degree 

Axis Determination Accuracy 

Sensor Elements to Sensor Axes: ±1 arc- 
second; Optical Reference to Sensor Axes: 
±1 arc-second 

External Field 

0, 5Y at 60 cm (24 inches) 

Boom Characteristics 

3-6 meter extensible 

Total Weight 

5 kg 

Sensor Volume 

30- cm diameter sphere 

Electronics Volume 

20 cm x 15 cm x 15 cm 

Power Consumption 

6 W 

Temperature 

Storage: -50°C to +85°C 
Thermal Shock: -50°C to +85°C 

• 

Electronics Operating: -10°C to +50°C 

Acceleration 

Sinusoidal Vibration: 20 £ (36-158 Hz) 


Random Vibration: 12. 9^-rms (200-400 Hz) 
Thrust Acceleration: 22. 5 £ 
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Section 2 - Magnetometer Design Requirements 


COMPONENTS OF THE VECTOR MAGNETOMETER 

A vector fluxgate magnetometer has three separate sensor channels, one 
for each axis. The characteristics of the sensor, signal electronics, 
feedback circuits, and digitizing circuits must be appropriately selected 
to meet the overall performance requirements. 

To design the Vector Magnetometer, it is necessary to translate the 
overall performance requirements just described into design require- 
ments for each component of the magnetometer. 

The requirements for klAGSAT can be met by use of either a vector 
helium or a three-axis fluxgate magnetometer. The vector helium mag- 
netometer provides better low-frequency noise performance at the cost of 
greater complexity; because this particular low-noise characteristic is not 
needed for the MAGSAT application, a fluxgate magnetometer has been 
selected for study. The components which typically make up a fluxgate 
magnetometer are shown opposite, for one of three sensor channels; the 
central power supply, sequencing logic, and other housekeeping and sup- 
port systems are not shown. 

How the Fluxgate Operates - The fluxgate sensor consists of a magnetic 
core consisting of either a pair of long rods or a tape-wound toroid of a 
high permeability square-loop core material. It is driven into saturation 
at a high audio frequency (typically 10 kHz) by current in a drive winding. 

In the presence of an external field, an even-harmonic flux change is gen- 
erated in the core as a result of a variation in the saturation time of the 
different parts of the core. This signal is amplified, synchronously 
demodulated, filtered, and fed back as a current input to a feedback coil 
which produces a field to null the external field. The sensor then func- 
tions only as a null- sensing device and does not have to meet high linearity 
or gain stability requirements. The signal which drives the feedback coil 
is digitized to provide the output data which is a number directly related 
to the value of the external field. 

To implement the block diagram shown, it is necessary to decide on 
the organization and performance characteristics for each block. 

Organization of the Feedback Arrangement - The field along each axis 
of the magnetometer is measured by a separate sensor with separate elec- 
tronics as shown. Often, the sensors for each axis are physically sep- 
arated, and use one coil both to sense the signal and to provide feedback 
to cancel one component of the applied field. Alternatively, a set of orthog- 
onal coils can be used to null the field over an entire volume occupied by 
three sensors. In either case, the three sensor electronics channels func- 
tion independently and have the same block diagram; but it is necessary to 
select the feedback arrangement which provides the best performance. 

Organization of the Digitizing System - For some applications the 
dynamic range of the applied field is very large. For MAGSAT, the ratio 
of the total field variation to the quantizing interval is 128000V/ 1Y, or 
slightly less than 2 17 . Covering such a wide dynamic range with one 
analog-to-digital converter can be avoided by using a sampling conversion 
within the feedback loop or by offsetting the feedback current by a switched 
constant current to optimize the overall digitizer performance. The selec- 
tion of the appropriate digitizing scheme is another question to be decided. 

In addition to these two organizational questions (considered in detail 
in the remainder of this report), one must specify the performance require- 
ments for each component of the magnetometer. This allocation of errors 
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can be done by identifying and specifying the error contributed by each com- 
ponent or each block in the block diagram. We now describe the allocation 
adopted for the present study. 
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Figure 2- 1. Block Diagram for One Axis of a Vector Magnetometer. 
Each of the three orthogonal sensors has the same block diagram. 
The feedback coils and the digitizing system may each be organized 
in several different ways. 




2- Magnetometer Design Requirements 
ALLOCATION OF THE ERRORS TO SOURCES 

A detailed error budget has been generated to control the design. Worst- 
case errors have been allocated to each error source, adding the effects 
of the separate errors by the root-square-sum method. 

As indicated in the discussion on meeting the performance require- 
ments, there are a number of error sources which contribute to each 
specified tolerance. To control our design, we have generated the error 
budget shown on the facing page. Because there are many error sources 
contributing, we have distributed the errors using the root-square-sum 
methods, taking the square root of the sum of the squares of the independ- 
ent worst-case errors to arrive at the overall tolerance. 

We now indicate briefly the rationale behind the allocation shown. 

Absolute Accuracy - The specification allows an overall error from 
all sources of ±5y for each axis. The most significant sources of error 
that attribute to this are, scale factor drift, non-linearity, and mis- 
alignment of the sensor. For these we have allowed the largest sensor 
budgets. In addition to these measurement errors result from offset 
uncertainity, noise, and quantization. We estimate the RSS total result- 
ing from these to be ± 4. 7y. Several of the above errors originate in 
both the sensor and the electronics. These are discussed below. 

Scale Factor - The contributions to the scale factor error attached 
to the sensor and the electronics are 1. 5y and 2. 5y for short term 
instability and 2. 5y for the long term drift of the electronics. Some of the 
sources causing this error are variations in the feedback coil constant; 
current source transimpedance, fluxgate translation error, forward loop 
amplification variations and feedback coil impedance variations. 

Linearity - The sources of non-linearities are primarily within the 
sensor"! Therefore of the allowed 1. 5*y, most has been allotted to the 
sensor. Only 0. 1 y has been allowed to errors caused by forward loop 
amplification waveforms. 

Zero Offset - There are many factors that contribute to the offset 
error. Some of these are sensor magnetostrictive, current source 
amplifiers, drive harmonics and supply voltage variations. However 
these are controllable, therefore, only ±0.4yhas been allotted, with 
±0.33\ being attributable to the electronics and ±0. 2y to the sensor. Both 
for short term enviromental stabilities and long term drifts. 

Given the allocation of allowable errors to the various sources in 
the sensor and the electronics we, now turn to a discussion of the over- 
all magnetometer systems design. 
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TABLE 2-2 DESIGN ERROR BUDGET FOR THE VECTOR MAGNETOMETER 


Absolute Accuracy Error 

Offset Uncertainty ±0.47 

Offset Long Term Drift ±0.47 
Scale Factor Drift 50 PPM 

Scale Factor Long Term 
Non-Linearity 50 PPM 

Misalignment 5 arc-sec. 

Noise ±1.07 

Quantization Error ± 1.07 

RSS Total 


Field 

Error (±7) 


Zero Offset 

Instability (±7) 

Agini 

k(± 

Sensor 

0. 20 

0 . 

20 

Electronics 

0. 33 

0, 

33 

RSS Total 

0. 39 

0 . 

39 


Scale F actor 


Le r actor Instability (±7) 

Sensor Scale Factor 1.5 

Electronics Scale Factor 2. 5 

2.9 


;ing (±7) 


Non- Linearity 

Sensor Heating 
Electronics 


Instability (±7) 
1. 5 
0 . 1 
1. 5 
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Section 3 - Magnetometer Systems Design 
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MEETING THE PERFORMANCE REQUIREMENTS 

The performance requirements can be met by the use of sensor feedback 
and thermally stable feedback circuits and coils. Both sensor and elec- 
tronics will require temperature control; low thermal gradients in the 
sensor can be attained by use of superinsulation and low coil power. 

Magnetometers with offset and noise performance equal to that required 
have been previously flown, but the requirements for absolute accuracy 
and alignment present new problems. The following qualitative discussion 
indicates what design approaches are necessary to meet the performance 
requireme nts . 

Gain Stability and Linearity - The accurate measurement of the field 
along each axis of the magnetometer requires very high gain and offset 
stability and linearity. To achieve the required gain stability and linearity, 
the magnetic field sensors (of whatever kind) must be enclosed in feedback 
coils so that the sensor operates only as a null detector. Since the field 
produced by a coil is a function of its size, thermal expansion of the coil 
becomes a critical matter, requiring temperature control of the sensors 
as well as the use of thermally stable materials. To generate the feed- 
back field, the feedback coils must be driven by a precision current source 
whose output impedance is much higher than the unstable coil impedance. 
The overall loop gain must be high. 

Offset Stability - The attainment of the required offset stability requires 
low offsets in the sensor, the ac drive and signal circuits, the dc feedback 
system, and the analog-to-digital converter (ADC). Fluxgate sensors with 
offset stabilities and noise performance comparable to these required have 
been built by Time-Zero Laboratories and others using core materials with 
low magnetostriction. The recommended sensors use the 6-81. 3 Mo- 
Permalloy recommended by the Naval Ordnance Laboratory. The elec- 
tronic offsets arise from several sources, namely, generation of second 
harmonic signals in the drive circuit and signal amplifiers, and dc offsets 
in the demodulator-filter, current source, and ADC. These can be con- 
trolled by the precautions indicated in the table. 

Axis Stability - Attainment of the required axis stability requires 
stability of the sensor mounting structure and internal orientation. To 
achieve such stability, the sensor subsystem must be kept at a constant 
temperature and must be designed for low thermal gradients in the struc- 
ture, which in turn implies minimum power input to the sensor. To avoid 
thermoelectric or induced currents, the sensor housing must be made of 
relatively non-conducting materials which possess sufficient thermal sta- 
bility to stay in alignment in spite of the internal thermal gradients pro- 
duced by the sensor. Finally, it is important to arrange the feedback so 
that the alignment depends only on the feedback coil geometry. The con- 
ventional single -axis fluxgate sensor is feedback-stabilized along only one 
axis. Motion of the sensitive magnetic element normal to this axis pro- 
duces a first-order misalignment error for fields normal to the feedback 
axis, even if the orientation of the feedback and pickup coil is stable. 

This difficulty can be circumvented by enclosing all three sensors in 
a single volume in which the vector field is nulled, rather than nulling the 
three components-separately at different locations. The only critical 
assembly is now the feedback coil structure which is less critical than 
for a single-axis system because there are fewer parameters to be con- 
trolled and the structural and thermal design is easier. The external field 
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TABLE 3-1. HOW THE PERFORMANCE REQUIREMENTS CAN BE MET 


Each performance 
requirements. 


requirement leads to a number of critical design 


• Scale factor stability 
and linearity 


• Offset stability 


• Axis stability 


• Axis determination 


Sensor feedback and sensor temperature 

- Thermally stable feedback coils 

- Low sensor power input 
Precision feedback current source 

feedback source impedance 

- High loop gain 

- Accurate ADC scale factor 


control 


- Low magnetostriction fluxgate cores 

- High symmetry core drive 
Narrowband second harmonic filters 

- High gain, linear AC amplifiers 

- Low offset current source 

- Low offset ADC 


- Triaxial feedback 

- Constant sensor temperature (1°C) 

- Low thermal gradients ( 1 °C) 

- Minimum sensor power 

- Stable graphite -fiber materials 

- Stable optical reference 

- Optical access to sensor reference surface 

- Precision calibration system 
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Section 3 - System Design 


FUNCTIONAL DESCRIPTION OF THE SYSTEM DESIGN 

The vector magnetometer design consists of three orthogonal fluxgate 
sensors mounted inside triaxial feedback coils. The sensor is enclosed 
in a housing insulated with a thermal blanket; active temperature control 
is accomplished with a printed circuit heater. A conventional one-loop 
system or two-loop system with an offset generator may be employed for 
the electronics subsystem. _ 

The design presented here is a result of analysis of the sensor and 
electronics design. The principal design decisions and their rationale are 
indicated in the table opposite, while the figure shows how the system is 

implemented. . 

Sensor Arrangement - As explained previously, it is necessary to 
enclose the sensor coils in the triaxial feedback coils to eliminate effects 
of sensor misalignment. This approach results in a true vector instru- 
ment. The alignment stability then depends almost entirely on the stability 
of the feedback coils. Analysis shows that the feedback coils can be built 
with the required stability. An extra set of coils is used to cancel the 
external field. 

Alignment - The alignment of the sensor will be measured by use of an 
alignment mirror or other optical alignment surface mounted on the sensor 

base plate. . 

Thermal Control System - It is exceedingly important to inhibit gradi- 
ents within the sensor assembly to avoid thermal expansion and distortions. 
To attain this uniformity, the sensor is insulated from the sun and from 
the cold of space by use of two multilayer super-insulated thermal blankets, 
separated by a graphite-epoxy housing which acts as a thermal shunt to the 
controlled baseplate. 

Electronics - The circuitry required is conventional. lhe second nar- 
monic signalTrom each fluxgate sensor is amplified by a tuned amplifier, 
synchronously demodulated, amplified and filtered. The amplified, fil- 
tered output is used to drive the feedback coil and is also converted to a 
digital signal in an ADC. An offset field generator (OFG) supplies discrete 
steps of additional feedback as determined by its logic controller. The 
entire process operates synchronously with signals from the spacecraft. 

The critical problems in the design of this system are the mechanical 
and thermal design of the sensor and the design of the stable feedback 
system. 
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TABLE 3-2. MAIN FEATURES OF VECTOR MAGNETOMETER DESIGN 


F eature 


Rationale 


Fluxgate sensor 
Triaxial feedback 
Switched offset generator 
Spherical housing 
Thermal blankets 

Printed Circuit Heater 


- Low power, weight, volume 

- Insensitive to fluxgate alignment 

- Stable offset over dynamic range 

- Reduced thermal radiation 

- Reduced effect of insolation and radiative 
cooling 

- Distributed to reduce gradients from mirror 
heat losses 



Figure 3-1. Block Diagram for Vector Magnetometer. Electronics 
for one of the three channels of a two-loop system are shown. 
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Section 3 - System Design 


PHYSICAL AND INTERFACE CHARACTERISTICS FOR THE 
RECOMMENDED SENSOR DESIGN 

The recommended triaxial feedback sensor is located at the end of a boom 
in a thermal enclosure mounted on a thermally shielded base plate carry- 
ing the alignment mirrors. The electronics package is mounted in the 
spacecraft and communicates with the sensor by means of a boom cable. 

The weight and power requirements are less than specified. 

Here we briefly summarize the physical and interface characteristics 
of the design. The vector magnetometer sensor has been designed under 
the assumption that it will be carried on the same boom as a scaler mag- 
netometer sensor mounted not closer than 24 inches away (center-to- 
center ). 

The envelope, weight, and power consumption characteristics are sum- 
marized in the table opposite. In general, we have tried to achieve the 
minimum weight by making the sensor as small as possible consistent with 
meeting the modified noise specification. (The noise determines the size 
because bigger fluxgate sensors have lower noise, and the size of the sen- 
sors directly impacts the size of the feedback coils and surrounding struc- 
tural and thermal shield.) The power is kept as low as possible, and the 
solar heating has been reduced by use of a low-absorptive finish. This 
approach not only reduces gradients but minimizes the difference between 
the solar insolation and eclipse conditions. The heater size has been 
chosen to compensate for the losses through the superinsulations and from 
the surface of the alignment mirrors. 

Electronics - The electronics characteristics are also tabulated. It 
appears that the electronics will fit in a housing approximately 15 cm x 
13 cm x 8 cm, weighing 1. 02 kg. Four printed wiring assemblies are 
required, including the power supply and logic required for sequencing. 
These numbers are estimates since a detailed electronics packing design 
was not attempted. The estimates are based on component counts and 
should be quite close. 

Boom Cable - The recommended triaxial feedback scheme requires 
separate feedback and sensor coils, so that the number of wires per sen- 
sor axis is six, including the drive (three shielded pairs). In addition, a 
pair of heater sensor leads and a pair of power leads are required. A 
total of 22 wires, minimum, is required for the vector magnetometer. 

Command Interface - Commands will be required to place the instru- 
ment in a calibration mode and to provide self-test modes for the ADC. 
These commanded modes and their operation were not considered in the 
study. 


3-6 




- TIME -ZERO 

TABLE 3-3. OVERALL WEIGHT AND POWER 


Specified Actual 

Total Weight 

5 kg 3. 53 kg 

Total Power 

6 W 5. 2 W 

Sensor Characteristics 


• Envelope 

20 cm dia sphere 

• Weight 

1. 82 kg 

• Heater Power 

165 mW 

• Sensor Power 

35 mW 

• Thermal Finish 

- White on sun side 

- Gold or white on space side 

• Thermal Shield 

- Gold kapton except over alignment 
mirrors — 200 mW peak 

Electronics Characteristics 


• Envelope 

15 cm x 13 cm x 8 cm 

• Weight 

1. 02 kg 

• Power (total instrument) 

5 W 

Boom Cabling Requirements 


• Total number of wires 

22 

• Wire size 

30 AWG shielded 



c 
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Section 4 - Sensor Design 
Subsection 4. 1 - Configuration Analysis 

SUMMARY OF OPTIMUM SENSOR CONFIGURATION ANALYSIS 

To meet the alignment stability requirement in an optimum way, the sen- 
sor should employ triaxial vector feedback. Helmholtz coil pairs produce 
adequate field uniformity, especially if the coil spacing is slightly modi- 
fied. A concentric larger set of Helmholtz coils will provide almost com- 
plete cancellation of the external field. 


Before plunging into the details of a specific sensor design and per- 
formance analysis, we will present the analysis which led to the selection 
of the design. In the following pages, the rationale for a design optimized 
for alignment stability is fully developed; here, we summarize the analy- 
sis briefly. 

Selection of Feedback Configuration - To provide the required stability 
and linearity, the fluxgate sensors used must be used in the null mode. 

The earth f s field is cancelled by a feedback field produced by coils whose 
coil constant and alignment must be stable. There are at least two ways 
to organize the feedback. Each orthogonal fluxgate sensor can have its 
own feedback coil — a single-axis system — or three fluxgates can be 
located inside three orthogonal feedback coils which provide a vector feed- 
back field. We call this system triaxial feedback. 

The analysis presented in this subsection shows that the sensitivity to 
misalignment of the fluxgate sensors is far less if triaxial feedback is 
used, not only because the feedback loop is closed in all three dimensions 
for each sensor, but also because of a simpler mechanical configuration. 
Because alignment stability is the most critical requirement for the 
MAGSAT Vector Magnetometer, the triaxial feedback system was selected. 

Fluxgate Sensor Configuration - The noise of a fluxgate sensor depends 
on its size and shape; the sensors can be designed to attain a specified 
noise level. Accordingly, we have selected a sensor size which is com- 
patible with a modified noise level requirement of ly peak, selecting the 
diameter -to -height ratio to give a minimum package size for the three 
orthogonal sensors mounted at the center of the triaxial feedback field. 

Feedback Coil Configuration - The three sensors cannot be all at 
exactly the same place in the feedback field. The field therefore must 
meet certain uniformity requirements so that errors do not arise from 
instabilities of the sensor position or orientation. 

Attainment of a sufficiently uniform field is simply a matter of using 
an appropriate number of feedback coils. The problem is to attain ade- 
quate uniformity with coils of minimum size, weight and complexity. In 
principle, the use of many coils allows the use of very small coils; but, 
in practice, the mechanical design to avoid interferences between coils 
prevents the coil system from being very small. Analysis shows that 
adequate field uniformity is attained with a single pair of coils per axis, 
with a coil spacing slightly less than their radius (modified Helmholtz 
configuration). 

Cancellation of External Field - Because the leading term in the 
external field is a dipole field, the external field can be greatly reduced 
by cancelling the dipole moment of the main feedback coils with an outer 
pair of coils. Because Helmholtz coils have no quadrupole field, the 
outer coils should also be Helmholtz coils. Less power is required if 
the coils are large — at least twice the size of the inner coils. 
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TABLE 4-1. SENSOR CONFIGURATION ANALYSIS SUMMARY 


The optimum configuration of the 
to the rationale indicated. 


sensor was selected 


according 


Selection 


Rationale 


• Triaxial feedback 

• Specific fluxgate size 

• Helmholtz inner feedback coil pair 

• Large outer Helmholtz 
Cancellation coils 


- Low misalignment sensitivity. 

- Structural advantages. 

- To meet noise requirement. 

- Adequate uniformity. 

- No external quadrupole moment. 

- Cancel dipole moment of inner 

coils. 

- Low power. 

- No external quadrupole moment. 


• Modified Spacing of Helmholtz 
coils 


Make field flat at sensor* 
Minimize required coil size. 


c 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

DEPENDENCE OF SENSOR AND FEEDBACK CONFIGURATIONS ON 
MISALIGNMENT SENSITIVITY 

To select the most promising sensor configuration for the Vector Mag- 
netometer, it is necessary to use one of the three fluxgate geometries — 
parallel-rod, racetrack, or ring-core — in a feedback coil system which 
assures low sensitivity to sensor misalignment. 

Fluxgate Sensor Characteristics - The fluxgate sensor is a remark- 
able device, capable of measurement of magnetic fields over a very wide 
dynamic range. Fluxgates have been flown to the moon and to all of the 
inner planets and Jupiter, and a fluxgate instrument is now en route to 
Saturn on Pioneer 11. Others will fly on the MJS (Mars - Jupiter -Saturn) 
spacecraft. 

The fluxgate is basically a transformer with a square loop core so that 
it acts as a magnetic amplifier. The flux is made to saturate alternately 
in each direction by an alternating current in an ac drive winding. In the 
presence of an external field, a flux change at twice the drive frequency is 
generated in the core. This flux change can be detected by a sense winding 
and amplified to obtain a measure of the applied field. The sense winding 
is wound differentially to reject the drive; in practice, this means that the 
sense winding is simply wrapped around the core. The noise is essentially 
Barkhausen noise, while the offset appears to be produced by magnetostric- 
tion. The two effects are related: materials with low magnetostriction 

have low noise as well. The theory of the signal generation, noise, and 
offset is discussed in detail in Appendix A. 

Several different core configurations are possible — a pair of straight 
rods, an elongated (racetrack) toroid, or a cylindrical ring core. His- 
torically, the racetrack sensor evolved from the parallel rod sensor, 
while the ring core evolved from the magnetic amplifier, being simpler to 
manufacture because of its cylindrical shape. The principle of operation 
is the same for all three. Selection of the fluxgate geometry depends on 
the characteristics of the feedback coil system. 

Feedback Characteristics - The important aspect of the sensor design 
approach is not so much the selection of the type of fluxgate to be used, but 
the way the fluxgate is used in an overall feedback loop. Feedback must 
be used to generate a magnetic field which nulls the ambient field to attain 
the required linearity and insensitivity to misalignment error. In this 
manner the sensor is used only as a null detector, so that the stability of 
its gain and linearity are not important. Such feedback is commonly used 
in flight magnetometers for precision measurements. 

To minimize the effects of sensor misalignment, it is important to 
organize this feedback in the proper way. Historically, vector magnetom- 
eters have often been made by placing three single-axis fluxgates — either 
parallel-rod, racetrack, or ring-core — at right angles, with a separate 
feedback coil for each sensor, as shown in Figure 4-1. The Pioneer 10/11 
Vector Helium Magnetometer, built for JPL by Time-Zero, is an excep- 
tion: since the helium-gas sensor is inherently a vector device, the feed- 

back was accomplished by a set of three Helmholtz coil pairs which pro- 
duced a vector feedback field (see Figure 4-2). Ring cores appear to be 
the best kind of fluxgate sensor for this triaxial configuration because they 
allow a more compact structure, making it easier to produce a feedback 
field uniform over the sensor volume. 
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Misalignment Sensitivity - Comparing Figures 4-1 and 4-2, one might 
expect the single-axis and triaxial-feedback arrangements to have different 
sensitivities to sensor and structural misalignment. For example, the single- 
axis sensors are mounted apart from one another on an extended structure to 
avoid sensor crosstalk, while in the triaxial geometry the sensors can be 
mounted in a compact assembly. Because of the stringency of the alignment 
requirement, the configuration which has the lowest sensitivity to misalign- 
ment errors must be selected for MAGSAT. The selection of the most prom- 
ising geometry thus requires an analysis of the misalignment sensitivity of 
each configuration. This analysis, which leads to the conclusion that the tri- 
axial configuration is far less sensitive to misalignment error, is reported 
in the following topics. 
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Figure 4-1. Single- Axis Vector Magnetometer 
Each fluxgate sensor is surrounded by an inde- 
pendent feedback coil which nulls one component 
of the ambient field. A vector field in an arbi- 
trary direction is not nulled. 



Figure 4-2. Triaxial -Feedback Vector Magnetometer. 
Three orthogonal fluxgate sensors are mounted within 
a triaxial coil structure which completely nulls the 
ambient field in the neighborhood of the sensors. Ring 
cores are preferred for symmetry. 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

FACTORS INFLUENCING MAGNETOMETER ALIGNMENT STABILITY 

The vector magnetometer is sensitive to misalignments of the optical 
mount and sensor structure, feedback coils, and the sensitive magnetic 
cores. 


To understand the alignment problem completely, it is necessary to 
understand how the instrument data must be reduced. The alignment as 
well as the calibration of the vector magnetometer is ultimately specified 

by giving the 9 components of a 3 x 3 matrix G~ 1 such that B = G~ A (N - N Q ), 

where N is a vector whose components are the digitized output of the three 

signal channels, and B is the measured field relative to the coordinates of 

the optical attitude determination system. The vector N Q is the offset. 

The scale factor matrix depends not only on the orthogonality and scale 
factors of the individual sensors but also on the stability of the rotation 
between the sensor axes and the axes of the optical determination system. 

Several alignment factors (formally described by rotation matrices) 
are involved in transferring between the coordinate systems. These fac- 
tors are identified in the accompanying table and are discussed below. 

Optical Reference Mount Alignment - The mirrors used by the optical 
reference system mount on a plane surface facing the spacecraft. This 
plane surface is part of the structure which also supports the sensor 
assembly. No matter what sensor arrangement is used, this structure 
must be rigid — for example, it must not deflect in the earth’s gravitational 
field as a result of its own weight or that of the sensors, thermal shield- 
ing, or mirrors. 

Sensor Structure Alignment - The structure on which the sensor feed- 
back coils themselves are mounted must be structurally and thermally 
stable, because this structure defines the relationship between the sensor 
axes and the axes defined by the mirrors and their mounting plate. For a 
single-axis system, the sensors are well- separated so that the fields of 
their independent feedback coils do not interact. For such an arrange- 
ment, the stability of the intervening structure is particularly critical. 

For the triaxial arrangement, the sensors are closer together, but the 
feedback coil system is larger, so that close attention is required to make 
the mechanical and thermal feedback coil structure rigid. The struc- 
tural rigidity requirement appears to be easier to meet for the highly 
symmetrical triaxial arrangement compared to the extended rod or plate 
structures used for single-axis magnetometers we know about. 

Feedback Coil Alignment - The orientation of the feedback coils is 
evidently critical for either sensor geometry. The coils, which are wound 
on forms and potted, must be designed to maintain the required alignment 
stability relative to the mounting structure. To achieve a stability of 
5 seconds of arc requires a mechanical stability of 25 microinches per 
inch. This is of the same order as the stability required of precision 
optical systems. It appears easier to define the coil orientation in the 
triaxial arrangement, because it is larger; if the structure is rigid, a 
small motion of the effective coil axis leads to a relatively smaller angu- 
lar misalignment than does the same linear motion of the relatively small 
feedback coil used in the single-axis sensors. 
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TABLE 4-2. ALIGNMENT STABILITY FACTORS 

• Optical Reference Mount 

- must be rigid 

• Sensor Structure 

must be structurally and thermally stable 

• Feedback Coil 

- must have mechanical stability of 25 microinches 
per inch 

• Magnetic Core 

must be relatively insensitive to core movement 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

CALCULATION OF THE MISALIGNMENT SENSITIVITY OF A SINGLE- 
AXIS MAGNETOMETER 


The single-axis fluxgate sensor is sensitive to misalignment of both the 
magnetic core and the feedback coil. . 

One way to construct a vector magnetometer is to use three separate 
single-axis sensors at right angles. Consideration of this approach shows 
that the overall misalignment sensitivity is not determined only by the 
misalignment sensitivity of the feedback coil, as one might hope. The 
alignment of the fluxgate core itself is also an important factor and is, in 
fact, extremely critical. 

Simplified Argument - To see that there is a problem with the single - 
axis sensor, consider the following simple example. Suppose a single- 
axis fluxgate is located inside a feedback coil. Suppose that the feedback 
coil produces a field precisely along the nominal magnetic axis, so that 
the feedback coil is not misaligned. Now suppose that the magnetic core 
itself is misaligned by an angle A in the x-y plane, and suppose that a 
field By is applied along the y axis, normal to the nominal magnetic axis. 
This field will have a component B y sin A along the magnetic axis. To 
null the field, the feedback loop will force the feedback coil to produce a 
field of strength B sec A sin A. This field appears as an error in the 
measurement of the x-component of the field. Thus, the feedback loop 
does not correct for magnetic core misalignments. 

Complete Analysis - We can readily extend this analysis to include the 
effects of finite loop gain and feedback coil misalignment. The scalar 
field measured by the sensor (that is, the sensor output signal) is 

B g = Rg ' (B + B f ) U) 

where Rs is a unit vector along the actual magnetic axis, B is the ambient 

field to be measured, and Bp- is the feedback field. Offsets are neglected 
(not relevant to the argument). 

The feedback field is determined by the loop equation: 

B f - -R f M B s (2) 


The vector Rp is a unit vector in the direction of the field produced by the 
feedback coil, while M is the loop gain. The telemetered field is a scalar 
variable for each sensor axis; it has the value 

B T = M Bg (3) 

and is assumed to correspond to the component of the ambient field along 

the nominal direction specified by the unit vector Rjq. 

We can now compute the deviation between the telemetered field com- 
ponent and the actual nominal value, 


B 


N 


= R 


N 


B 


(4) 
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Eliminate Bp from Equation (1) by use of Equation (2): 

B s - R s • (B - R F M Bg) (5) 

This leads to the result 


R, 


B 


B S = 


1 + M(R g • R f ) 

The error in the telemetered component is then 

B 


AB, 


[MR g - M(R g • 

1 + M (R. 


V 

rV 


( 6 ) 


(7) 


This is the complete formula, showing the effect of finite loop gain and both 
sensor and feedback coil misalignment. 

_ , -F ^ 

In the case in which Rp = Rj^ (feedback coil aligned), while R is not 
parallel to the nominal axis Rjvj* there is an error proportional to v 

[R s /(Rs • Rn)] • B - ^We can plug in the assumptions of the simple example 
above, by assuming B = B y , R y , Rg = R x cos A + R y sin A, R N = R*; the 

resulting error for infinite loop gain is, as before, B y sin A sec A. 

Example - The accompanying figure shows the arrangement of a typical 
fluxgate sensor (ring core or racetrack — the problem is the same). The 
permalloy tape is mounted on the bobbin without excessive constraints to 
avoid stress-induced offsets. A motion of this tape by only 5 seconds of arc, 
or 25 microinches per inch, is enough to equal the error allowed for MAGSAT. 
This amounts to 0. 12 pm/inch or 1200 A — an optical tolerance exceedingly 
difficult to meet in an electrical device which is essentially a transformer 
whose core is concealed by windings. 

The single-axis magnetometer thus does not benefit very much from the 
use of feedback, as far as core misalignment errors are concerned (it bene- 
fits greatly, of course, in linearity and stability). We can correct for these 
errors by using the sensors in a manner which closes the loop in all dimen- 
sions so that the sensor acts as a true null sensor. 
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Figure 4-3. Typical Fluxgate SenBor. A motion of only 5 seconds of 
®^c (25 microradians) of tape, bobbin, or feedback coil is sufficient 
to produce all the error allowable. Rotations about the sensitive axis 
are also equally critical for either the ring-core or racetrack sensor. 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

USE OF TRIAXIAL FEEDBACK TO REDUCE SENSOR MISALIGNMENT 
SENSITIVITY 

By using a triaxial sensor within a triaxial feedback coil, the sensitivity to 
sensor misalignments is greatly reduced, so that the triaxial geometry is 
the best choice for the MAGSAT Vector Magnetometer, 

The basic problem with the single-axis sensor is the fact that it has no 
feedback for field components perpendicular to its axis. It seems intui- 
tively obvious that this difficulty can be corrected by use of a three dimen- 
sional feedback field, using three orthogonal sensors. This supposition 
is readily confirmed by analysis. 

Consider now a vector equation for the three sensor outputs 

b s = r s (b + b f ) (1) 

As before, B is the ambient field and B-p is the feedback field, but Rg is 
now a 3 x 3 matrix (an orthogonal rotation matrix if all three sensor scale 
factors are equal). 

As before, the feedback field is determined by an equation of the form 

B f = -R f M B g (2) 

except that Rjr and M are now 3x3 matrices: Rp is an orthogonal rotation 
matrix, corresponding to the misalignment of the feedback coils, while M 
is a diagonal matrix, whose diagonal components are the three respective 
loop gains, which of course may differ somewhat. The telemetered field 
is 

B T = M B s = -R f _1 B f (3) 

while the error vector is AB = Bt - B. 

This equation is solved by eliminating Bg and Bj, to obtain the result 

B t = (1 + M R s R f ) _1 M Rg B (4) 

The error in the telemetered field is then given by 

AB = (1 + M RgR F ) _1 [M Rg(l - R f ) - 1] B (5) 

For infinite loop gain along all three axes, the result is 

AB = R f ' 1 ( 1 - R f ) B = (R F 1 - 1) B (6) 

This error is proportional to the differences between the feedback coil rota- 
tion matrix and the unit matrix, which means that misalignments of the 
feedback coil by an angle A produce error signals in orthogonal channels 
of order B sin A. 

The error is independent of sensor alignment (but not sensor offset 
and noise, of course). This arrangement then operates like a true feed- 
back system in that only the properties of the feedback element are sig- 
nificant in determining the transfer function of the system. 
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E xample - Suppose that one allows for reasonable, easily-inspected 
tolerances for the fluxgate sensor shown in the previous discussion. One can 
readily control the tape angular alignment to 1 milliradian or 200 seconds of 

Th C a . t ° le u rance /° j t L imes less stringent than that for the single-axis sensor 
Thus it takes a feedback factor, M, of 40 or more to bring the effect of a 

Wk V C ? 8e mi8 *? l * nm « nt within specification. However, a much larger feed- 

' fx 1t 0rder ° f 20 ' 000 ’ is "«<ied ‘o linearity and gain 

of the 1 If' , ha , result of usin 8 tr iaxial feedback is to make alignment 
of the fluxgate itself totally uncritical. g 

Summary - The analysis just presented has shown that: 

• There is a serious problem in making stable fluxgate sensors for 
single axis system. 

# J„ hl f probl ® m can be eliminated by using triaxial feedback - that is, by 

enclosing the sensors in a three-dimensional feedback field. Y 

is s ”n U cHt 1 !?kl th b u ft i h, n . ment M nd Stabll,t1 ' of lhe ! ee( i! j ack coil structure itself 
same problem. Pnoblem .. not new - the single-axis system has the 

sho l: 1 .^. i The°MA, n sAi C v nfi , fi " r ?: i0n ‘ The analyais thal has been Presented 
f Vector Magnetometer should employ triaxial feedback 

to reduce the sensitivity to fluxgate sensor misalignment: There are some 

minor new problems resulting from the fact that the sensors are not all at the 

stringed Bu? as' ^ fV" 1 '""’? requirements are somewhat more 
readny solvedl’ " ‘ n th " n ' X ‘ “ >pic ’ ,he8e Problems are 

TABLE 4-3. HOW^TRIAXLAL FEEDBACK REDUCES MISALIGNMENT 

The mathematics in the tex t expresses the following physical facts: 

• Each sensor responds to fields along its own axis. ~ 

• Each sensor produces a feedback field to null that axial field. 

• This nu n s th e total (vector) field, even if sensors are not orthogonal. 

• Since the total field is zero, its alignment is irrelevant. 

• approximation* Se " S ° r mi8ali ® nm ‘ at ia irrelevant to the first order of 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

UNIFORMITY REQUIREMENTS FOR THE TRIAXIAL FEEDBACK FIELD 

Non-uniformity of the field and the physical separation of the sensors lead 
to a gain correction which must be held stable. For variations in sensor 
position and orientation of 0.001 inch and 0.001 radian, the feedback field 
must be flat to about 5 parts in 10 3 to reduce errors to the 0. 5 Y level. 

In actual practice, the three sensors enclosed in the triaxial feedback 
volume have finite dimension, so that the places at which they measure the 
field are different. If the feedback field is not uniform, the field meas- 
ured by the sensors will not be the same. The sketch opposite shows the 
geometry. The three fluxgate ring- core sensors are mounted orthogonally 
inside the feedback coils. The arrows shqlw the sensitive directions. The 
z-axis sensor is offset in the x direction, and so on in cyclic order. 

We can see how this complication affects the analysis in the previous 
topic by adding a term in Equation (1) to represent an offset in the feed- 
back field. Equation (1) of the previous topic now reads 

S s * R S (6 + S F + R 0 S F> (1A) 

The added term involving the matrix Rq accounts for the difference between 
the feedback field at the center of the feedback coils and the average field 
sensed by each sensor. 

Following through the analysis as before, we now obtain the following 
equation in place of Equation (6) of the previous topic 

AB = (R F 1 - Rq' 1 ) B (6A) 

There is now a correction to the output proportional to the field, so that 
the added term - Ro" 1b represents a gain or alignment variation from the 
ideal situation. 

Translational Error - Instability of the sensor position in a non- 
uniform feedback field leads to an error which appears as a variation in 
the coil constant and hence the closed-loop scale factor. In keeping with a 
total error budget of 1. 5Y, we assign 0. 5Y to this error. The instability 
of the feedback field must be less than 

^ = 0. 5Y/64, 000 Y = 7. 8 ppm (7) 

Referring to the accompanying figure, we see that the critical direc- 
tion of motion for each sensor is along the direction of the offset (x-direc- 
tion for the z-axis sensor). Motion along the other axes is far less criti- 
cal by symmetry (the field has an extremum). We expect to be able to hold 
tolerances of 0.001 inch worst-case; therefore, the field should have a 
gradient of at most 8 ppm/mil. 

Orientation Error - A field along the x axis is cancelled by feedback 
action at the location of the x sensor. Because of field non-uniformity 
there is a small residual field AB X at the y and z sensors; the error in the 
measured z-axis field produced by a misalignment of the z-axis sensor by 
an angle A 0 is then 

AB = AB A6 (8) 

z x 
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Hft« 8 ^A\ Wh T h iS ’ ° f C ° UrSe ’ a s P e cial case of the general result in Equa- 
tion (6A). To ensure an overall error less than 5 seconds of arc, we assign 

fie"[d r error S of 0 r 8 O v r °‘ *T t °! his Source ' corresponding to a 

field error of 0. 8Y, or 12 ppm cpmpared to 64,000Y. It is easy to hold the 

f r orientation staWe to 1 milliradian (200 seconds of arc), which requires 
that the field should vary by less than 1. 2% between sensors, or about 770Y 
This analysis is very conservative. Precision single-axis sensors such 
J S , h ° Se used on MJS (in Project MAGNET) and in earthquake prediction 
exhibit axis stabilities in the 10-50 arc-second range, a number not good 
qHqqq, ° r , ^ consistent with dimensional stability in the 

. 001 -inch range. We therefore expect that the actual errors will be 
maiier than allowed when a field meeting the above uniformity requirements 

gate sensors used® The*’/'. feed . ba, ; k coi ' s compatible with the specific flux- 
considered ne« determination of the fluxgate sensor configuration is 


TABLE 4-4 . 


REQUIREMENTS FOR FIELD UNIFORMITY 


Change in field resulting from 0. 001-inch sensor motion < 0. 78V 
Difference in axial field between sensors < 770 Y 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

DEPENDENCE OF NOISE AND OFFSET ON FLUXGATE DIMENSIONS 

The noise and offset of the fluxgate sensors depend on the sensor core 
dimensions. Larger sensors exhibit better offset and noise performance. 

Now that the triaxial-feedback system has been selected, we need to 
establish the characteristics of the orthogonal fluxgate sensors and the 
corresponding feedback coils. 

To decide how large the uniform-field region must be, we must decide 
how large the fluxgate sensor array must be. We can do this by relating 
the physical size of the sensors to the required offset and noise perform- 
ance using the simplified theory presented in Appendix A, where the 
dependence of these parameters on the sensor dimensions is derived. The 
results of this theoretical analysis are summarized briefly here. 

Signal Flux and Capture Cross Section - The field inside the iron is 
related to the field outside by a formula involving the demagnetizing factor, 
which depends on the geometrical shape of the sensor core. While it is 
conventional to discuss the behavior of the fluxgate in terms of the demag- 
netizing factor, this factor is not the relevant figure of merit. What mat- 
ters is the "capture cross section" — the effective area which multiplies 
the external field to obtain the sensed flux. The voltage sensed by the 
pickup winding is proportional to the rate of change of flux in the core. 

This flux can be expressed as the product of the external (ambient) field 
and the capture cross section — the effective projected area of the core. 

For a long thin rod with high permeability, this cross section varies as 
the square of the rod length and is only very weakly dependent on rod 
diameter. For a ring-core of finite tape width, the capture cross section 
is approximately equal to the product of the sensor diameter and tape 
width, or to the projected area of the tape viewed along the sensitive axis. 
Thus, the second harmonic signal is* essentially determined by the core 
dimensions . 

However, it is not the signal level, but the signal-to-noise ratio — or 
noise-equivalent field — which matters. 

Noise-Equivalent Field - The noise in magnetic amplifiers appears to 
be produced by domain reversal (Barkhausen noise). The mean-square 
voltage noise on the second harmonic output is directly proportional to the 
area of the core and inversely proportional to the effective length of the 
magnetic circuit. Combining the noise with the sensitivity, one can derive 
an expression for the noise-equivalent field of a ring-core (Table 4-5). It 
shows that the noise -equivalent field is proportional to (t/wd3)l/2. where 
t is the core tape thickness, w is the tape width, and d is the diam- 
eter of the ring core. 

Offset-Equivalent Field - The offset appears to be produced by mag- 
netostriction, which can be reduced by a proper selection of the core 
alloy, since permalloys with 81. 3% nickel in iron have zero longitudinal 
magnetostriction. The estimate derived in Appendix A is shown in 
Table 4-6; the offset-equivalent field varies as t/d, the ratio of the tape 
thickness to tape diameter. 

Both of these estimates assume that the demagnetizing factor for a 
ring core of finite height can be calculated using an infinite- cylinder 
approximation. This approximation should be good as long as the core 
tape width is much greater than the thickness. 

Sensor Dimensional Tradeoff - Because the noise depends on both the 
diameter and width of the core, there is a dimensional tradeoff or 
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optimization involved when the sensors are required to fit within a specified 
volume. This optimization is considered in the next topic. 


TABLE 4-5. NOISE EQUIVALENT FIELD OF A RING-CORE FLUXGATE 


The noise appears to be generated by domain reversal. 
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AB^ = noise - equivalent field 

B = saturation field 
s 

n^ = number of elementary dipoles per domain 

v , = domain volume 
a 

f = drive frequency 
T = signal integration time 
t = tape thickness 
d = core diameter 
w = tape width 


TABLE 4-6. OFFSET EQUIVALENT FIELD 
The offset is generated by magnetostriction. 


AB q = offset-equivalent field 

X g = saturation magnetostriction 
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Section 4 - Sensor Design 
Subsection 4. 1 - Configuration Analysis 

OPTIMIZATION OF FLUXGATE GEOMETRY UNDER THE NOISE AND 
SIZE CONSTRAINTS 

To minimize the noise-equivalent field with a specified sensor offset, the 
ring-core tape width should be one-third of the tape diameter. With this 
shape, one- centimeter diameter sensors should be large enough to meet 
the noise specification. 

To apply the triaxial feedback concept, it is necessary to place three 
orthogonal fluxgate sensors at the center of a field-free region generated 
by action of the feedback coils. To simplify the feedback coil design and 
minimize the overall weight and volume, it is desirable to minimize the 
sensor size consistent with the noise and offset requirements. 

Fluxgate Mounting - The general arrangement of the three orthogonal 
sensors is shown in the sketch opposite. The field at the location of the 
sensor will depend mainly on the offset, R, between the center of the flux- 
gate tape and the center of the feedback field, which is at the center of a 
cube on which three orthogonal fluxgates are mounted. 

The radial offset, R, can be expressed in terms of the diameter and 
height of the core as well as the buildup in diameter, d Q , and height, w Q , 
resulting from the finite dimensions of the bobbin and sense winding. 

(The sense winding can be wound directly on the bobbin to minimize the 
height buildup w Q . ) The radial offset is given by 

R=4(d+d +w+w) (1) 

2 o o ' 1 

Optimization - The noise equivalent field is inversely proportional to 
d3/ZwT7T To minimize this noise subject to the constraint that the radial 
offset be constant, we maximize the quantity 

S 2 = d 3 w - X.R (2) 

where \ is a constraint parameter (Lagrange multiplier). 

By setting the derivatives of this expression with respect to d and w 
equal to zero, and eliminating the constraint parameter from the two 
equations which result, one obtains the optimum condition 


For optimum noise performance, the tape width should be one-third of the 
diameter. 

Determination of Core Size - Now that this geometric optimization has 
been performed, one can estimate the sensor size needed to meet the noise 
requirement. The MAGSAT study specification requires a 0. 1 zero-to- 
peak noise for a 10-second scan with a 25 Hz bandwidth. But, as dis- 
cussed earlier, it has been agreed that this specification is unnecessarily 
stringent. Because the digitizing resolution is 2Y, a peak-to-peak noise 
of 2 Y should be acceptable — a factor of 10 greater than the specified value. 

Now, we know that NOL-type one-inch ring cores, with a smaller 
width-to-diameter ratio, exhibit peak noise levels of better than 0. 1Y. 

Since the noise -equivalent field, by our optimizing constraint, varies as 
d2, we should be able to reduce the diameter of the sensor by a factor 
of /TO and still obtain adequate performance. A further safety margin 
is obtained by using the optimum width. Thus, a core whose diameter 
is 0. 3 inch or greater should provide more than adequate performance. 
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widtTof o e i?f ed i di J, meter of a ^ out 1 cm - actually 0. 375 inch, with a tape 
width of 0. 125 inch. We expect that the rms noise will be about 0. 1Y. 

rp, ° f // et P erf °rmance - Will this sensor exhibit adequate offset stability? 
The offset drift because the core is smaller will be larger than that of a 
one-inch NOL sensor. However, the typical sensor exhibits offset stability 

MAGSAT vlrf ° f M gam f ma over a ver Y la rge temperature range, while the 
^ Vector Magnetometer Sensor will be temperature-controlled. Off- 
set drift is thus not perceived as a problem. 

Now having selected the appropriate fluxgate size, we need to select the 
imp lest feedback coil configuration which provides the desired feedback field 

un.forrmty and yet provides low power and ease of assembly. Th” selection 
is discussed in the following topic. selection 
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Fluxgate Sensor Mounting. To minimize the height 
sense winding is wound directly on the bobbin with- 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

FINDING THE SIMPLEST ARRANGEMENT FOR GENERATING A 
UNIFORM FEEDBACK FIELD 

Because of requirements for structural stability and ease of assembly, 
the simplest feedback coil configuration which provides acceptable uni- 
formity should be used. A system of Helmholtz coils appears best. 

Theoretically a perfectly uniform field can be generated by spacing a 
large number of windings uniformly along the polar axis of a sphere (or, 
in fact, a spheroid, either prolate or oblate). Good approximations can 
be produced by using a finite number of windings . . . the more the better, 
from the standpoint of theoretical uniformity. Yet, there are other 
requirements. Overkilling the uniformity problem can result in a diffi- 
cult assembly, with a great many critical windings whose relative align- 
ment mO.st be stable. Then, too, one can show that the uniformly wound 
sphere requires more drive power than a simpler configuration. A dis- 
proportionate amount of coil power is required to assure the ultimate 
in field uniformity. Minimizing coil power is important to reduce mis- 
alignments caused by self- heating, and concomitant thermal expansion of 
the feedback coils. 

Increases in uniformity are attained by increasing the number of wind- 
ings, spacing them nearly uniformly along the polar axis. The simplest 
configuration is of course a single coil; the next simplest is a pair of coils 
separated by their radius (Helmholtz coils). As is shown in elementary 
textbooks, even this simple configuration already provides a fairly remark- 
able degree of uniformity, because the first three derivatives of the axial 
field vanish at the center of the coil system. 

The problem with adding more coils is that coils of increasingly 
smaller diameter are needed near the poles of its sphere. Yet, access 
along the polar axis is needed to insert the orthogonal fluxgate assembly 
and to provide a mounting flange of adequate diameter. (See Figure 4. 6. ) 

It is also important to use the minimum number of coils, in view of the 
need for a simple, rigid structure with a larger number of coils. If many 
coils are used, the mechanical interferences between coils becomes a 
real problem; (four coils interfere in 48 places), so that a larger struc- 
ture must be used to handle the interference problem. One might as well 
use a simpler configuration of larger diameter to get the same degree of 
uniformity. 

Thus, although the use of four or more coils can theoretically produce 
a field of spectacular uniformity, mechanical considerations motivate the 
search for the simplest acceptable coil configuration. The analysis pre- 
sented subsequently shows that one pair of feedback coils is sufficient to 
meet the uniformity requirements. 

This Helmholtz coil configuration is sketched opposite. Each pair of 
coils is separated by the mutual radius. While the sketch shows that the 
coils have the same diameter, the coils for different axes must in fact 
have slightly different diameters to avoid the mechanical interference 
between the coils. This complication is neglected in the conceptual dis- 
cussion which follows but is addressed later. 

Another complication is introduced by the need to have a small exter- 
nal field at the scalar magnetometer (0. 5 at 24 inches). If the proposed 
coil system does not provide the required far-field characteristics it will 
be necessary to add cancellation coils. The best cancellation system 
appears to be another set of larger Helmholtz coils, as shown in the Figure. 
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avmg initially selected this design for investigation (for the reasons 
summarized in the list below) we need to determine the uniformity of the 
leld to see if such a design will produce a sufficiently uniform field with 
co, s of ac.eptable .l«, weigh,, end power. As we /hall see the Helmholtz 
*! d S' 1 appears to proride the required performance, especially if the 

of the o P ffs“ 1 /e„s e or JUS ' SUgh ‘ ly maka " eld « th * ^cation 


TABLE 4-7. 


REASONS 

POSSIBLE 


FOR DESIRING THE SIMPLEST 
CONFIGURATION 


• Maximize access to sensor volume. 

• Provide strongest possible mounting flange. 

• Provide simplest possible structure. 

• Reduce coil power, 

• Reduce number of critical tolerances. 

• Reduce number of mechanical interferences. 



^- Complexity and Access - If many coils are used 
(the dashed lines show the location of a four-coil array) the 
mechanical complexity increases and the size of the access 
port along the coil axis is reduced. access 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

MAGNETIC FIELD OF COIL PAIRS 

For purposes of analysis, the magnetic field of a pair of coils can be cal- 
culated from a spherical harmonic expansion. If the spacing of the coil is 
equal to the coil radius, the field between the coils becomes quite uniform. 

Having tentatively selected the two-coil Helmholtz configuration for the 
feedback coils, we must determine the coil size which produces accept- 
able field uniformity according to the sensor position stability require- 
ments already established. For analytical purposes, the most convenient 
form for the field seems to be a spherical harmonic expansion. By using 
this expansion, one can readily estimate the field for coils of small cross 
section and select the arrangement most likely to work. The selection 
can then be verified for an actual coil of finite cross section by computing 
the field directly by numerical quadrature of the integral for the vector 
field; this procedure is described later. 

Spherical Harmonic Expansion for Coil Fields - To keep things simple, 
we first give the formulas for a single pair of coils. The complication 
resulting from the addition of cancellation coils will be treated later. 

The expressions given opposite are obtained from the single-coil field 
given by W. R. Smythe (Static and Dynamic Electricity , 3rd edition, 

McGraw Hill).* All formulas assume the use of SI units, although conven- 
tional English units have been used to describe mechanical dimensions. 

The formulas give the radial and polar components of the field, B r and 
B e in terms of the field point as shown in Figure 4-7. The coils are sym- 
metrically placed along the polar axis of a sphere of radius a at polar 
angles a and tt - a , as shown. 

The first term in the internal field expansion represents a uniform 
field along the polar axis (z-direction), which has the radial and polar 
components By = B 0 cos 0, Bq = Bq sin0. Externally, this term pro- 
duces a dipole field. The axial field at the center of the coil pair is: 

B 0 = P Q Nl/a sin 2 a (see footnote), while the external diple moment is: 
m = irNIa^ sin^a. 

Helmholtz Coils - If we select the coil spacing so that the second term 
in the expansion vanishes, we have the so-called Helmholtz configuration. 

The vanishing of the second term requires that P 3 (coso) = 1/2(5 cos 2 o- 1) sinCX 
= 0, or that cos & = 1//5. The radius of the coils is then sin a = a/ 4/5; the 
spacing is equal to the radius. 

Since the second term in both the internal and external field expansions 
vanishes, the Helmholtz coil configuration has an internal field flat to 
order r 4 at the center of the field, and also no external quadrupole (r~ 5 ) 
field. This fact suggests that cancellation of the dipole moment of an inner 
pair of coils with an outer pair which are also Helmholtz coils will lead to 
a very rapid falloff of the external field (as r -7 ). Thus, larger coils may 
be used to obtain internal field uniformity and — at the same time — a low 
external field at the nearby scalar magnetometer. 


*These formulas can be used to show directly that the field inside a 
sphere wound uniformly along the polar axis is uniform. The current den- 
sity of uniformly spaced coils has only one term in a spherical harmonic 
expansion, proportional to P 0 (cos a ) = 1. Integrating over this current 
density leads to a uniform field because all the higher-order turns in the 
expression to the right drop out because they are orthogonal to P Q (cos a). 
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• Sha11 now see , how this works and then address the problem of select- 

mg the minimum coil size. ieci 


TABLE 4-8. SPHERICAL HARMONIC EXPANSION FOR FIELD OF A 
COIL PAIR 

The field is given in spherical polar coordinates for N turn coils of 
infinitesimal cross section wound in the azimuthal direction on a 

afona / adlUS a and carrying a current I. The position of the coils 
ong the axis is specified by an angle a (see Figure). The functions 

P n (x) are associated Legendre functions. The value of p is 
4tr x 10 H/m. ° 


Field Inside Sphere 
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Figure 4-7. Symmetrical Coil Pair. 
The coils are symmetrically located 
along the polar axis of a sphere of 
radius a. 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

CANCELLATION OF THE EXTERNAL FIELD 


Cancelling the dipole field by use of an external pair of coils with twice the 
radius of the inner pair leads to a very low external field, less than 0. 18Y 
at a distance of 24 inches for eight-inch diameter outer coils, or 0. 007 Y 
for 5-inch coils. 

As presently conceived, the MAGSAT has a single boom, which sup- 
ports a magnetometer mounting shelf carrying a scalar magnetometer as 
well as the vector magnetometer described in the present study. The 
external field from the vector magnetometer should not disturb the scalar 
magnetometer which may be placed close to the vector instrument. Vari- 
ous specifications for the external field have been discussed: we settled 
on a 0. 5Y requirement at a distance of 24 inches, or 60 cm from the cen- 
ter of the vector magnetometer. 

Since the main term in the external field is a dipole term, the obvious 
thjng to do is to provide external coils which cancel the dipole moment of 
the inner coils. If the inner and outer coils have radii Rj and R 2 and 
ampere-turns Njlj and Njl2* the dipole moment cancellation requirement 
is 

= »n 2 i 2 r 2 2 (1) 

The coils are wound so that the fields oppose one another. It is desirable 
that the coil pairs be connected in series, so that 1^ = I 2 ; then. 



The field inside the coils is inversely proportional to the radius. If both 
coils are Helmholtz, 



Thus, even if the ratio of the coil radius R 2 /R 1 is only 2, the internal 
field is only reduced by one-eighth, a reduction easily compensated for by 
increasing the drive current by the factor 8/7. 

Suppose that the external coils have the maximum diameter permitted 
by the study specification (about 8 inches or 20 cm, allowing something 
for thermal shielding). Assuming that the inner coils have half that diam- 
eter, we can estimate the external field from the next term in the external 
field expansion. The dipole and quadrupole field are both zero; the octupole 
field, calculated from the formulas given opposite is only 0. 18Y, much 
less than the 0. 5 Y requirement. The external field decreases very rapidly 
if the coils are made smaller (as the seventh power of the diameter). For 
the five -inch outer coils actually recommended, the external field at 
24 inches is only 0.007Y. 

Notice that the dipole moments must cancel to within about 1 percent, 
since the dipole field of the inner coil is 51. 8Yand the field is to be reduced 
to 0. 5Y. That level of cancellation accuracy is consistent with the required 
accuracy with which the absolute value of the coil constant must be 
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( controlled. While this cancellation accuracy requirement can be met, it 

would be easier if the coils and the external field were smaller. Thus, for 
this reason as well as to minimize the weight and thermal gradients, we wish 
to make the coils as small as possible consistent with the required uniformity 
of the inner field. We look into this next. 


TABLE 4-9. STRENGTH OF EXTERNAL FIELD FOR 
HELMHOLTZ COILS 

The radial component of the external field along the 
polar axis is given in terms of the internal field at the 
center of the coils. The expansion is carried out to 
order r~ 7 (octupole term). Numerical values apply 
for four-inch diameter inner coils, eight-inch diameter 
outer coils, and a 64000 V central field. 


Field of a Single Helmholtz Coil Pair 

3 



Residual Field of Two Helmholtz Coil Pairs of 
Opposite Dipole Moment 
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Numerical Values for Field at 24 Inches 

• Dipole field of uncorrected inner coil 51.8V 

• Total residual field after correction 0. 18 Y 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

MODIFYING THE FEEDBACK COIL SPACING TO MINIMIZE COIL SIZE 

The uniformity of the field at the offset fluxgate sensor can be improved 
by decreasing the coil spacing, making the axial field flatter at the loca- 
tion of the sensor. 


All of the conventional arrangements of coils used to produce uniform 
fields, including the Helmholtz arrangement, produce the greatest uni- 
formity at the center of the coil system. However, in the triaxial feed- 
back arrangement the sensors must be offset. It makes sense to modify 
the coil arrangement to flatten the field at the location of the sensor rather 
than at the center of the coils. Then smaller coils can be used and a 
lower weight can be attained. 

As noted previously, the critical direction for translation is in the 
direction of the existing sensor offset. In the other directions the sensor 
is symmetrically located, and translational motions give rise only to very 
small second-order changes in the average field. Single-axis sensors 
symmetrically located in fields of much poorer uniformity provide ade- 
quate performance with respect to sensor translation, so we expect no 
problems with motion along the directions in which the field is symmetric. 
Thus, we wish to make the field more or less symmetric in the offset 
direction as well. 

The actual way in which the sensor responds to the details of the field 
is very complicated. To solve the problem in detail, it would be necessary 
to obtain a numerical solution of the magnetic boundary value problem in 
a non-uniform field — a very respectable exercise in numerical analysis, 
even if one assumes constant permeability in the sensor. A more prac- 
tical approach is to adjust the field for a maximum at the center of the 
sensor for a first approximation and then to optimize the coil separation 
experimentally to yield zero first-order tolerance to sensor position. The 
design to be presented has so much margin such that further optimization is 
probably unnecessary. 

We wish first to establish the spacing which provides a suitable 
extremum in the field at the offset position of the sensor. If the variation 
of the field averaged over the sensor is then small enough so that sensor 
motion gives rise to an acceptably small change in the average field, we 
can expect that the actual performance will also be acceptable. 

Field Calculation - To evaluate the variation of the field as a function 
of coil spacing, it is more convenient to solve for the field directly by 
numerical integration than to use the spherical harmonic expansion. The 
axial and radial field components of the field for a single coil are given in 
cylindrical coordinates (p, z) by the integrals opposite (Table 4-11) for a 
coil of radius a and axial distance z Q . These integrals can be shown to be 
expressible in terms of elliptic integrals, which have known numerical 
approximations. For practical purposes it is easier, if somewhat slower, 
to do the integrals by numerical quadrature. We used Simpson’s rule to 
obtain the plot opposite (Figure 4-8), which shows the radial dependence 
of the axial field on the central plane for various inner coil spacings with 
two outer cancellation coils of twice the radius and separation. 

From these curves, we can estimate the sensitivity to sensor transla- 
tion and select the diameter of the inner coils. 

To make the coil size selection, we selected one of the curves in the 
plot, located the sensor on the peak of the field curve, and calculated the 
corresponding coil radius and location of the sensor tape edges. The 
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slope of the field at the edge of the sensor is taken as a very conservative 
measure of the sensitivity to sensor translation. 

Following this procedure, one finds that if 2. 5-inch diameter coils are 
wound on a sphere, spaced by 93. 5 percent of their radius, the field error 
per mil of translational sensor motion will be less than 0. 7Y. This design 
is adopted as the study baseline. 

Since we have pushed the inner coils together so that the configuration is 
no longer exactly Helmholtz, the quadrupole moment of the inner coils is 
now non-zero. In principle, the outer coils can be slightly pushed together 
also, but by a smaller amount, to compensate for the added field. However, 
a calculation shows that the added field at 24 inches is only 0.0013Y, so that' 
one might as well leave the outer coils in the Helmholtz configuration. 

TABLE 4-10. GENERAL EXPRESSION FOR THE FIELD 
PRODUCED BY A CURREN T 

The following integral must be performed for each current- 
carrying element of length ds at a distance r from the field 
output: 


2^/jdsx, 


TABLE 4-11. COMPONENTS OF THE FIELD OF A CIRCULAR COIL 

The radial and axial components, B p and B z , are given in terms of 
the coil radius a axial position z D at the field point (p, z). 


p NI 
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(a + p + (z - z ) - 2ap cos <p ) 


a (a - p cos <p ) 

(a + p + (z - Z Q )^ - 2ap cos <p ^ 


372 d< ? 


d <p 


Figure 4-8. Variation in Axial 
Field for Various Radial Off- 
sets in Central Plane. The 
field variations apply for a 
64000 y central field. For the 
sensor position indicated the 
field can be flattened by using 
2. 5-inch diameter inner coils 
separated by 93. 5 percent of 
the radius so that variation of 
the field resulting from sensor 
radial motion is negligible. 
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Section 4 - Sensor Design 

Subsection 4. 1 - Configuration Analysis 

DERIVATION OF STRUCTURAL, STABILITY REQUIREMENTS 

The requirements for structural stability follow from the field equations. 
The most critical parameter is the coil alignment. The dimensional sta- 
bility required to ensure scale factor stability is equally critical. 

To do the structural and thermal design, it is necessary to know the 
sensitivity of the field to the coil parameters. These sensitivities can be 
readily calculated from the field equations or from geometrical 
considerations. 

Coil Alignment Stability - Misalignment of the magnetometer can arise 
from structural or thermal instability of the feedback coil orientations 
with respect to each other or with respect to the alignment surfaces; the 
total angular instability from either of these sources is not to exceed 
5 seconds of arc. This requirement is very stringent, since it corre- 
sponds to a dimensional change of 25 ppm (25 microinches per inch). Note 
that if one coil of a pair tilts with respect to the other, the net misalign- 
ment is half the tilt. Pincushion or other axially symmetric distortion of 
the coil will cause a scale factor error but not a misalignment. 

Scale Factor Stability - The required scale factor stability is equally 
difficult to attain. On the opposite page (Table 4-12), the central field has 
been expressed in terms of the diameters and spacings of the inner and 
outer cancellation coil pairs. By differentiating this formula, one obtains 
the sensitivity of the central field to changes in the dimensions. For the 
nominal design in which both the inner and outer cancellation coils are 
Helmholtz coil pairs, with the outer coils having twice the diameter of the 
inner coils, the numerical values shown are obtained. The most sensitive 
parameter is the inner coil spacing, which produces a fractional field vari- 
ation which is 24/35 the fractional variation in the spacing. The absolute 
sensitivity is 35Y/mil for a 1. 25-inch spacing. A tolerance of 1. 5Y was 
assigned to the feedback coil constants in the error allocation presented 
earlier, or 23 ppm in terms of the 64000Y full-scale field. Assuming that 
the structure expands uniformly, the overall sensitivity to thermal expan- 
sion can be calculated by adding the coefficients opposite algebraically. 

The result is unity (35/35) so that the regional structural dimensional sta- 
bility is 23 ppm — about equal to the requirement for alignment stability. 

Fluxgate Translation and Orientation - Errors from these sources are 
controlled by triaxial feedback with proper shaping of the feedback field, 
as explained in detail previously. The field has been designed to produce 
acceptable scale factor stability and alignment if the fluxgate position is 
stable to 0.001 inch and the orientation is stable to 1 milliradian. These 
requirements are summarized in Table 4-13. Only the overall require- 
ments are given in the table; the mechanical design must consider other 
details, since the coil system may very well expand non-uniformly (distort). 
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TABLE 4- 12. 


SENSITIVITY OF THE FIELD TO DIMENSIONAL VARIATIONS 
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• Alignment Stability 

• Scale Factor 

• Sensor Translation 

• Sensor Orientation 
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■*5 ppm 
0. 00 1 inch 
0. 00 1 radian 
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Section 4 - Sensor Design 

Subsection 4. 2 - Sensor Design and Performance Analysis 
DESIGN OF THE SENSOR FOR HIGH STABILITY 

The sensor stability requirements can be exceeded by winding the triaxial- 
feedback coils on the spherical graphite epoxy forms and by keeping the 
coil powe r low enough to red uce the thermal effe c ts. . 

The design of the sensor depends on a combination of magnetic, ther- 
mal, and mechanical factors. After so much analysis it may be helpful to 
review the overall design approach by way of a summary. As indicated 
previously, our analysis and that of others shows that high axial stability 
is much more easily attained if the three fluxgate sensors are enclosed in 
a single triaxial (vector) feedback field, rather than being enclosed in 
three separate single-axis feedback solenoids. A detailed analysis of the 
triaxial feedback scheme and a comparison of triaxial and single axis feed- 
back has been presented. The principal error source in the single -axis 
system is the instability of alignment of the sensitive core; this sensitivity 
is reduced by a large factor by triaxial feedback. Other errors, such as 
the axial instability of the feedback coils, are comparable for either system. 

Design of Feedback Coil/Sensor System - If the field inside the tri- 
axial feedback coils were perfectly uniform, the position and orientation 
of the sensors would be almost irrelevant. We have selected the simplest 
feedback coil system which produces adequate uniformity for ring-core 
sensors large enough (3/8 inch, or about 1 cm, in diameter) to meet t e 
noise specification — namely, a Helmholtz coil pair along each axis. A 
larger pair is used to cancel the external field, as discussed previously 
and shown in the figure opposite. 

The ring-core fluxgate sensors are mounted orthogonally on a cube. 

They are made as small as possible consistent with meeting the noise 
specification. Their dimensions were adjusted to minimize the noise- 
equivalent field subject to the geometric size constraint. The inner feed- 
back coils are pushed together slightly to make the field flat at the sensors 
rather than at the geometric center of the coils. 

Mechanical and Thermal Design of the Sensor Assembly - The primary 
consideration in the design of the sensor assembly is mechanical and ther- 
mal stability of the feedback coils, which is the main factor in determining 
the gain, linearity, and alignment stability of the instrument. It has been 
determined that a graphite -fiber epoxy-resin material will provide the 
required mechanical and thermal stability. 

The sensor assembly is composed of a base, a sensor tower, and two 
spheres of graphite -epoxy on which the coils are wound. In order to 
achieve thermal stability, a thermal resistance and coil power must be 
selected that will result in an acceptable thermally induced variation in the 
feedback field. The proposed design produces a thermal variation of the 
field of less than 0. 17 with a coil power of 20 mW. 

Reiec tion of Thermal Gradients - The sensor enclosure should be iso- 
thermal^ since gradients of 1°C lead to field variations in the 0. IV region 
(because of imperfect cancellation of the copper coil thermal expansion). 

A thermal shield is required to reduce the internal temperature gra- 
dients. We use a shield which is a "tee attenuator" - it consists of two 
layers of multilayer low-emis sivity superinsulation separated by a graphite- 
epoxy housing which performs the dual functions of thermal shunt and cover. 

We will now describe the design in more detail, starting with the flux- 
gate assembly. 


4-28 



( 


time-zero 

TABLE 4-14. FEATURES OF THE RECOMMENDED SENSOR DESIGN 


Design Feature 

• Triaxial Feedback 

• Small Ring Cores (3/8 inch, ~1 cm) 

• Helmholtz Coils 

• Large Helmholtz Cancellation Coils 

• Spherical Coil Supports 

• Thermal shield 


Motivation 

- Greatly improves axial stability. 

- Meet noise specification. 

- Reduced alignment sensitivity. 

- Adequate uniformity; even 
better if adjusted to optimize 
flatness . 

- Exactly cancel dipole field. 

- Very low internal gradients. 

- No quadrupole field, 

- Lower power. 

- Single rigid structure. 

- High thermal conductivity. 

- Acceptable weight. 

- Isothermal enclosure. 



Figure 4-9. Conceptual Sensor Arrangement. The three orthogonal 
fluxgates are mounted at the center of two coil systems: one to pro- 

duce the feedback field, and one to cancel the external field. 
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Section 4 - Sensor Design 

Section 4. 2 - Sensor Design and Performance 
DETAILED DESIGN OF THE FLUXGATE SENSORS 

The 3/8-inch ring-core fluxgates use half-mil 6-8L3 Mo-Permalloy tape 
wound on an Inconel bobbin. The tape height is 1/8 inch. The sensors are 
mounted orthogonally on a cubic mount. 

A specific design for the fluxgate sensor to be used has been devel- 
oped as a baseline for the study. The sensor described here has been 
made as small as possible consistent with the budgeted noise. 

The fluxgate sensor is essentially a magnetic amplifier. It consists 
of a square loop core which is made to saturate alternately in each direc- 
tion by an alternating current in an ac drive winding. In the presence of 
an external field, the two halves of the core saturate at slightly different 
times, giving rise to a flux change at twice the drive frequency which is 
sensed and amplified to obtain a measure of the external field. 

The sensor exhibits both noise and offset; the theory of these effects 
is discussed in Appendix A. The considerations required to obtain low 
noise and offset are somewhat different than for a magnetic amplifier 
because the signal field is external; it is important to design the sensor 
to intercept as much of the external flux as possible using a minimum 
amount of noisy core material. 

Design Considerations - To achieve a high signal-to-noise ratio with 
low drive power one must use a core material with a high permeability 
and square hysteresis loop such as molybdenum Permalloy. The best 
results have been obtained with an alloy consisting of 6% Mo, 81. 3% Ni, 
and 13.7% Fe; this alloy has zero longitudinal magnetostriction. The 
method by which the material is prepared and annealed is important. 

Because the Permalloys are very sensitive to stress, extreme care 
must be taken in mounting the core to avoid offsets. The procedure used 
at NOL and TZL employs a single spot-weld to hold the tape on an Inconel 
bobbin whose thermal expansion matches that of the tape. It is not a good 
idea to clamp or pot the tape, so it is difficult to guarantee second-of-arc 
alignment stability by design — hence, the use of triaxial feedback. 

Specific Design for MAGSAT - One-inch ring-core sensors built in the 
above fashion exhibit typical noise levels of 0. IV peak for a 25-Hz band- 
width (about O.OlYrms). The theory in Appendix A shows that the noise- 
equivalent field scales as the 3/2 power of the diameter and the square 
root of the tape height. Thus, the diameter can be reduced by a factor up 
to 10^/3 if the noise can be as large as IV peak— or to about 0. 3 inch. 

The allowable size of the fluxgate is set by the size of the flat-field vol- 
ume in the triaxial feedback coils (see following discussion); a diameter of 
about 1 cm, or 3/8 inch, is acceptable.* We predict a noise level of about 
0.8Ypeak for this size sensor, within the 1 V specification. The offset 
variation will be close to the ±0. 4Y specification, so one would not want to 
make the sensors any smaller. 

The tape is mounted on a thin-walled Inconel bobbin. A thin bobbin 
wall thickness is desirable to reduce the drive- circuit leakage inductance 
and also to reduce the eddy currents in the Inconel. 


Because tape and other mechanical dimensions are still specified in indus- 
try in English units, we have used these units rather than SI units in what 
follows. 
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The drive winding is then wound on the toroidal bobbin and the sense 

uu er 11 With ° Ut a f ° rm to minimize the height. A plug is 
used to fill the bobbin to provide support using a flexible epoxy to avoid 
stressing the bobbin and tape. 1 

The electrical characteristics of the sensor are shown in Figure 4-10. 

cnh ^Th g - The sensors are mounted on the surface of a 

* , The hei g ht -to-diameter ratio of the sensor core tape has been selected 

ontiri mm f 1Ze the sl § nal - to - noise ratio with a specified size constraint; this 
optimization leads to the conclusion that the tape height should be 1/3 of the 
tape winding diameter. 

S ensor Manufacture - TZL has fabricated such sensors in-house using 
a released manufacturing procedure based on the NOL procedure, and has 
aiso recently purchased similar sensors. The procedure is not difficult pro- 
viding the material and the assembled sensor are annealed properly. 
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, 6-81.3 MoPERMALLOY TAPE 
ON INCONEL BOBBIN 

TAPE DIMENSIONS 0.375 INCH DIAMETER 
0.125 INCH HEIGHT 
0.5 MIL THICKNESS 

- RECTANGULAR SENSE WINDING 
200 TURNS #36 5 LAYERS 
L - 400 jiH, R - 10ft, S = 5 ^V/7 


■ TOROIDAL DRIVE WINDING 
100 TURNS #30 

L - 2MH-12mH, R - 0.4ft. I = 2 mA 


Figure 4-10. Fluxgate Sensor. The inductance and resistance are given 
for each winding as well as the sense winding sensitivity, S, and the drive 
winding saturation current, I s . y ' e arive 



Figure 4-11. Sensor Mounting on a Plastic Cube. The sensors have 
been designed to provide maximum signal to noise ratio with this 
mounting configuration. The sensors are spot-bonded to the cube. 

tr1axil° r f. a e5L 8 ack 8r “ ° f m °“ ntin 8 P” cisl °" « eliminated by use of 
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Section 4 - Sensor Design 

Subsection 4. 2 - Sensor Design and Performance Analysis 

MAGNETIC CONFIGURATION OF THE FEEDBACK COILS 

The simplest feedback coil system which provides the required uniformity 
and external field cancellation consists of two pairs of Helmholtz coils per 
axis with a slight modification of spacing to flatten the field at the sensor. 

As explained previously, the problem in designing the feedback coils 
is to provide a sufficiently uniform feedback field with a simple, mechan- 
ically stable configuration. The previously presented configuration analy- 
sis is reviewed briefly here. 

Requirements for Field Uniformity - According to the adopted error 
budget, translation of the sensor should give rise to repeatability of less 
than 0. 5Y, while rotation of the sensor should give rise to worst-case 
repeatability of less than 2. 5 seconds of arc or 0. 78Y; both values apply 
for a full-scale feedback field strength of 64000 Y. 

Methods for Generating the Field - As we have indicated previously, 
an almost perfectly uniform field can be generated by distributing many 
windings in the surface of a sphere uniformly spaced along the polar axis; 
the question is simply how many windings to use to achieve a given uni- 
formity. If the number of windings used is 2N, the field inside varies near 
the center as r^N when the windings are positioned correctly. One can do 
even better by adjusting the radius of each coil as well as the spacing. 

The simplest configuration is the so-called Helmholtz coil pair, consisting 
of two coils separated by their radius. For this spacing, the field at the 
center of the coils is flat to the fourth order in the distance from the axis 
(that is, the first three derivatives of the field strength vanish at the cen- 
ter of the field). A consequence of this constraint is that the quadrupole 
moment of the coils vanishes; there is no term in the far field varying as 
r - 5, a consideration important in meeting the external field specification. 
The external field has an r -3 dipole term, but the next (octupole) term 
varies as r"^. 

Cancellation of the External Field - By the arguments presented pre- 
viously, the outer coils must be designed to cancel the external dipole 
field of the inner coils, which means that the dipole moments must be 
equal and opposite. The dipole moment of a coil is equal to the product 
of the ampere-turns and area of the coil. Thus, if the outer coils are 
twice as large as the inner coils, they need only have one-fourth the num- 
ber of turns. The internal field of the outer coils is then one-eighth the 
field of the inner main coils. The external field is then essentially an 
octupole field, which is only 0. 016Y at a distance of 60 cm (24 inches) 
from the center of the coils when the central field is 64000Y. 

Adjusting the Field for Flatness at the Sensor - The sensors are not 
located at the center of the field but are offset from it. Motion in the direc- 
tion of this offset is important but motion in the other directions is neg- 
ligible by symmetry. 

One can adjust the coil spacing to reduce the gradients in the neighbor- 
hood of the sensor. Figure 4-12 shows the variation of the axial com- 
ponent of the field, in Y, as a function of radius in the offset radial direc- 
tion along the central plane. The calculation assumes that the outer 
cancellation coils have twice the radius of the inner coils and have a fixed 
spacing equal to their radius. The calculation was performed by a num- 
erical quadrature of the vector field integral. The sensor location is 
indicated. 
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The field gradient can clearly be eliminated by making the inner coil 
spacing slightly less than 0. 94 times the radius so that the field is sym- 
metrical across the sensor. An uncertainty of 0. 004R (0. 005 inch) in the 
feedback coil spacing leads to a sensor gradient of about 20Y/0. 04R, or 
4007 per inch. This gradient implies a sensor position repeatability of about 
0. 001 inch which is an easily met tolerance. 

Sensor Alignme nt Requirement - If the sensor normal to the field is 
cocked by a small angle A 0, the sensor output will be in error by the amount 
ABA6, where AB is the difference in the field seen by the sensor parallel to 
the field and the one perpendicular to it. Figure 4-13 shows the field as a 
function of the coil spacing. If the spacing is chosen to minimize the transla- 
tional error (the larger effect) the field difference is only 120Y, which leads 
to a sensor core alignment repeatability requirement of 0. 3 degree, a require 
ment which is 200 times easier to meet than 5 seconds of arc. 

This analysis shows that the performance of a Helmholtz coil system is 
more than adequate to provide a feedback field allowing good margins for 
alignment stability of the fluxgate sensors. 



Figure 4- 12, Variation in Axial Field for Various Radial Offsets in 
Central Plane. The sensor position is indicated. The field variations 
apply for an axial field of 64, 000Y. The field can be flattened by mak- 
ing the inner coil separation about 0. 935R. The field of the cancella- 
tion coil is included. 



Figure 4-13. Radial Field Difference Between Parallel and Perpendicular 
Sensors as a Function of Coil Separation. The difference of about 120 Y 
means that the sensor alignment need only be stable to about 0. 3°. 


4-33 




Section 4 - Sensor Design 

Subsection 4. 2 - Sensor Design and Performance Analysis 
ELECTRICAL DESIGN OF THE FEEDBACK COILS 


The main feedback coils will be about 2. 5 inches in diameter with a 2 mm 
cross section. The outer coils will be slightly larger than 5 inches in 
diameter to allow for mechanical clearances. The total coil power will 
vary from 5 to 20 mW, 


Having selected the basic coil geometry, we must now determine the 
coil power, winding cross section, and winding requirements. The design 
objectives are to minimize the coil power to reduce self-heating, and to 
minimize the coil cross section so that mechanical interferences will be 
reduced. The values calculated here are nominal, because in practice the 
coil diameters have to be adjusted slightly for mechanical clearance of 
coils on the three axes. 

Coil Power - The power allowable in the coils was estimated by a 
thermal calculation. The thermal expansion of the coils due to self-heating 
effect produces a non-linearity which must be less than 1. 5Y/64000Y, or 
23 ppm. One would like to see something a good deal smaller, to avoid 
thermal lags. We took 0. 5Y as a design goal, and the temperature coeffi- 
cient of copper (17 ppm/°C) as the relevant sensitivity parameter, imply- 
ing that the temperature rise should be less than 0. 5°C. A simple thermal 
model for the coil support was used to select a design goal of 5 mW for 
each inner coil, neglecting the complication of the outer coil. 

Determination of Coil Cross Section - For a given field-filled volume, 
the coil power depends only on the copper cross section. For a Helmholtz 
coil pair, the power is given in terms of the internal field by the following 
equation, in which p is winding resistivity, A w is winding cross section, 

R is coil radius, and is 4tt x 10"^ H/m: 



The formula is writen this way for a reason: the second expression involv- 
ing is proportional to the magnetic energy in the field, while the term 
involving p is dimensionally a frequency which depends only on the ratio of 
resistivity to the winding cross section. From this formula, one finds that 
for P = 5 mW per coil, A w should be 3.46 mm^, corresponding to 273 turns 
of No. 36 wire on each coil with a resistance of 71. 4f2, and an inductance 
of about 8 mH per coil. 

Drive Current - The drive current can be calculated from the formula 
for the field of a Helmholtz coil, which is 


B = 



3 

2 p NI 
*o 

R 


( 2 ) 


For B = 64, 000 Y and R = 1.25 inches (0. 03125 m), the value of NI is 
2. 26 A-turns, giving a value for the current of 8. 27 mA. Smaller wire 
can be used to raise the coil resistance and reduce the drive current, 
keeping NI constant. 

Correction for Cancellation Coil - The cancellation coil reduces the 
internal field by one -eighth, so that the drive current must be increased 
by a factor 8/7 to 9.45 mA. The cancellation coil power is half the inner 
coil power; the total power is 
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P = 


(2 x 5 + 2 x 



19. 6 mW 


(3) 


The coil constant is thus 6. 77Y/pA if #36 wire is used. The coil winding 
parameters are summarized in the schematic opposite. 

Power Variation in Orbit - The power input to the triaxial coil assembly 
never goes to zero in orbit as it does for a single-axis sensor. Since the 
power is proportional to the squared magnitude of the vector field, one may 
use the dipole equations to estimate that the power for an orbit crossing the 
magnetic pole varies as 


P ~ 


P 

— — ( 1 + 3 sin 
4 


2 


X) 


where X is the magnetic latitude, and P Q = 19. 6 mW. The maximum rate of 
change of power input occurs for a magnetic latitude of 45 degrees, when it 
is 5 pW/second. This result can be used to estimate the self-heating lag. 
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Figure 4-14. Coil Schematic. The inner and outer coil diameters 
are 2. 5 inches and 5 inches respectively. This design is nominal, 
since it does not allow for the slight differences in coil sizes 
needed to allow for mechanical clearance between coils. 
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Section 4 - Sensor Design 

Subsection 4. 2 - Sensor Design and Performance Analysis 

MECHANICAL AND THERMAL DESIGN OF THE SENSOR ASSEMBLY 

Designing with graphite -fiber epoxy-resin material produces a mechan- 
ically rigid and thermally stable structure. 

Mechanical and thermal stability of the feedback coils is a principal 
factor in determining the gain, linearity, and alignment stability of the 
instrument. The coil supports should also be light and easily assembled. 

The design must be rigid enough so that no permanent deformation is 
experienced after environmental testing. The thermal deformation result- 
ing from operating temperature change, self-heating, or external gradi- 
ents must be small enough to produce acceptable errors in the magnetic 
field. 

To meet these requirements, the structure must be designed with a 
minimum number of pieces, so that tolerance buildup at mechanical joints 
is reduced. 

The structure must be made with adequate cross section to provide the 
desired rigidity, and the pieces must be made from a light, thermally 
and structurally stable material with a close-tolerance construction. 

The design which meets these requirements is shown opposite; the 
table summarizes how each of the design requirements is met. 

Mechanical Design - The structure is composed of four pieces: base, 

sensor tower, and two spheres. Assembly of each sphere to the tower, 
and the tower to the base employs an interference fit. An adhesive is used 
as an additional retaining feature to provide a backup during dynamic 
environments . 

Each element of the structure is designed to optimize its individual 
requirements. The base plate provides a flat surface for mirror mounts 
on one side and a platform for the sensors, feedback coils, and field can- 
cellation coils. 

The sensors are mounted on a stepped cylindrical tower which is 
machine-finished to provide the necessary flatness for axis alignment. 

The feedback and cancellation coils are molded on spherical forms which 
provide a significantly higher degree of structural stability than could be 
attained by attempting to assemble individual coil forms into a triaxial 
arrangement. The spherical form also provides an inherent structural 
symmetry which provides resistance against the tendency to warp or twist. 

Selection of Material - An investigation of the nonmetallic materials 
capable of meeting design requirements leads to use of a graphite -fiber 
epoxy-resin system such as the GY-70 fiber-reinforced epoxy. When 
properly fabricated (proper fiber orientation), the material exhibits prop- 
erties such as: 

• Temperature coefficient less than 1 ppm/°C. 

• Tensile modulus >40 x 10& pounds per square inch. 

• Tensile strength >100, 000 pounds per square inch. 

• Thermal conductivity >most metals. 

Thermal Considerations - The feedback coils require high stability 
with regard to thermal effects. By using a thick-walled sphere as the 
coil form, three positive thermal benefits are obtained: first, the thermal 
conducting path is good, minimizing gradients; second, a large surface 
area is available for radiative heat transfer; and third, thermal stresses 
can be more uniformly distributed within the structure. 

To further minimize inputs to the feedback coils, an active /pas s ive 
thermal control system is provided. The system contains multilayer 
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insulation, a heating element, and a thermal shunt. Design of the thermal 
shield and controller is discussed in another topic. 

This design has been analyzed to determine its mechanical and thermal 
performance; an analysis of the structural stability is next presented. 

In addition, a detailed evaluation of the recommended structure and mate- 
rial selection was performed. The results of this evaluation appear in 
Appendix B. 


TABLE 4- 15. DESIGN REQUIREMENTS AND APPROACHES 


Requirement 

Approach 

• Lightness 

- Simple, rigid structure. 

- Graphite fiber material. 

• Ease of Assembly 

- Four-piece assembly. 

• Rigidity 

- Adequate cross sections. 

- Graphite fiber material. 

- Close tolerances. 

- Press-fit construction. 

• Thermal Stability 

- Thick-walled spherical forms 

- Graphite fiber material. 


- Thermal control system. 



Figure 4-15. Sensor Mechanical Design. A high degree of stability is 
attained through the use of a simple spherical form and graphite-fiber 
epoxy-resin material. 
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Section 4 - Sensor Design 

Subsection 4. 2 - Sensor Design and Performance Analysis 
ANALYSIS OF STRUCTURAL STABILITY 

Analysis of static one-g loads and dynamically induced deflections and 
stress shows that the sensor will exhibit excellent repeatability and very 
low static misalignments. 

The analysis was performed to evaluate the effect of all mechanical 
and thermal factors which could have an impact on sensor stability. Empha- 
sis was placed on factors affecting the alignment of the field inside the tri- 
axial feedback coils, as these are the principal factors affecting sensor sta- 
bility. The details of the analyses presented here appear in Appendix C. 

The analysis examines two areas: stability in zero- to one-g fields 

requiring stability of ±0.5 arc-second per degree Celsius, and dynamic 
loads of thermal cycling, shock, acceleration and vibration which require 
cumulative stability of ±2.5 arc-seconds (repeatability). The stability fig- 
ures are those assigned by the error budget discussed previously. The 
analysis also considers effects of the use of adhesives and creep, a trouble- 
some property of viscoelastic materials. 

Static Analysis - Static feedback coil stability is affected by mechanical 
deflection of the spherical coil form. The coil form nominal diameter is 
2.5 inches with 0.2-inch thick walls and weighs 0.43 pounds. Assuming a 
tensile modulus of 10 x 10^ pounds per square inch, the maximum static 
deflection which occurs in test at 1 g (in space it is 0 g) is less than 
4. 0 x 10-7 inches/inch (corresponding to 0.01 arc-second). This value is 
better than the requirement value by more than an order of magnitude. 

Another effect on static feedback coil stability is deflection of the coil 
with respect to the mirror mount. This value is 1.0 x 10" 7 inches/inch, 
more than an order of magnitude better than the required value. 

Dynamic Analysis - The design approach used for the structure pro- 
duces low stresses during maximum dynamic inputs. The low stresses, 
coupled with the high tensile modulus (42 x 10*° psi parallel to the fibers) 
permit repeatability within ±0.2 arc- second (the mechanical budget is 
±2.5 arc-seconds). The lowest resonant frequency of the structure exceeds 
4000 Hertz, precluding amplification of dynamic inputs during vibration 
(maximum vibration test is 2000 Hz). 

The maximum dynamic deflection (assuming a 28 g input) within the 
structure occurs at the location of the 3 axis sensors, and is 9.52 x 
10-6 inches per inch. This maximum deflection is far less than even the 
allowable static deflection of 0.001 inch. All other deflections contribute 
less than 0.2 arc-second error. The maximum dynamic stress of the 
assembly which occurs in the graphite-epoxy baseplate, is 189 pounds per 
square inch, which is insignificant when compared with the baseplate 
material yield strength of about 50,000 pounds per square inch. 

Dynamic Thermal Stress - The temperature coefficient of the graphite- 
epoxy coil .form is negligible compared to that of the copper winding. 
Because the coils will be wound at room temperature (~21°C), the coil 
form diameter will try to contract or expand about 0.005 inch when exposed 
to the specification temperature extremes of -50°C and +85°C. Analysis 
shows that at -50°C, stress in the coil wires will exceed the wire yield 
stress (-10,000 psi) and actually exceed the maximum allowable tensile 
stress (35, 000 psi) if the coil diameter is prevented from shrinking. This 
condition can be avoided by providing a thin flexible elastomeric layer as a 
base for the coil windings, allowing diametral changes without exceeding 
the yield stress. Elastomeric adhesives will be used to hold the individual 
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coil wires in place at temperatures above 21°C so that the wires will return 
to the same relative positions at the normal operating temperature of the 
6ensor. By design, the ±2.5 arc-second repeatability is met. 

Adhesives . - Adhesives are used only as a backup to the assembly, which 
utilizes interference fits. The adhesives do not contribute to stresses or 
deflections. 

Creep Effects - Creep can be expressed as continuous deformation of a 
viscoelastic material under load beyond the initial deformation, measured 
immediately after the load is applied. Because creep rate is proportional to 
elastic modulus, a material with a high modulus exhibits a minimum tendency 
toward this effect. By comparing the graphite-fiber epoxy-resin system with 
other viscoelastic materials, a conservative estimate of the maximum creep 
rate can be determined. Using a value of creep rate equal to 0.002, we can 
equate this to a decay in the elastic modulus of 2% in a 5-year period. It fol- 
lows that the effect on stability will be 2% additional error. Because stability 
is well within acceptable limits, a 2% additional error is negligible. 

This analysis shows that the structure is far stronger than required to 
withstand the environmental stress and to exhibit negligible differences in 
alignment between 0 and 1 £. 


TABLE 4-16. ANALYSIS OF CRITICAL ELEMENTS 


Element 

Analysis Type 

Maximum 
Deflection 
or Stress 

Sgec. 

Result 

F eedback Coils 

Static 1 £ 

0.02 arc- second 

±0. 5 

Meets spec. 

Feedback Coils 

Dynamic (28 £) 

0.56 arc-second 

N/A* 

Meets spec. 

Sensor Cores 

Dynamic (28 jg) 

1,63 arc-seconds 

N/A* 

Meets spec. 

Baseplate 

Dynamic (28 £) 

189 pounds/ 
square inch 


Does not 
yield 

F eedback Coil 
Wires 

Thermal 

«2. 5 arc-seconds 

— 

Meets spec. 

F eedback Coil 
Mount 

Static 1 £ 

0.34 arc-second 

±0.5 

Meets spec. 


*The instrument does not have to operate in vibration, but note that 
the dynamic misalignments during vibration test are less than the 
specified misalignment tolerance of 2.5 seconds of arc. 
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Section 4 - Sensor Design 

Subsection 4.2 - Sensor Design and Performance Analysis 
STRATEGY FOR CONTROLLING THERMAL EFFECTS 

Mean temperature variations are reduced by use of active temperature con- 
trol, while the effects of gradients are eliminated by the use of isotropic 
structural material, thermal insulation, and a thermally conductive struc- 
ture with low internal power dissipation. 

To provide the absolute accuracy and alignment stability required for 
MAGSAT, it is necessary to avoid thermal expansion of the feedback colls, 
particularly distortions arising from unequal expansion of different parts 
of the structure. Even if the structure is monolithic, so that all parts of 
it have the same expansion coefficient, such differential expansions can 
result from temperature gradients in the structure. 

Need for Active Temperature Control - Even when thermally stable 
structural materials are used, there is still a requirement for temperature 
control. The overall temperature coefficient of the coil structure will have 
a non-zero thermal expansion coefficient because of an imperfect match 
between the thermal expansion of the copper and that of the coil forms; we 
estimate that the coefficient could be as large as 2 ppm/°C, leading to a 
field temperature coefficient somewhere between 1 and 2 ppm/°C by our 
previous sensitivity analysis. Because an overall stability of about 20 ppm 
is needed, with errors arising from several additional sources, mean tem- 
perature control to about ±1°C is desirable. Several methods of achieving 
this control were considered; as discussed later, it appears that an elec- 
trical heating scheme fits the requirements best. 

Material Selection - The graphite -epoxy material recommended has two 
virtues: very high elastic modulus and hence high rigidity, and high ther- 

mal conductivity — much higher than that of most plastics and at best com- 
parable to such metals as titanium. However, this material is often fur- 
nished as a laminate having high elastic and thermal conductivity only in 
the plane of the laminations. Use of such an anisotropic material for the 
coil forms and supports might lead to trouble as a result of unequal ther- 
mal expansion in different directions. For this reason, the recommended 
material is isotropic with randomly- oriented graphite fibers. 

Shielding for Elimination of Temperature Gradients - Temperature 
gradients arise from the heating of the magnetometer by the sun on one 
side and by cooling to space on the other side, as well as from self-heating 
of the coils and sensors. 

The elimination of gradients implies that the thermal model of the 
important parts of the instrument — namely the feedback coil structure and 
baseplate - should be represented to a high degree of accuracy as a single 
thermal node. This requirement implies careful insulation of this node 
from external influences, the use of a highly conductive structure, and the 
avoidance of self-heating by keeping the coil power input low. All of these 
requirements are met by the design. 

Solar Heating Effect - The orientation of the instrument and satellite in 
orbit is indicated in Figure 4-16. In the nominal near polar orbit, one side 
of the instrument will be continuously illuminated by the sun, while the 
other side will be continuously exposed to space (the eclipse case is dis- 
cussed later). By using a surface coating with a low solar absorptivity and 
high emissivity on the sun side (e.g., white paint), and (if possible) a sur- 
face with a low emissivity on the space side, the temperature gradient 
across the instrument can be minimized. Then, a superinsulation and a 
conductive housing and coil structure can provide the rest of the isolation 
required. 
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The main effects then become heat losses fromthe first- surface alignment 
mirrors mounted on the baseplate and losses to the uncontrolled instrument 
mounting shelf. Gradients produced by this cooling effect can be reduced by 
applying heat behind the mirrors using an appropriate printed-circuit heater, 
and by isolating the baseplate from the instrument mount. This thermal 
design is discussed in more detail later. 

Self-Heating - The coil and sensor power was selected to provide accept- 
able gradients from self-heating effects. The analysis presented later shows 
that in the recommended design such effects produce negligible changes in the 
measured field. 

We now turn to the details of the thermal design and performance analysis. 

TABLE 4-17. THERMAL CONTROL STRATEGY 


Effect 

Control 

• Mean Temperature 

- Active thermal control with 

Variation 

electric heater. 

• Solar Heating Gradients 

- White surface coating. 

- Super insulation. 

- Conductive housing, structure 
and baseplate. 

• Backplate Gradients 

- Heater location. 

- Housing insulation. 

• Self-Heating Gradients 

- Low coil power. 

- Conductive coil structure 
and baseplate. 


NAS 5 23661/4 16 



Figure 4-16. Orbital Heating. In a dawn-dusk orbit, one side 
of the instrument will be continuously illuminated. 


\ 
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Section 4 - Sensor Design 

Subsection 4.2 - Sensor Design and Performance Analysis 
ANALYSIS OF SELF-HEATING EFFECTS 

Analysis shows that self-heating errors are negligible, being less than 

5 x 10“ 3 Y for the largest effect, due to the use of high thermal conductivity 

graphite composite structural materials. 

We now try to estimate the thermally induced errors for the most cri- 
tical sensor component, the feedback coils and the supporting structures. 
Since 7/8 of the feedback field is produced by the inner set of Helmholtz 
coils, their behavior is seven times more critical than that of the outer 
cancellation coils, and primary consideration is given to the inner coils. 
Operating temperature, gradients, and self-heating can affect the gain, 
linearity, and alignment through geometric changes due to thermal 
expansion. 

Temperature Rise Due to Self-Heating - An approximate thermal model 
shown opposite has been analyzed by linearizing the radiation conduction 
equations between the three major components of the sensor: flux gate 
cores, feedback coils, and field cancellation coils. The cores, which gen- 
erate 15 mW, maximum, are conductively coupled through their support to 
the baseplate reference with a resistance of 40°C/W and radiatively coupled 
to the surrounding feedback coil sphere with 224°C/W. The feedback coils, 
generating 13 mW maximum, are conductively coupled to the reference via 
4°C/W and radiatively to the field cancellation coils by 13°C/W. The field 
cancellation coils with 6. 5 mW also in turn conduct through 4°C/W, and we 
neglect for the moment any radiation terms to the outside. The solutions 
for this network are T^ = 0. 5°C, T^ = 0. 06°C, and T^ = 0. 04°C due to 
maximum self-heating. 

Thermal Expansions - We use an effective coefficient of thermal expan- 
si orTofT^prnT^C^for^he support structure calculated from the coefficients 
for copper and graphite-epoxy and their relative tensile strengths. Trans- 
lation of the sensor cores, with a sensitivity of 0. 5 Y /mil, changes in feed- 
back coil constant (gain and linearity) with AB/B ~ 1 ppm/°C, and alignment 
sensitivity of 1 ppm/°C yield errors as summarized in the table for the 
maximum internal field of 6. 4 x 10^Y # The effects of self- heating are seen 
to be well within the budgeted errors and are in fact negligible. 

Thermal Lag and Other Effects - The heat capacity of the structural 
material used is high and results in an estimated thermal time constant of 
300 seconds. Since the total self-heating errors are so small, errors due 
to thermal lag in a changing field are also negligible. Also, with an elec- 
trical resistivity of ~3000 p ohm-cm, yielding a surface resistance of 

6 x 10“ ^ ohms /square, both eddy current heating in the spherical shell of 
the coil supports and field screening from the interior are negligible for 
rates of 250 Y /sec. 

With virtually no self-heating errors as determined by the above anal- 
ysis, we next consider design of the thermal control system. 
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TABLE 4- 18. THERMAL ERROR SUMMARY 


Source 

Self-Heat 

Error 

Temperature 

Sensitivity 

Budgeted 

Error 

Sensor Core 
T ranslation 

2. 5 x 10 " 4 Y 

5 x 10" 4 Y/°C Tj 

±0. 5 Y 

Gain and 
Linearity- 

5. 1 x 10" 3 Y 

6.4 x 10" 2 Y/°C T 2 

±0. 5Y 

Feedback 
Axis Alignment 

5 . 1 x 1 0 ~ 3 Y 

0. 2"/°C gradient 

±0. 5 arc- 
seconds 



Figure 4-17. Approximate Thermal Model. T,, To, To are the tern 
peratur.s of the fluagate cores, feedback coils,* and^eld cInc.U,u”n 

coive'rtTf r eCl ‘ V ! I ly ' *" * he matrix «q»«ion resistances have been 
converted to conductances with units of milliwatts and degrees Celsius. 
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Section 4 - Sensor Design 

Subsection 4. 2 - Sensor Design and Performance Analysis 
DESIGN OF THE THERMAL CONTROL SYSTEM 

A combination active/passive thermal control system prevents excessive 
heat interchange with the sensor feedback coils. 

The thermal interchange between the sensor and its external environment 
must be minimized to keep misalignment and gain errors within limits. The 
major problem is heating by the sun on one side and cooling to space on the 
other side. To stay within the error allocations specified, the design must 
include the following features: 

Super insulation - This recommended insulation is actually multilayer 
aluminized or goldized mylar or kapton. The purpose of the insulation is to 
provide a very high thermal resistance path from the external environment 
inputs (solar and space shadow) to the sensor field cancellation coils and 
feedback coils. The insulation is about 0. 5 inch thick and will be placed on 
the outside of the sensor and in the space between the field cancellation coil 
and the thermal shunt. The insulation is specially designed for spheri- 
cal shape to avoid the wrinkling effect obtained if a flat shape is wrapped 
around a sphere. The outer surface, on the sun side, is painted white 
to reduce the temperature difference across the instrument. If the orienta- 
tion can be guaranteed, the finish on the space side should be gold or alu- 
minum. If no such guarantee is possible, painting th entire outside sur- 
face white is acceptable. 

Thermal Shunt - The thermal shunt is a graphite-fiber epoxy-resin 
sphere. The principal purpose of the shunt is to equalize thermal inputs to 
the coil portions of the sensor. This is possible because the shunt provides 
a good heat conduction path so that the temperature gradient across the 
shunt is minimized. A thermal control finish is applied to the exterior of 
the shunt. 

Base and Mirror Shield - The mirror and part of the sensor base 
should not be covered by the superinsulation for obvious reasons. The mirrors 
must be exposed for alignment by a spacecraft light source and the base 
must interface with the boom. The rest of the base plate is covered by 
superinsulation. To minmize sun and shadow radiative heat transfer, a 
shield is provided to limit the field of view to primarily the spacecraft. 

Thermal insulation between the sensor and boom interface will minimize 
conduction heat transfer. It is recommended that the portion of the boom 
interface with the sensor also be insulated and that the boom thermal con- 
ductive area be small. 

Base Heating - Because the passive thermal protection provided by 
the superinsulation and thermal shunt is only partially available at the base, 
active thermal control is used (see next topic for design of thermal con- 
troller). The mirrors will lose heat by radiation and hence an active heat- 
ing element is necessary. A printed circuit type film heater and tem- 
perature sensors attached to the sensor base can provide the required base 
temperature. (Other heaters attached to the thermal shunt can be used to 
reduce the gradients even further. ) 

Performance - The multiple-layer approach attenuates the solar flux 
so that the direct effect on the feedback coils is negligible. The main 
source of temperature gradients then results from heat flux on the highly 
conductive base plate, as we will now see. 

) 
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The fir st- surface mirror losses are a potentially serious source of loss 
even though its emissivity is low — 100 mW for a 10 cm x 10 cm mirror with 
an emissivity of 0. 02. By placing the heater behind the mirror, the unbal- 
anced heat flux due to this concentrated source of heat flux are reduced to 
levels of about 1 mW/cm^, leading to low baseplate gradients. 

TABLE 4-19- THERMAL PROBLEMS AND CONTROL APPROACH 
Problem 

• Radiation to Sensor Coils 

• Sun/Shadow Gradients 


• Exposed Base and Mirror 
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Figure 4-18. Thermal Control System. A thermal shield protects 
the sensor from external thermal effects. 
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Control Approach 

- Passive, Wrap with Superinsulation 

- Passive, Use a thermal shunt with 
thermal control finishes. 

- Passive, Shield Base and mirror to 
minimize thermal inputs. 


Active, Provide heating element to 
keep base at required operating 
temperature. 






Section 4 - Sensor Design 

Subsection 4. 2 - Sensor Design and Performance Analysis 
DESIGN OF THE THERMAL, CONTROLLER 

Use of a printed circuit heater is recommended to warm the instrument 
baseplate. The heater should be designed for low stray fields and driven 
by an ac current, using a thermistor sensor and control amplifier to pro- 
vide mean temperature control to ±1°C. 

Temperature control of the sensor is needed to prevent variation of the 
scale factor as well as instability of the sensor axis alignment. The con- 
trol requirement is about ±1°C, although one should be able to do better 
than that without much trouble. 

Selection of Temperature Control Scheme - Three different control 
schemes were considered — gas flow, passive control by a moving louver, 
and electrical heating. The gas flow method, using cooling and heating of 
gas piped up the boom, was rejected because it was difficult to insulate the 
gas-carrying tubes. The use of a thermostatic mechanism to control a 
moving louver is an attractive alternative, but the best place to locate such 
a mechanism is on the baseplate, which is reserved for alignment mirrors. 

As discussed previously, a principal mode of heat loss is radiation from 
the alignment mirror surfaces. It makes sense to distribute heat behind the 
mirrors, in effect guarding them and eliminating the gradients which the 
concentrated load would otherwise produce in the baseplate. This uniform 
heat distribution can most effectively be provided by use of a printed cir- 
cuit heater. 

Avoiding Magnetic Fields - Use of a dc heater drive would obviously be 
asking for trouble, because of the stray fields that might be produced. An 
ac drive should be used, at a frequency which is midway (geometric mean) 
between the 25-Hz low frequency response and the 10-kHz drive frequency — 
or about 1 kHz. 

Use of such an ac drive does not guarantee lack of stray field effects, 
because the sensor is after all a highly non-linear device and cross- 
modulation products may result if the stray field is too large, whether it is 
ac or not. To reduce the stray field of the heater, it should be laid out using 
a bifilar pattern with high spatial frequencies in two dimensions. 

The formula opposite gives the spatial dependence of the vector poten- 
tial produced by a current distribution periodic in the x-y plane. To satisfy 
the vector Laplace equation, the fall- off of the z direction must be exponen- 
tial. One readily estimates from this formula that if the printed circuit 
conductor spacing is periodic with a typical easily-attained spacing 
s/2 = 0.050 inch, the field at the sensor location (2. 5 inches away) will be 
less than the field near the heater by a very large factor (10^). The field 
produced by such a periodic conductor layout will thus be very small. The 
primary source of stray field will be the dipole field produced by unavoid- 
able current loops occurring where connections are made; estimates of 
these stray fields lead to field strengths of 0. 2 V at the sensors. Such a 
field strength is acceptable, especially for ac fields, so we conclude that 
a bifilar pattern heater will indeed work. 

Controller Design - The design of the control circuit is straightforward 
and so has not been worked out in detail. Use of a GSFC PPL thermistor 
is proposed for the sensor. 

It is recommended that the controller be a variable duty-cycle type to 
achieve efficient operation at constant frequency. The heater can be 
switched directly across the spacecraft supply with a switch conduction 
interval determined by the controller error signal. Chokes in series with 
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the line and load can be used to eliminate harmonic currents which might 
enter the signal band* A maximum power of 200 mW can be delivered by 
switching a current of 12. 5 mA with a 16 V power supply. 

This design appears to offer no difficulties in implementation. 

TABLE 4-20. SPATIAL DEPENDENCE OF VECTOR 
POTENTIAL FROM PERIODIC CURRENT 
DISTRIBUTION 


If the current distribution is periodic in the x-y plane with 
periods, the two components of the vector potential fall off 
exponentially in the z direction, perpendicular to the heater 
plane. 


-k z 

A ~ A ~ £e sin k x sin k y 
x y xx 7 
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Figure 4-19. Temperature Controller Block Diagram. 
The heater is driven by an ac current to avoid dc fields. 
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Section 4 - Sensor Design 

Subsection 4. 2 - Design and Performance Analysis 
PERFORMANCE OF THE THERMAL. SHIELDING 


* 


The successive layers of radiative and conductive shields eliminate tempera- 
ture gradients across the feedback coils. Temperature variations across 
the baseplate must be controlled by appropriate thermal design of the heaters 
and the boom-baseplate interface. 

We can estimate the performance of the thermal shield by a simple 
successive -approximations calculation. 

Assume that the thermal shunt is controlled to 300 K, and that the 
effective emissivity of the surrounding thermal blanket is 0, 005 -- equivalent 
to an ideal 4-layer blanket with surface emissivity 0. 02. The outer surface 
of the blanket is painted white, with a solar absorptivity 0. 10 and a thermal 
emissivity of 0. 60. 

The temperature Tj at a point on the illuminated surface at an angle 
0 from the subsolar point can then be found from the equation 

0£ s T l 4 + < *b< T l 4 - T 2 4) = a S F S ‘ :099 

which expresses the radiation balance on the outer blanket surfaces, neglect- 
ing conduction parallel to the blanket. In this equation, O is the Stefan- 
Boltzmann constant, e g is the surface emissivity, is the effective blanket 
emissivity, T£ is the inner (shunt) temperature, a is the solar absorptivity, 
and F fl is die solar flux (137 mW/cm^). Numerical evaluation yields a value 
of Tj equal to 253°K at the subsolar point and 90. 5°K on the shadowed side. 

The flux into the shunt (assumed spherical for purposes of the estimate) 
can then be estimated by assuming that the heat transfer is essentially radial 
and that the shunt is nearly isothermal (it is designed, of course, to make that 
assumption valid). The net heat transfer on the shadowed side is, for an 
average blanket radius a = 7. 5 cm (3 inches). 

jFdA = 2m 2 a€ b (T 2 4 - Tj 4 ) = 80. 5 mW 

On the sun side, one must perform an integral: 

JFdA = Joe b LT 2 4 - T 1 4 (0) ]2TTa 2 sin0d6 

Using Simpson 1 s rule with interval 0 = 45°, one obtains the result 

JFdA = 57. 3 mW 

The flux on both sides is outward, so that the net flux across the shunt is 
23.2 mW. 

Feedback Coil Gradient - We can then estimate the temperature gradient 
across the shunt and the inner and outer coils using the thermal model shown 
opposite. The values of the thermal resistances in the model were estimated 
from the dimensions and thermal conductivity of the graphite -epoxy shunt 
and coil forms. If the shunt has series and shunt resistances of 100°C/W, 
as shown (obtained with a thickness of about 1 mm and a thermal conductivity 
of 0. 02 W/cm^ - °C), the net gradient across the shunt will only be about 
7°C. Then, the high thermal resistance of the inner radiation shield and 
low thermal resistance of the outer coil form combined to attenuate this 
gradient to 0. 1°C, assuming only a single-layer inner radiation shield. 

This difference is so small that the gradients must be dominated by gradients 
in the baseplate. 
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, f ^f^ Se ? late ? r < adl tu tS ’ Slnce the thermal control strategy is to shunt 
e to the baseplate, there will inevitably be gradients in the baseplate. 
However, since the component must be thick and heavy for structural reasons 
ermal flux will not lead to large temperature gradients. The estimated 
temperature o gradient is on the order of 2°C across the entire baseplate or 
perhaps 0. 5 C across the feedback coil mounting flange. Assuming that this 

tui a COil “Alignment error with f £mpe~- 
second of arc. PFm/ ^ C » the corresponding misalignment would be 0. 2 

To summarize, there appears to be no difficulty in shielding the critical 
components of the vector magnetometer so that the effect of gradients pro- 
duced by sunshine are negligible. Careful attention must be paid to control 
of baseplate gradients and design of the boom interface. The^nterface must 
be designed to produce low gradients in the baseplate. 

Because of the high degree of thermal isolation provided the Derfor 

E TnSct rb blc W a ith " P artial f de « ree of shadowing ( P .clip“ wfub^ similar - 
Detter, m fact, because of greater symmetry. 

Although the calculations presented are based on a simple model and 

rLiTiremlnT 1 !?^ 16 ! 8 ° f therma . 1 resistance, the estimated performance exceeds 
equirements by a large enough margin so that no difficulties are to be exnect^H 
when the calculation is refined to improve its accuracy. expected 


TABLE 4-21 ESTIMATED TEMPERATURE GRADIENTS ACROSS COMPONENTS 


These values assume the use of 
the thermal shield. 


a distributed heater on the mirror but not 


on 



o Shunt 7°C 

o Outer coil form 0. 1°C 

o Inner coil form 0. 02 °C 

o Baseplate 1-2°C 


OUTER BLANKET 



Figure 4-20. Thermal Model 
for Temperature Gradient 
Analysis. Successive radia- 
tive and conductive layers 
suppress temperature gradients. 
To reduce the gradients to the 
maximum extent, the mirror 
losses and surface losses should 
be supplied by distributed heaters. 
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Section 4 - Sensor Design 

Subsection 4. 2 - Design and Performance Analysis 
ESTIMATES OF SENSOR STABILITY 

The estimated scale factor stability of the sensor is ±0.067, with an align- 
ment stability of ±0. 6 arc-seconds. The sensor offset is expected to exceed 
the ±0. 27 budget slightly. 

For convenience in comparing the estimated and specified performance, 
we have collected the results of our analysis here. 

Our performance estimates are based on simple mechanical and ther- 
mal models using typical values for the thermal and mechanical parameters 
of the material used. The estimates will need to be refined in a hardware 
development program using more sophisticated calculational tools. However, 
the estimates presented show that the predicted performance is far better 
than required. We then expect the triaxial feedback sensor to exceed the 
MAGSAT requirements. 

The various specifications are compared with the estimated perfor- 
mance in the table to the right (Table 4-22). 

Offset - The fraction of the budgeted offset drift ascribable to the 
sensor is ±0. 27. As noted earlier, smaller sensors have larger offset 
drifts. We did not build and test a sensor of the appropriate size, but do 
not believe that meeting the offset requirement will be difficult in spite of 
the fact that the sensor is smaller. 

This judgement is based on the NOL test results, in which offset 
variations of±0. 17at constant temperature for a one inch sensor were 
reported. We predict ±0. 37 for our sensor. This level of offset drift 
should not be too significant; it is less than the noise, and an insignificant 
contributor to the overall inaccuracy. 

Scale Factor - The scale factor error caused by thermal and stru- 
ctural effects is quite small, even though the estimates are very much worst- 
case. 

The principal source of sensor non-linearity, assuming adequate loop 
gain in the electronics (treated elsewhere), is self -heating. The use of thick 
walled graphite -epoxy coil forms makes this effect almost negligible. 

Alignment Instability - Alignment instabilities arise from structural 
and thermal distortions which are not easily calculated. The estimates 
shown are again very much worst-case assuming that a given thermal or 
structural expansion translates directly into a misalignment with unit 
efficiency. The most important effects appear to be baseplate deflections 
from thermal distortion and, as a result of the finite weight, of the align- 
ment mirrors, but the reduced finite effects are well within the allotted 
error budget. 

Summary - The design appears to yield performance as follows: 


• 

Offset 

- ±0. 37 

• 

Scale Factor and 



Linearity 

- ±0.067 

• 

Alignment Stability 

- ±0.6 arc-seconds 


\ 

I 
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TABLE 4-22 COMPARISON OF ESTIMATED 

TiME-zeaa 

AND REQUIRED PERFORMANCE 

Parameter 

Effect 

Estimated 

Performance 

Required 

Performance 

Offset 

Fluxgate Offset 

±0.3 7 

±0. 37 

Scale Factor and 
Linearity 

Coil Mean Temperature 
Solar Heating 
Self -Heating 

±0. 067 
±0. 0017 
±0. 0057 



RSS 

±0. 067 

±1. 57 

Alignment 

Coil Mean Temperature 
Solar Heating 
Self-Heating 
Baseplate Thermal 
Distortion 
1 -g Deflection 

±0. 2 arc -sec 
±0. 004 
±0. 01 

±0.4 

±0.4 



RSS 

— 0.6 arc - sec 

±2. 5 arc-sec. 
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Section 5 - Electronics Design 


SUMMARY OF ELECTRONICS ANALYSIS 

When implemented with carefully selected available components, a two- 
loop system employing a 12-bit analog-to -digital converter in the inner 
loop and an offset field generator in th$ outer loop enables meeting the 
MASAT Vector Magnetometer performance specifications. 

Two topologies were analyzed — a single-loop system employing a 
17-bit analog-to-digital converter (ADC) and a two-loop system employing 
an inner loop with a 12-bit ADC and an offset field generator (CFG) in the 
outer loop. The inner loop of the two-loop system is similar in many 
respects to the single loop system. 

Because of their similarities, initial performance analyses were per- 
formed on both systems. However, after the frequency and non-steady 
state field error analyses had been completed, and the two systems had 
been compared with respect to implementation difficulties and systemic 
characteristics, it was determined that the two-loop system was the best 
candidate for detail design and performance analysis. 

Study Tasks - The electronic study tasks fell into four general areas: 

• Analyzing the system characteristics of the single-loop and the two- 
loop system, as well as their frequency response and stability. 

• Selection of one system for detail design. 

• Selecting suitable components and topologies for the critical subsys- 
tems, which include the ADC, a digital-to-analog converter (DAC), 
and the feedback current source. 

• Synthesizing a complete system to enable a complete weight and 
power analysis. 

A summary of the study results in each of the above areas follows: 

System Tradeoff Considerations - Simplified block diagrams of these 
two systems are shown in Figures 5-1 and 5-2. The single-loop system 
is simple in concept, but difficult to mechanize as a result of the need for 
a 17-bit ADC (16 bits plus sign). The two-loop system is somewhat more 
complex functionally than the single loop because it requires an offset 
field generator (OFG) with its control circuitry. However, it has the off- 
setting advantages of putting the critical components in the low bandwidth 
OFG, which minimizes the effects of systemic noise and enables the use 
of less telemetry bandwidth. It also has better overload recovery charac- 
teristics due to a much shorter integrator time constant. 

Frequency Analysis - The frequency characteristics of the two-loop 
system are stable and easily mechanized. 

Critical Components - Both the two-loop system and the single-loop 
system require ultra- stable components (primarily resistors and 
operational-amplifiers) which additionally require active thermal control. 

System Implementation - An estimate of the weight and power of a 
totally implemented system shows that the specific weight, volume, and 
power requirements can all be met. 

The general conclusion of the results of each of the above areas is that 
the critical MAGSAT specifications are achievable. The analyses support- 
ing this conclusion are presented in the following pages. 
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Figure 5-1. Single-Loop System. The single-loop system requires a 
difficult to mechanize 17-bit ADC which is susceptible to systemic 
noise. 



DIGITAL 

OUTPUT 


Bi = INPUT MAGNETIC FIELD 
Ao = FORWARD LOOP GAIN 
0 = FEEDBACK GAIN 


ADC = ANALOG TO DIGITAL 
CONVERTER 

OFG = OFFSET FIELD GENERATOR 


Figure 5-2. Two-Loop System. The two-loop system has 
advantages overriding its additional complexity. 


5-3 







Section 5 - Electronics Design 
Subsection 5. 1 - Configuration Analysis 


CHARACTERISTICS OF THE SINGLE- LOOP SYSTEM 


Analysis shows that the single- loop system can achieve the required per- 
formance, but the wide bandwidth of the critical ADC and current source 
cause it to be susceptible to systemic errors (pickup, cross coupling, EMI 
impedances, etc.). These errors are difficult to analyze but are present 

in all systems. 


The block diagram of the single-loop system is shown in Figure 5-3. 
The transfer function of each block is defined in the following. 

Fluxgate Sensor - The nominal gain and expected uncertainty of the 
fluxgate sensor (which has been previously described) is 


K 

s 


3. 9 x 10‘ 6 


±0.2 x 10 


volts 

gamma 


a) 


Tuned Bandpass Amplifier - The minimum gain of the bandpass ampli- 
fier is determined by the allocated offset error from the total equivalent 
input drift of the demodulator and the integrator. Assuming that this drift 
is less than 2 x 10" 3 V, and using an error allocation of 0.07Y, then 


I » 10' 3 V 

K a > 0.07YK s 


3.66 x 10 


3 V 
V 


( 2 ) 


The actual design value for was chosen to be 

4. 07 x 10 3 V (31 

K - 71 v 

a (1 + 0. 5 X 10 s) 

Since the amplifier is tuned to the second harmonic of the drive signal, 
the equivalent pole must be included in its response. The anticipated 
worst-case variation in K a and in its associated time constant are both 

five percent. ... , . , 

Demodulator - The output of the ac amplifier is converted to a dc sig- 
nal by use of a synchronous demodulator. The dc transfer function of this 

circuit is 

K , = 0.63 cos <j>, 
d 


where <|> is nominally 90 degrees. . 

Some variation in phase will exist primarily due to variations in the 
tuned circuits of the ac amplifier. There is also a pole associated with 
the demodulator, whose value is selected to trim the closed loop fre- 
quency response, giving an overall nominal transfer function of 


k 


d 


0, 63 

(1 + 1. 84 x 10~ 3 s) 


V dc 
V peak 


The gain stability of the demodulator (including phase variations) will 
be better than two percent and its time constant will be stable to less than 

three percent. . , . . 

Integrator - The output stage of the forward loop is an integrator which 

provides the required high loop gain and also conveniently stabilizes the 
loop frequency response. The approximate response for the integrator is 
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1 V 

I T.s V (5) 

1 


where is selected to be 0. 582 second with an uncertainty of one percent. 
This time constant was selected to provide the desired closed loop response 
as described in the following section. 

Feedback, ..Networ k - The feedback network contains the feedback coil 
previously described and a current source for driving each coil. The band- 
width of the current source was made as wide as practicable. One pole, how- 
ever, does exist which is due to the coil time constant. The transfer function 
is given by 


K 

g 


1. 17 x 10~ 3 
(1+1.1x10 4 s) 


a_ 

V 


( 6 ) 


The dc errors of the current source are discussed in a following section. 
For the purpose of the frequency analysis of a following section, it is assumed 
that the dc scale factor is constant and that the worst- case drift in the time 
constant is one percent. When combined with the coil constant, the feedback 
factor becomes 


6 = K K 
r g c 

= 6. 67 x 10^ y 

which is the reciprocal of the assumed scale factor of 150 x 10" 6 volts per 
gamma. K 


" ■ ' NA S 5 23661/5 3 

FLUXGATE TUNED SYNC. I 



Figure 5-3. The single-loop system is characterized by the need 
for wide bandwidth in the critical 17-bit ADC and the current 
generator. 
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Section 5 - Electronics Design 
Subsection 5.1- Configuration Analysis 

CHARACTERISTICS OF THE TWO- LOOP SYSTEM 

In the two-loop system the inner loop provides the high resolution and band- 
width required, while the outer loop — with its offset field generator — 
reduces susceptibility to systemic errors. 

The two-loop system is characterized by an inner loop, which is 
similar to the single-loop system, and an outer loop or offset field 
generator (OFG). The block diagram of the two-loop system is shown 
in Figure 5-4. 

Inner Loop - Table 5-1 shows the overall characteristics of the inner 
loop compared to that of the single-loop system. The most significant 
parameter variation is that of the output voltage resolution, which is 
16 times greater for the two-loop system. 

Noise Considerations - The two-loop system can accept systemic noise 
i n the wide bandwidth loop of more than ten times that of the single-loop 
system. This is a distinct advantage during implementation when dealing 
with digital tran sients , power supply no ise, spacecraft interference, and all 
the other sources of noise which are always present and not usually amen- 
able to analysis. 

Outer Loop - The outer loop, or OFG, must provide the additional 
dynamic range to allow the two-loop system to operate over the specified 
range of ±64, 000Y. It accomplishes this by sampling the converted digital 
output of the inner loop to determine whether it is within either the bottom 
or top one-eighth of full scale (or in between) at a sample rate of approxi- 
mately 0.5 Hz. If within either of the outer octants, the sample pulse 
either increments or decrements a six-bit reversible counter in the OFG 
logic. If within the interior 75 percent of the ADC range, the counter 

state is not changed. , 

The output of the reversible counter is input to a six-bit DAC where the 

output is converted to a feedback current in the current generator. 

The system parameters are selected so that each step of the OFG is 
equivalent to 2500Y. Thus, after each increment, a hysteresis of approxi- 
mately 500Y exists to ensure that the occurrence of OFG step changes is 
minimized. 

Typical input/output curves for increasing and decreasing inputs are 
shown in Figure 5-5. An obvious secondary advantage of this two-loop sys- 
tem is that a cross calibration between the ADC and the OFG is provided. 

Operating Range - With the selected parameters, this system is capable 
of operating over a range of -80, 000Y to +81, 500Y; however, if desired, the 
logic could limit the up/down counter range to a limit of approximately 
±64, 000 Y. 

Summary - We have now described the general characteristics of the 
single and the two-loop systems, which show that each has its advantages 
and disadvantages. Further analysis of frequency characteristics and 
errors for a non- steady state field follow in order to select a system for 
detail design. 
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TABLE 5-1. COMPARISON OF INNER-LOOP OF TWO-LOOP 
SYSTEM WITH SINGLE- LOOP SYSTEM 


Parameter 

Two -Loop 
Inner Loop 

Single- Loop 
System 

Input Range 

0 to 4000Y 

±64, 000Y 

Output Range 

0 to 9. 6V 

±9. 6V 

Input Resolution 

±1. 0Y 

±1. 0Y 

Output Resolution 

±2. 4 x 10" 3 V 

±0. 15 x 10" 3 V 

Bandwidth 

25 Hz 

25 Hz 

Feedback Gain 

4. 167 x 10 2 Y/V 

6. 67 x 10 3 Y/V 


NAS » 23M1/I 4 



Figure 5-4. Two-Loop System. The two-loop system requires 
a 12- bit ADC with a 6- bit DAC in the outer feedback path. 



Figure 5-5. Typical Input/Output Curves for Increasing and Decreas- 
ing Inputs. The OFG has switching levels of 500 and 3500V and 
hysteresis of 500Y. 
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Section 5 - Electronics Design 
Subsection 5. 1 - Configuration Analysis 

FREQUENCY ANALYSIS OF THE SINGLE- LOOP SYSTEM AND INNER 
LOOP OF TWO- LOOP SYSTEM 

Although the single-loop system transfer function is similar to that of the 
inner loop of the two-loop system, because of a shorter integrator time 
constant the two-loop system recovers more quickly from overload 
conditions . 

The frequency analysis of the single-loop system is also applicable to 
the inner loop of the two-loop system. Its response is critically damped 
and is -3 dB at 25 Hz. The overload recovery time of the single-loop sys- 
tem is much longer than that of the two-loop system. 

The equations expressing the frequency response characteristics of the 
single-loop system are shown in Table 5-2. Equation (8) expresses the 
basic transfer function of any closed-loop system, while Equations (9) and 
(10) relate the forward and feedback gains to the specific blocks described 
in the previous topic. 

After making the notational simplifications shown. Equation (11) gives 
the exact closed-loop response of the system. At this point in the analysis 
a simplifying assumption was made based on the assumed values of the time 
constants shown. The selected value gives a ratio of this inequality of 
approximately one to 1600, which certainly justifies its use. 

The resulting Equation (12) can be factored, giving the final result 
shown in ( 1 3) . 

Analysis Methodology - It should be noted that the actual sequence of 
synthesizing the loop parameters was approximately the reverse of the 
order of the equations. The system response was selected for critical 
damping with a -3 dB frequency response at 25 Hz, giving Equation (13). 

The topology (shown previously) was selected based primarily on experi- 
ence with other similar systems. The values of K s , K a , Kj, and T a were 
set by considerations described previously. The factors p and t c were 
dictated by the required dc response and coil design, respectively. 

The two free parameters left were the integrator and the demodulator 
filter time constants. The final exercise left was to calculate these values. 

Inner Loop of Two-Loop System - The transfer function of the inner 
loop of the two-loop system is exactly the same form as that of Equation 
(13). The only difference in values is of Kf, the dc feedback gain, which, 
in turn, gives a different value of t^, the integrator time constant, since 

ra- 
the 
ng a 

feedback gain of 417Y /V for the two-loop system, its integrator time con- 
stant is 0. 036 second. 

These values indicate that the two-loop system with its much shorter 
integrator time constant will recover from overload conditions much more 
quickly than the single-loop system. 

We now analyze the errors of the described system for non-steady 
state input fields. 

) 


t x and K x are constants for both systems. 

Comparison of Loop Time Constants - Based upon the preceding pa 


graph and using values of 8. 73 x 10--* seconds for x x and 10' 2 for K x , 
calculated value of t • for the single-loop system is 0. 582 second. Usii 
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TABLE 5-2. EQUATIONS FOR EXPRESSING THE FRFnrrrMrv 
RESPONSE CHARACTERISTICS OF SINcS LOOpTySTEM 


The loop parameters are 
of Equation (13). 


selected from the desired closed-loop 


response 


A = o 

c 1 + P A 


( 8 ) 


K K K, 

A - sad 

° T ; s (1 + t s) (1 + r S ) 


P = 


= demodulator gain 
= demodulator time constant 
K g = current generator gain 
A = closed-loop gain 
A 0 = forward open-loop gain 


P - feedback gain 


K K 
g c 

(1 + r c s) 


K g = sensor gain 
K a = ac amplifier gain 
T a = ac amplifier time constant 
T £ = integrator time constant 
T c = coil time constant 
K = feedback coil gain 


(9) 


( 10 ) 


Letting: K f = K g K,, Kjf = K s K;> K d , r = 


T . 
1 


K K. 


A - -L 
c - K f 


(1 + T.S) 


1 + v +T x (T a + T d + T c )s + T *< Vd + Vc + T d + T x T a T d T c s ' 


Under the conditions that: t ^<< T (t + \ 

c x v a d ; 


A c * K 


L[ l 

f L 1 + (T x - T c )s + T x (T a + T d )s2 + T x T a T c 

K f [(1 + T lS ) 2 (1 + r 2 S>] 

Where Tj = 4. 07 x 10' 3 sec and t 2 = 0.48 x 10' 3 sec. 


( 11 ) 

( 12 ) 

(13) 
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Section 5 - Electronics Design 
Subsection 5.1- Configuration Analysis 


ERROR ANALYSIS OF THE SINGLE- LOOP SYSTEM AND 
INNER LOOP OF TWO- LOOP SYSTEM 

The errors for non- steady state fields were not specified but were consid- 
ered worthy of analysis. These errors were the same for both systems. 


Errors for non- steady state fields dominate the dc error for frequencies 
above 3 Hz. The lag for ramp inputs is also analyzed. Both errors 
are applicable to the single-loop system and the inner-loop of the two-loop 
system. 

As the frequency of the input field increses, the loop gain (pAo) 
decreases, giving rise to increasing errors due to uncertainties in the 
open-loop parameters. 

Equation (12) from the previous topic was used for the error analysis. 
The magnitude of the gain (IA| ) was calculated, from which the partial 
derivatives of Equations (14) were derived. 

dlAl f 1 [Id IAI J \2 . / 9|A[_ 


IAI 


IAI 


J^H 2+ (^: dT l + (^ dT l 


1/2 


(14) 


The RSS error was then calculated as a function of frequency. Equation (15) 
was required to relate dT x to the circuit parameters from which it was 
derived. 


dT 




ft ft 


1/2 


(15) 


Each of the partial derivatives, which are all functions of powers of w, 
were calculated for the frequencies shown in Table 5-3. The individual 
error terms were then derived and RSS*d to give the ac error terms shown 
in the table. 

The dc error term is simply the expected error in Kf, the dc component 
of the feedback gain, P . The derivation of this error is discussed in a later 
topic describing the current feedback amplifier. 

It can be seen from the data in Table 5-3 that the dc error term domi- 
nates the total error through a frequency of w = 1.0. The ac error term 
then becomes predominant and increases to about 6 percent at 25 Hz 
(w = 157) and to almost 10 percent for higher frequencies. 

The largest factors in the ac error terms are those caused by the esti- 
mated 5 percent uncertainties in the sensor and ac amplifier gains. 

Ramp Input Error - Using Equation (13) from the previous topic and let- 
ting the field input be in the form of a ramp (Bi(s) = B/s ), one can calcu- 
late the offset error for any ramp input. The result is that 

( 16 ) 


B(offset) Y = 8.6 x 10" 3 x Bi 


Using this equation, the offset error due to the maximum expected ramp 
input of 250 Y/sec is 2. 15Y. 


5-10 






time-zgoo 

i The i re ? UCn f y and ramp errors described above also apply to the inner 

loop of the two-loop system. y inner 

TABLE 5-3. PARTIAL DERIVATIVES FOR VARIOUS FREQUENCIES 
The ac error dominates for frequencies above 3.0 radians/sec 

F r equency 

Radians Per Sec AC Gain Error 

DC Gain Error 

Total RSS Error 

10" 2 

4.2 x 10" 10 

2 x 10" 5 

2. 0 x 10" 5 

10" 1 

4. 2 x 10" 8 

2 x 10' 5 

2. 0 x 10" 5 

1 . 0 

4.2 x 10" 6 

2 x 10" 5 

2. 0 x 10" 5 

3.0 

3. 8 x 10 ^ 

2 x 10" 5 

4. 3 x 10 -5 

10 

4.2 x 10' 4 

2 x 10" 5 

4. 2 x 10' 4 

30 

3. 7 x 10" 3 

2 x 10' 5 

3. 7 x 10' 3 

io 2 

3. 3 x 10" 2 

2 x 10" 5 

3. 3 x 10" 2 

157 

5. 8 x 10" 2 

2 x 10" 5 

5.8x10 2 

300 

8. 7 x 10' 2 

2 x 10" 5 

8.7x10 2 

103 

8. 4 x 10" 2 

2 x 10' 5 

8.4 x 10" 2 

3 x 10 3 

8.9 x 10~ 2 

2 x 10' 5 

8. 9 x 10" 2 

104 

9. 5 x 10' 2 

2 x 10" 5 

9. 5 x 10" 2 

4 

3 x 10 

9. 6 x 10" 2 

2 x 10 -5 

9. 6 x 10" 2 
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Section 5 - Electronics Design 

Subsection 5. 2 - Electronics Design and Performance Analysis 

SUMMARY OF TWO- LOOP SYSTEM DESIGN PARAMETERS 

The two-loop system design parameters are established to enable alloca- 
tion of errors. The estimated steady-state errors from the resulting 
analyses meet the requirements. 

The frequency stability of the two-loop system has been discussed 
along with its general topology and operation. The ability to meet the 
specified performance must now be established. 

Electronic Error Budget - Table 5-4(a) shows the overall budget of 
errors allocated to the electronics subsystem. The following topics of 
this subsection will address the sources of these errors and estimate their 
magnitude. The analysis will be performed assuming steady-state (or dc) 
conditions, since the unspecified frequency-dependent instabilities were 
addressed previously. 

Summary of the Error Analysis - The summary of the following error 
analyses is shown in Table 5-4(b). As indicated, the errors were calcu- 
lated for long-term (one year) drifts and for temperature drifts assuming 
both no control (±30°C) and for a controlled environment of ±5 C for cri- 
tical components. 

Comparison of the results versus the requirements shows that the over- 
all electronics system is capable — but just barely— of meeting the require- 
ments without a temperature controlled environment for the critical com- 
ponents. Therefore, it is strongly recommended that a thermal controller 
be used. This will provide a margin of safety and will also ensure that 
moisture is not alllowed to affect the Vishay resistors, which could cause 
an additional 20 ppm per year of long-term drift. 

Data Requirements- It might appear that the data requirements of the 
two- loop system (12 + 6 bits) are greater than that required for a single - 
loop system with the same resolution (17 bits); however, the six-bit DAC 
need only be sampled for incrementing and read out approximately once 
every two seconds, assuming a maximum input rate of 200 T/second. Thus 
if the desired sampling rate were 8 samples per second, the single-loop 
system would require 8 x 17 = 136 BPS while the two-loop system would 
need only 8x12 + 6/2 = 99 BPS. This is a significant saving in telemetry 
space and an additional argument for the two-loop system. 
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TABLE 5-4. SUMMARY OF TWO- LOOP SYSTEM ERROR BUDGET 

All error budgets can be met with temperature control except for a 
slight over condition on long-term zero offset. 


(a) Error Budget 


Source 

Long-Term 

Temperature 

Total 

Zero Offset 

0. 33 y/year 

0. 33y 

0. 47 

Linearity 

— 

o. iy 

0. 10 

Scale Factor 

2. 5 7/vear 

2 . sy 

3. 54 

RSS Total 

2. 52 y 

2. 52Y 

3. siy 


(b) Error Analysis (Errors inT ) 


Source 

Long Term 

Temperature 

±30°C 

Temperature 

±5°C 

Total 

±30°C 

T otal 
±5°C 

Zero Offset 

Demod/Int 

0. 05 

0. 05 

0. 05 (1) 

0. 07 

0. 07 (1) 

DAC 

0. 25 

0. 30 

0.05 

0. 39 

0. 26 

Current Gen . 

0. 24 

0. 27 

0. 05 

0. 36 

0. 25 

Zero Offset Total 

0. 35 

0. 41 

0.09 

0. 54 

0. 37 

Scale Factor 

Loop Non-lin. 

0. 001 

0. 001 

0. 001 (1) 

0. 002 

0. 002 (1) 

Loop Gain 

0. 008 

0. 003 

0. 003^ 

0. 009 

0. 009 (1) 

ADC 

0. 15 

0. 08 

0. 08 (1) 

0. 17 

0. 17 (1) 

DAC 

0. 9 

2. 1 

0. 35 

2. 3 

0. 97 

Current Gen 

0. 26 

2. 1 

0. 35 

2. 1 

0. 44 

Scale Factor Total 

0. 95 

2. 97 

0. 50 

3. 0 

1 . 1 

Total Accuracy 

1. 0 

3.0 

0. 51 

3. 0 

1. 2 


(1) These circuits not temperature controlled. 
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Section 5 - Electronics Design 

Subsection 5. 2 - Electronics Design and Performance Analysis 

DESIGN OF THE REFERENCE GENERATOR AND THE PROBE DRIVE 

The recommended probe drive circuit uses a crystal oscillator with a 
digital count -down chain and reference phase adjustment. The probe drive 
provides high efficiency operation. 

The probe drive circuit and reference generator were determined to be 
areas where high confidence in many proven solutions existed. The probe 
design calls for an excitation frequency of 10 kHz, which is comparable to 
many existing designs. 

Requirements - The general requirements of this subsystem include 
good frequency stability, easy phase adjustment between the probe drive 
and the demodulator, 2fo, reference signal, and an efficient probe drive 
with good symmetry to minimize even harmonics. 

A Solution - The schematic of Figure 5-6 shows one method for gener- 
ating the probe drive signal and the demodulator reference voltage which 
meets these requirements. 

As shown, a 2-MHz crystal oscillator is used to generate a stable fre- 
quency. Its output is counted down in devices which provide decoded deci- 
mal outputs. Proper selection of one of the R lines and one of the C lines 
enables setting the relative phases of drive and reference signals to about 
two degrees. The properly phased signals from the countdown chain are 
divided once more to generate the probe drive reference. 

The circuit used to generate the probe drive current is one which has 
been used with the MJS Spacecraft Magnetometer System described by 
M. Acuna*, and is frequently used to drive magnetic amplifiers. This 
circuit has the advantage of providing a large ratio of peak-to- aver age 
current resulting from the discharge of the capacitor into the core when 
it saturates; it exhibits high efficiency because the drive circuit supplies 
only circuit losses. The nonlinear resonant circuit provides energy stor- 
age. The digital drive generator has the advantage that the required odd 
symmetry of the square wave voltage is easily accomplished at these fre- 
quencies with CMOS logic elements. 

Redundancy - It is planned that each of the three axes will have iden- 
tical circuit boards, including the reference generators. Additional 
studies should be performed on the total system reliability to determine 
the best way of using the available redundancy of the three oscillators. 

Next we describe the second harmonic amplifiers of the inner loop. 


*M. H. Acuna, IEEE Transactions on Magnetics , Vol. MAG- 10, 519 (1974). 
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Figure 5-6. Reference Generator and Probe Drive. This design uses 
a crystal oscillator for stability and a high- efficiency drive circuit. 
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Section 5 - Electronics Design 

Subsection 5. 2 - Electronics Design and Performance Analysis 
DESIGN OF THE SECOND HARMONIC AMPLIFIER 

As with the probe drive circuit, the implementation of the second harmonic 
amplifier is amenable to proven solutions. 

The second harmonic signal generated in the fluxgate probe contains 
the magnetic field information and must be amplified in the presence of 
large odd harmonic components which are coupled across the core. 

Requirements - The amplifier must-be tuned to the second harmonic 
signal, providing a gain of about 4 x 10 3 and a high rejection of the funda- 
mental and third harmonics; it must not create significant noise relative 
to that generated in the probe (less than 0. 03 Y rms); and it must be 
designed so that overload conditions (which will certainly exist at turn-on) 
do not cause excessive phase shifts due to hard limiting. 

A Solution - As might be expected, a large number of solutions to 
these requirements have been implemented for previous space programs 
(including those designed at TZL). A recommended solution is shown in 
the simplified schematic of Figure 5-7. This circuit provides passive 
gain and tuning in the input circuit. Cl is selected to provide series reso- 
nance with the probe secondary winding inductance. The values shown 
provide a Q of about 5 for this resonant circuit and are large enough to 
swamp the stray cable capacity. 

The input transformer is selected to provide a voltage gain of about 20. 
This value was selected on the basis that the secondary capacitance, C2, 
should not be less than 30 pF to ensure minimum changes due to stray 
capacities, including the amplifier input capacity. The capacity reflected 
into the primary is therefore about 0.012 pF. The transformer inductance 
is tuned for parallel resonance with the total secondary capacitance. 

Having achieved a passive gain of 20, the additional gain required is 
about 200. This can easily be provided by two stages of gain with a tuned 
filter coupling them. 

Both amplifiers should contain controlled limiting characteristics to 
ensure, under overload conditions, that the phase shift through the ampli- 
fier does not become excessive. 

Noise - A reasonable estimate of noise at the input of the amplifier is 
about 30 nV / \'Hz. Assuming a bandwidth of about 600 Hz, this represents 
an input noise of approximately 0. 73 pV which, when reflected through the 
input transformer and the sensor scale factor, results in an equivalent 
magnetic field noise of somewhat less than 0.0 1Y rms, which meets the 
system requirements. 

Next we describe the demodulator and integrator which convert the 
second harmonic signal to dc and provide additional gain. 



TIME-ZERO 


( 



Ls = SENSOR INDUCTANCE (=&376x 10* 6 ) 

Rs = SENSOR RESISTANCE (^9.2 fi) 

Cs = CABLE CAPACITANCE (^300 pf) 

Cl = CAPACITANCE TO RESONAT WITH Ls (=fc0.168*;f) 
C2 = SECONDARY TUNING CAPACITOR (=s30 pf) 


Figure 5-7. Second Harmonic Amplifier. Provides passive gain 
and timing in the input circuit. 
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Section 5 - Electronics Design 

Subsection 5.2- Electronics Design and Performance Analysis 
DESIGN OF THE DEMODULATOR AND INTEGRATOR 


The demodulator synchronously converts the second harmonic output to a 
dc signal which is amplified by an operational amplifier configured as an 
integrator. 

The simplified schematic of Figure 5-8 shows the basic configuration 
of a demodulator and integrator which meet the system requirements. 

Demodulator Configuration - The demodulator configuration shown, 
although utilizing only one switch, has a gain given by 


e dc ^ R-, 

o _ Z 2 

e. peak 77 X R, + R 0 

l 12 


under the conditions that 


R 


1 


» t 


and 


R 2 c 2 


» t 


( 1 ) 


where t = — 


Thus, it provides the same output voltage as most full-wave demodulators. 
Many switches are available, including the Harris HI 5049, which has pre- 
viously been successfully used in this type of application. 

Integrator - Following the demodulator is an integrator which is mech- 
anized with an LM- 108 operational amplifier. This circuit provides the 
high loop gain necessary to ensure minimal non-linearity and scale factor 
drift. 

The transfer function of the integrator is 


where 



T. 

1 


T 

a 



( 2 ) 


= DC gain of the op-amp ( ~ 3 x 10^) 

= Integrator time constant ( 5:5 0. 036 sec) 

= Corner frequency of the operational amplifier compensation 
(= 0.2 sec) 


As indicated, the transfer function is really that of a dc amplifier with 
two poles. One pole is due to the external components, R 3 and C 3 , and 
the other pole is due to the amplifier compensation. 

Since has a nominal value of 3 x 10^, using the time constants 
shown, the corner frequencies equivalent to the poles are about 1. 5 x 
10“5 Hz for Kj t ^ and 2 x 10^ Hz for T a /K c j. The approximation shown 
is therefore valid for any analysis except for dc stability. 

Errors - The required 1. 0 x 10“^ Vdc stability at the input of the 
demodulator and integrator is easily achievable with the selected 
components. 
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An assumed short term drift in K d of ±50 percent and a long term drift of 
-300 percent give rise to scale factor errors of 3 x 10' 3 Y and 8 x 10' 3 Y, 
respectively. Assuming a non-linearity of 20 percent in K d results in an 
error (for the minimum expected value of K d ) of less than 1 x 10' 3 Y in the 
system closed-loop gain* 

In summary, the integrator and demodulator mechanization is straight- 
forward. The errors produced are negligible compared to other sources 
such as the feedback current generator* The output of the integrator is con- 
verted to a digital signal by the ADC, which is described in the following 
topic. 6 



Figure 5-8. Demodulator- Integrator Design. This simple, 
straightforward design meets all system requirements. 
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Section 5 - Electronics Design 

Subsection 5.2- Electronics Design and Performance Analysis 
DESIGN OF THE 12- BIT ADC 

The ADC design is a dual-slope integrating type. This is a proven tech- 
nique which results in small error contributions and also eliminates offset 
errors. 


The ADC converts the analog output from the loop integrator to a digi- 
tal signal. The overall system was designed to simplify this circuit by 
necessitating that only positive inputs from zero to 10 V be converted. 
Twelve-bit conversion was selected so that the digital output would contain 
a fair measure of the system noise to allow its evaluation. This approach, 
of course, also reduces the resolution error from the required ±1. 0Y to 
half that value. 

Recommended Solution - The schematic shown in Figure 5-9 shows the 
recommended implementation of the ADC. This system is a dual- slope 
integrating type ADC with the feature that an additional cycle is added 
which measures the system offset and subtracts it from the digital output 
data. 

Four subcycles of the system operation are used — each requires a time 
equivalent to 4096 clock periods, where the clock period is approximately 
2 ns. The conversion thus requires about 32 ms. 

At the beginning of the first period, S4 is opened and SI closed. The 
integrator output, at the end of 8 ms, is proportional to the offset voltage, 
Vq|. This offset is adjusted to ensure that the polarity of the offset is 
equivalent to a positive input. 

During the second period, SI is opened and S2 is closed, applying the 
reference voltage to the input. Clock pulses are then counted into the 
logic system until the comparator fires. This count, of Equation (1), 
is a measure of the total offset voltage. 

During the third period, S3 is closed, thus integrating the total value 
of (V^ + Vqi) for 4096 counts. 

Finally, in the fourth period, S2 is again closed while clock pulses are 
counted in the logic until the comparator threshold is again reached (note 
that the comparator fires twice — both times in the same direction, elimi- 
nating hysteresis and timing problems). The count stored in N 4 is shown 
in Equation (2). 

The logic subtracts N 2 from N 4 , giving a total count N^, as shown in 
Equation (3). The two important characteristics of this equation are that 
the offset voltage of the comparator, Vq£, has been eliminated from the 
output count and that the integrator offset has been converted to a scale 
factor error. 

System Errors - Table 5-6 shows an estimate of the system errors, 
including those of the reference voltage (which is described in the DAC 
section) and the integrator. The total RSS error expected for both tem- 
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TABLE 5-5. DETERMINATION OF OUTPUT COUNT 

The equations show that the comparator offset voltage is cancelled and the 
integrator offset causes a scale factor error. 


V 01 N 

T 

v oz 

( - V R + V 01> 

At 

( - V R + V 01» 

( V i+ V 01 )N 

T 

V 02 

( - V R + V 01> 

At 

<- V R + V 01> 


N ( = N 4 


<- V R + V 01> 


N^ = Counts for offset 
N^ = Counts for input plus offset 

N = Output count 

N = Full scale count (4096) 

At = Converter clock period 


V. = Input voltage 
V Q1 = Integrator Offset 
Vq 2 = Comparator offset 
V R = Reference voltage 

t = Integrator time constant 


TABLE 5-6. ADC SYSTEM ERRORS 
The ADC errors are minor with no temperature control. 
(a) | (b) ~~ 


Source 


Scale Factor Error 
PPM/°C PPM/yr 


Voltage Ref 

1. 0 

12. 5 

Integrator 

V 

os 

0 . 1 

2. 0 

Integrator 

i 

o s 

0.4 

12. 0 

Sw Res 

Y 

o s 

0. 5 

10. 0 

Sw Time 

t 

sw 

0. 01 

0.2 

RSS Total 
(PPM) 


1.2 

20. 1 

RSS Y 


0. 005 Y 

0. 08 


Error Source 

Scale Factor Temp 

Scale Factor Long 
Term 

RSS Total 


Error Y 

No Temp 
Control ±30°C 

0. 15 
0. 08 

0. 17 V 
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Section 5 - Electronics Design 

Subsection 5. 2 - Electronics Design and Performance Analysis 
DESIGN OF THE FEEDBACK CURRENT GENERATOR 

The outputs of both inner and outer loops are summed and converted to a 
current for the feedback coil. Temperature control of the critical compo- 
nents is required. - 

The magnetometer system stability is directly proportional to the 
stability of those elements in the feedback loop which include the current 
generator. A simplified schematic of this circuit is shown in Figure 5-10. 

The system consists of a two-stage voltage summing amplifier which 
accepts the outputs of the inner and outer loops and scales them properly, 
and a two- stage current generator which converts its input voltage to a cur- 
rent through a scale factor equal to the resistor, Rl. 

Two stages of gain were used for both the summing amplifier and the 
current generator to ensure negligible errors due to gain drift and non- 
linearities in the operational-amplifiers. Precision Monolithic Op-07 ampli- 
fiers are used in both input stages due to their extreme stability. This 
device was also used in the current generator output stage due to its cur- 
rent output capability (8 mA full-scale current is assumed). 

The precision resistors required are Rl through R4. These would be 
mechanized using Vishay precision film components, which experience has 
shown to be highly reliable and stable. 

Standard non-precision resistor-capacitor networks are used to stabil- 
ize both loops, with one series network used across the coil. The only 
significant pole in this circuit is that due to the L/R feedback coil constant, 
which is 0. 11 x 10" 3 second. 

Component Instabilities - The component instabilities which were util- 
ized in the error analysis are shown in Table 5-7. 

Error Analysis - Using the value of Table 5-7 an analysis was performed 
to determine offset and scale factor errors, which are temperature depen- 
dent both with and without a temperature controlled environment, and long- 
term errors. The results of this analysis are shown in Table 5-8(a) and (b). 

As indicated, the total inaccuracy of the current generator, both long 
and short term, for scale factor and zero offset, are reduced by a factor 
of four if the critical components are kept in a temperature controlled 
environment of i5^C* This is highly recommended to provide a margin of 
error and also to eliminate moisture from entering the components if the 
thermal environment is hermetically sealed and evacuated. 

The other critical element in the feedback path is the DAC of the outer 
loop, which is described next. 



Figure 5-10. Feedback Current Generator. The circuit sums the voltage 
outputs of both loops and converts the sum to a feedback current. 
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TABLE 5-7. COMPONENT INSTABILITIES USED IN ERROR ANALYSIS 


Instability- 


Source 
OP -07 

Offset Voltage 
Offset Current 
Gain Error 

Vishay Resistor 
T racking 
Absolute 


Temperature 

0. 6 MV/°C 
25 P A/°C 
10% over ±30°C 

0. 5 ppm/°C 

1 . 0 ppm /°C 


Long-Term 

12 liV /year 
700 pA/year (1) 
50%/year 

. / ( 2 ) 

3 ppm/year 

3 ppm/ year ^ 


(1) This was based on an assumed doubling of the offset current 
over a 5-year period, 

(2) This assumes that the resistors will be hermetically sealed 
to eliminate a possible 25 ppm/year instability due to 
moisture. 

TABLE 5-8. ANALYSIS OF OFFSET AND SCALE FACTOR ERRORS. 
The DAC errors are significantly reduced by a crude temperature 
controller. 


(a) 


Scale F actor 

Zero Offset 

Source 

ppm/°C ppm/yr 

ppm/°C 

ppm/yr 

DAC Input R4/R2 

0.47 2.8 



F. B. Input R3/R2 

0.03 0.2 



S. F. Res. R 1 

1.0 3.0 



Zj Vos 


0. 2 

3. 0 

Z j Ios 


0.01 

0. 28 

Z y Z Gain 

0,007 0.5 


_ 

Z 3 Vos 


0. 12 

2. 4 

Zj Ios 


0. 001 

0. 03 

Z^ Gain 

0, 001 0. 15 


- 

RSS Total (ppm) 

1.1 4. 1 

0. 23 

3. 9 

RSS Total? 

0.07? 0.26? 

0. 01? 

0. 25? 

(b) 


Errors (?) 


Error Source 

No Temperature 

Temperature 

Control 

Control ±30°C 

±5°C 


Scale Factor Temperature 2. 1 

0. 35 


Scale Factor Long-Term 0.26 

0. 26 



RSS 2. 1 

0. 44 


Zero Offset Temperature 0. 30 

0. 05 


Zero Offset Long-Te 

rm 0. 25 

0. 25 



RSS 0. 39 

0. 26 


Total Accuracy RSS 

2. 14? 

0. 51? 
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5. Z - Electronics Design and Performance Analysis 


DESIGN OF THE DIGIT AL-TO- ANALOG CONVERTER (DAC) 


As with the current generator of the previous section, the DAC is also in 
the feedback path and, therefore, must exhibit high stability. Simple tem- 
perature control of critical elements en hances performance. 


The DAC accepts inputs from the Up /Down counter in the OFG logic 
and converts the 6-bit counter state to a corresponding precision voltage. 
The DAC output is then fed back to the feedback coil of the sensor roug 
the previously described current generator. Each step is equivalent to 


2500V 

Several DAC configurations were examined (see the next topic) and the 
one recommended is shown in the simplified schematic of Figure 5-11. 
consists of two precision reference supplies (one of each polarity), a pr 
cision R-2R network, precision switches, and an output buffer amplifier. 
The principal reasons for selection of this topology are simplicity, sym- 
metry, and the availability of suitable components. 

Comp onent Instabilities - As with the previously described current 
source, Op-07 ampliiiers and Vishay resistors are used. The previous 
table showing the instability of these components is applicable. In addi- 
tion, Harris HI 5049 switches and National LM-199 voltage references are 
used. Their instabilities are shown in Table 5-9 below. 

Error Analysis - Using the value of Table 5-9, an analysis was^per- 
formed to determine their effect on the subsystem performance. The 
results of this analysis are shown in Table 5-10 (a) and (b). As with th 
current generator, the results of temperature control are significant. 
Again it is recommended that this be used to maximize the probability of 

a successful mission. , .. , ■ . 

The DAC described above is only one of several configurations w 
were investigated. Two alternative methods of providing this function are 
described in the next topic. 

NAS 5 23661 5 1 1 



Figure 5-11. Digital-to-Analog Converter (DAC). The DAC configuration 
uses precision R-2R network and two Reference Supplies. 


) 
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TABLE 5-9. COMPONENT INSTABILITIES USED IN ERROR ANALYSIS. 
The instabilities of the additional components used in the DAC are shown. 


Instability 


Source 


Temperature 


Long Term 


Switch 

ON resistance 0.05 57/ C 1.0 £2/year 

OFF leakage 1. 8 na/°C 31 na/year 

Zener voltage 0. 5 ppm/°C 12 ppm/year 


TABLE 5-10. ANALYSIS OF OFFSET AND SCALE FACTOR ERRORS. 
The error analysis shows that temperature control reduces errors 
about in half. * v 




Scale Factor 

Zero 

Offset 

Source 


ppm/°C 

ppm/yr 

ppm/°C 

ppm/ yr 

Switch on Res. 

r DS 

0. 29 

5.8 

— 

— 

Switch Leakage 

*D off 

0. 07 

0.9 

— 

— 

R/ 2R Network 


0. 29 

1. 73 

— 

— 

z i 

V 

os 

— 

— 

0. 06 

1. 2 

z i 

i 

os 

— 

— 

0. 125 

3. 5 

Voltage Ref 

V R 

1. o 

12. 5 

— 

— 

RSS Total 

1. 1 

13. 9 

0. 14 

3. 7 

RSS Total 7 


o. 07 y 

0. 89Y 

o. 009 y 

0. 24 Y 


(b) 

Errors (y) 

No Temperature Temperature Control 
Error Source Control ±30°C ±5°C 


Scale Factor Temperature 

2. 1 

0. 35 

Scale Factor Long Term 

0.9 

0. 9 

RSS 

2. 3 

0. 97 

Zero Offset Temperature 

0. 27 

0. 05 

Zero Offset Long Term 

0. 24 

0. 24 

RSS 

0. 36 

0. 25 

Total Accuracy RSS 

2. 3 y 

1. 0Y 
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Section 5 - Electronics Design 

Subsection 5. 2 - Electronics Design and Performance Analysis 
ALTERNATIVE DAC CONFIGURATIONS 

Two alternative methods of implementing the DAC were investigated and 
are worthy of consideration for use — switched coils and pulse-width 
modulation. 


Switched Coils - One configuration uses a constant current source which 
is passed through multiple coils, as shown in Figure 5-12. The output of the 
UP/DOWN counter is used as the control inputs for the coil switches via a 
ROM code converter. The key to this implementation is that the current 
can be passed through each coil in either direction or may bypass the coil 
completely. This ternary type system can be implemented with four coils 
of ratios 1, 3, 9, 9 to provide plus and minus 22 steps of OFG flux. 

The switches can be mechanized by devices such as the Harris HI5049. 
The ON resistance of these devices is nominally 30 ohms, for a total of 
240 ohms for the eight switches. A preliminary error analysis indicates 
that this resistance is non- critical, as in the leakage current of the OFF 
switch sections. 

The primary disadvantages of this system, which precluded its recom- 
mendation for use, is that the number of wires required complicates the 
boom cable and also complicates the fabrication of the feedback coils. 

One method of easing the cabling problem would be to transfer serial data 
down the boom to a non-magnetic logic hybrid in the sensor assembly which 
would perform the required storage and decoding. 

Pulse-Width Modulator - A second configuration which has merit is a 
DAC which utilizes pulse-width modulation. Figure 5-13 shows a simplified 
schematic. 

This system generates a dc pulse of the proper polarity with constant 
amplitude and a width proportional to the contents of the UP /DOWN counter 
at a sample rate of about 700 Hz. This frequency was selected as the geo- 
metric mean of the loop bandwidth (25 Hz) and the 2fo frequency (20 kHz) to 
minimize the effect of ripple on system operation. 

The errors for this system are very comparable to the recommended 
DAC with the additional error of about 1. 0 ppm/°C due to drift in the switch- 
ing transients. The primary reason, however, for not recommending this 
method is lack of experience with its operation. 

Next we describe the auxiliary circuits, which are not part of the basic 
magnetometer but would be required in a full system implementation. 
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Figure 5-12. Constant Current Source. The current is applied 
through switched coil segments. 
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Figure 5-13. This DAC is basically a pulse-width modulato: 




Section 5 - Electronics Design 

Subsection 5.2- Electronics Design and Performance Analysis 
AUXILIARY SUBSYSTEMS AND POWER ANALYSIS 

The basic magnetometer studied has been described. To perform overall 
weight and power estimates, the auxiliary s ubsystems must be included. 

The simplified block diagram of Figure 5-14 shows the auxiliary cir- 
cuits which will be required to mechanize the magnetometer system. They 
include the digital subsystems, interface circuits, a low voltage power sup- 
ply, and temperature controllers. Each of these circuits has been investi- 
gated in sufficient detail to enable realistic power estimates and to allow an 

estimate of a typical mechanical configuration. 

Logic and Telemetry Subsystem - The logic and telemetry blocks were 
sized assuming the use of CMOS logic elements throughout It was 
assumed that the telemetry interfaces would be typical of other spacecraf 
and that the data would not have to be double buffered. Elements for sys- 
tern initialization were included. A total of 50 CMOS packages are required. 

Temperature Controlling - One temperature controller is required for 
the sensor and is described elsewhere. An additional temperature con- 
troller for the critical components of the DACs and the current generators 
will require an estimated 0. 7 W of power. The temperature control 
requirements are easily met since ±5°C is sufficient. An operating tem- 
perature of +85°C is assumed with a basic construction consisting of dis 
Crete elements mounted on a hybrid substrate which, in turn, is mounted 
in a hermetically sealed, evacuated housing to minimize heat leaks 

Low Voltage Power Supply - Many designs are available to meet the ^ 
LVPS requirements. The LVPS should use switching techniques to provide 
efficiencies in the order of 75 percent. Post regulators are assumed on 
±15 V analog lines to provide better than 0. 5 percent regulation. 

' Power Estimate - Table 5-7 shows the estimated power breakdown. 
These estimates are considered realistic and, as indicated, enable the 
system to conform to the specified 6 W limit. 

We next describe a method of packaging the entire system and give 

results of a weight analysis. 
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TABLE 5-11. POWER ESTIMATE. 
The power estimate meets the 


specification. 

Subs ystem Power (mW) 

Sensor Drive (3) 216 

AC Amplifier (3) 60 

Demodulator (3) 30 

Integrator (3) 30 

ADC (12-Bit) (3) 525 

DAC (6-Bit) (3) 405 

Current Generator (3) 705 

Reference Supplies (3) 195 

Logic and Interface (1) 352 

Temperature Controller (2) 1000 

3518 

Low Voltage Supply at 75% Efficiency 117 3 

4691 


NAS 5 23061/5 14 



Figure 5-14. Auxiliary Circuits. Auxiliary circuits were included in 
weight and power analysis. 
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Subsection 5. 2 - Electronics Design and Performance 

PACKAGING AND THERMAL CONTROL OF THE ELECTRONICS 

A one-piece aluminum housing with a top access cover provides a struc- 
turally rigid, lightweight housing with good heat transfer characteristics. 

The electronics must be designed to withstand shock and vibration and 
operate in a vacuum. The electronics housing- must be lightweight, use 
space approved materials, and meet electrostatic and RFI requirements. 

Structural Rigidity - The electronics housing is a high strength, alumi- 
num alloy structure with a resonant frequency at least one octave higher 
than the circuit cards. The mounting of the circuit cards provides addi- 
tional stiffening to the chassis. Resonant frequency of the circuit cards is 
about 200 Hz, with a Q of 14. With a 5 jg input, this results in a maximum 
stress in a card equal to 1200 pounds per square inch, which is well below 
the allowable yield stress of 35, 000 pounds per square inch for epoxy glass 
material. 

Heat Transfer - Heat transfer within the electronics housing is by a 
combination of radiation (between circuit cards and between circuit cards 
and the aluminum case) and conduction (through the circuit card mounts 
and the aluminum case). The principal method of heat transfer from the 
electronics housing is through the aluminum case to the spacecraft shelf. 
With normal spacecraft surface specifications (63 microinch finish and 
0.001 inch per inch flatness), a thermal resistance between the aluminum 
case and the spacecraft shelf of Ri = 10 o C-inch^/W can be expected. With 
an interface area equal to 30 square inches and a maximum dissipation of 
5 W, the spacecraft to aluminum case temperature rise will be about 2°C. 

Material Cleanliness Requirements - Materials selected for use will 
meet the low out-gassing requirements of 1% total weight loss and 0. 1% 
volatile condensibles. All materials, with the exception of certain com- 
ponent leads and magnetics, have magnetic permeabilities of less than 1.02. 

Electrostatic and RFI - The aluminum case will be chromate conversion 
coated, providing an electrically conductive mounting surface. The cover 
will fit in a groove in the case to provide a lap joint with a long path to 
avoid radiated noise. An internal RFI filter is used at the spacecraft inter- 
face connector for conducted noise. 

Acces sibility - Assembly of each circuit element into the case is facil- 
itated by use of an internal harness. Because each circuit card size is 
identical, test setup is simplified. 

Size a nd Weight - Overall housing size is 6 inches (15. 24 cm) long by 
5 inches (12. 7 cm) wide by 3 inches (7. 62 cm) high. Total weight is 
2. 25 pounds (1. 02 Kg). 

The electronics housing is designed to meet the requirements normally 
imposed on space flight hardware. 
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TABLE 5-12. DESIGN REQUIREMENT AND APPROACHES 
Requirement Approach 

Structural Rigidity _ High cross section modulus. 

- High strength-to-weight ratio 
materials. 

- Additional support from circuit 
cards. 


- Good thermal path from dissipating 
elements to heat sink. 

- Low thermal resistance from instru 
ment base to spacecraft shelf. 

- Low outgassing. 

- Low magnetic permeability. 

- Electrically conducting case. 

- Line filters. 

- RFI-tight joints. 



Good Heat Transfer 


Material Cleanliness 


Electrostatic and RFI 
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Section 6 - Calibration and Alignment 


NATURE OF THE CALIBRATION AND ALIGNMENT PROBLEM 

The total specified overall accuracy of 5Y per vector field component for 
this instrument requires a stringent error budget limiting the contribution 
of each of the many potential sources of error. With individual sources 
necessarily small, their accurate determination becomes difficult. 

The projected performance of this magnetometer substantially exceeds 
that of previous designs in stability and accuracy of alignment of the mag- 
netic axes. This advance in the state of the art in magnetometers will 
require a commensurate advance in the techniques and instrumentation 
used to characterize and calibrate such sensors. The degree of precision 
sought in the determination of sensor parameters is in many cases as 
great or greater than that of the available measuring instruments. 

The table opposite enumerates the primary sensor parameters that 
the calibration program must determine and the allowable errors attributed 
to each of them expressed in units of magnetic field. 

Offset - The zero applied magnetic field output of the sensor is called 
offset and its maximum value and drift per year is specified at 0. 47. Due 
to the presence of the earth's magnetic field, measurement of this param- 
eter requires the generation of an opposing vector field to obtain a field- 
free volume near the sensor. Since any such cancellation procedure is 
likely to be imperfect, a gimbal mount for the sensor is also required to 
allow 180-degree rotations to remove the effects of any residual fields on 
the measurement. 

Noise - The random fluctuations in output from a sensor in a steady 
field can readily be measured for a given bandwidth provided that a suffi- 
ciently quiet field environment can be supplied. Measurement of sensor 
noise of 0. 37 will require an environment that is an order of magnitude 
(0.037) quieter. The earth's magnetic field variations range in amplitude 
from a few 7 at periods of 10-50 seconds to £17 at periods of 2-8 seconds. 
A mu-metal flux tank can reduce ambient fields by several orders of mag- 
nitude; since variations scale the same way as amplitudes, adequately 
quiet fields can be achieved over periods of time of the order of 

one minute. . 

Scale Factor - The magnitude of voltage output due to a magnetic field 
applied parallel to a given sensor axis is subject to variation as a function 
of sensor operating temperature. The allowable error due to both initial 
uncertainty and stability over the operating environment is budgeted at 
2.67. Since the scale factor determination involves knowledge of an output 
voltage and a magnetic field, both quantities must be measured to the 
required accuracy. For a maximum applied field of 6. 4 x 10 7 the bud- 
geted error represents 20 ppm which, equally distributed, allows 10 ppm 
errors from each source. Currently achievable accuracies in magnetic 
field measurement at GSFC are in the neighborhood of 16 ppm, and avail- 
able potentiometric voltmeters have absolute accuracies of ~20 ppm. Sta- 
tistical averaging of many readings will be needed to improve the precision 
of the measurements. 

Linearity - In addition to varying with ambient temperature, the scale 
factor for each axis can vary as a function of the applied magnetic field 
magnitude. Such variations are budgeted to less than 2. 67, the same as 
the accuracy specified at a given field, and the same considerations apply 
to this parameter. They can best be determined simultaneously with the 
magnitude of the scale factor by a least squares fit of data to a straight 
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line. The slope is then the scale factor and the deviation from it is the non- 
linearity for each magnetometer axis. 

Axis Stability - The position of each magnetic axis may vary due to 
thermal and/or mechanical effects. The cone of uncertainty allotted to each 
axis is 5 arc-seconds, which, at maximum applied field, translates into an 
allowable error of ±1. 6'V. In order to determine changes of orientation less 
than 5 arc-seconds, the precision of field and voltage measurements must be 
better than the 10 ppm specified for the scale factor determination. In order 
to eliminate the effects of changes in sensor output due to rotations of the 
complete sensor assembly relative to the applied field, the changes in orien- 
tation of the whole assembly must be monitored to the same or better preci- 
sion as that desired in the final result. A system consisting of fixed autocol- 
limators (theodolites) and optical reflectors mounted on the sensor assembly 
can fulfill this requirement. 

Axis Orientation - In addition to measuring changes in orientation, the 
absolute position of at least two of the magnetometer axes must be defined to 
within one arc-second with respect to an optical reflector on the sensor pack- 
age, This extremely stringent specification requires not only great accuracy 
of magnetic field and voltage measurements, but also knowledge of the abso- 
lute position of the test facility field axes to the same degree of accuracy 
(better than 1 arc-second). When last measured several years ago, the posi- 
tion of the GSFC test facility axes were determined to an accuracy of ~10 arc- 
seconds. It may be possible to improve on this measurement with a precision 
two- axis gimbal sensor mount. Such mounts presumably exist in high preci- 
sion theodolites, but whether they can be modified to properly mount the mag- 
netometer sensor is not known (they are likely to be made of ferrous metals). 

The philosophy of defining separable sources of error and specifying the 
contribution of each to the desired total seems most likely of assuring repro- 
ducible performance within the overall requirements. Representing the 
response of the magnetometer to applied magnetic fields as a vector-matrix 
equation 

V = S B + V 

o 

the various contributions to the matrix S must be identified and separated. 
These contributions are discussed next. 


TABLE 6-1. BUDGETED ERRORS 


Parameters to be Determined 

Budgeted Error 

Offset 

0.4Y 

Noise 

0. 3Y 

Scale Factor 

2. 6Y 

Linearity 

2. 6Y 

Axis Stability 

1. 6Y 

Axis Orientation 

0. 3 Y 


4. 07 Sensor 
Total 

A-D Converter 

3. oy 


5. 07 Overall 
Total 


6-3 



Section 6 - Calibration and Alignment 


IDENTIFICATION OF CONTRIBUTIONS TO RESPONSE MATRIX 

The sensor response matrix operator can be decomposed into a sequence 
of linear transformations consisting of axis distortions and rotations. 

From these, the sensor scale factors and orientation can be determined. 

The three output voltage components of the magnetometer can be con- 
sidered to form a vector response to an applied magnetic field vector. For 
each sensor axis the output minus the offset voltage is proportional to the 
field component directed along that axis, with the "constant 11 of proportion- 
ality equal to the corresponding scale factor. That is V. - V. = X. (B • 1) 

1 10 1 

where i is a unit length vector directed along the i-th sensor axis and X. is 

l 

the i-th axis scale factor. Expressed in vector form, this is the instrument 
response equation shown in the table; the response matrix S contains the 
scale factors as well as the direction cosines describing the relative orien- 
tation of the sensor axes in the field. From measurements of sensor output 
*♦> 

voltage V and knowledge of a pplied field B, S may be determined, but the 
various contributions to S must be identified and separated since the scale 
factors and axis orientations must be individually determined. The contri- 
butions are enumerated opposite. 

Coordinate Systems - For the instrument response equation to be mean- 
ingful, vector components on both sides of the equation must refer to the 
same coordinate system. The choice of coordinates is arbitrary, but ulti- 
mately the calibration program must allow the determination of magnetic 
field components in the orthogonal "mirror" coordinate system from mea- 
sured sensor output voltages. The "mirror" coordinate system is defined 
by sensor mounted mirrors which form a part of the optical attitude trans- 
fer system that will monitor boom twist and deflection in flight. Taking V 
to be composed of the three measured sensor output voltages and B the 
three known field components applied along the test facility axes, the 
matrix S is composed of a coordinate transformation from field axes to 
sensor axes followed by a scale factor transformation. The decomposition 
of the coordinate transform depends on the orthogonality and orientation of 
the two sets of axes. 

Scale Factor Matrix - If sensor and field axes are both orthogonal and 
aligned with each other, then (V - V ) = E B where E is a diagonal matrix 
whose elements are the three scale factors i = 1, 2, 3 for the three 
axes. 
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Distortion Matrices - Neither the sensor axes nor the test facility field 
axes are necessarily orthogonal. If we denote the matrix operator that trans- 
forms an orthogonal set of axes into a distorted one by D, then the inverse 
transform, D , transforms non- orthogonal axes into orthogonal ones. With 
the matrix for^the sensor axes and field axes denoted by D y and D , respec- 
tively, D y (V - V q ) = E D b 1 B if the resulting orthogonal axes are aligned. 
The matrices D contain terms d. . = i • j where T and J are vectors directed 
along the i-th and j-th axis, respectively. Since vector dot products are com- 
mutative, d.j = d„ and the matrices are symmetric. If the distortions are 
small, the matrices are also nearly diagonal. 

Ro tation Matrices - The orthogonalized sensor axes are in general not 
aligned with the orthogonalized field axes. Representing the transform that 
expresses field components in the rotated orthogonal sensor coordinates by 

R BV’ D v (y ' y “ E R BV D B B * Multiplying both sides of this equa- 
tion by D y , (V -^ V Q ) - D yE R bv °B B and the response matrix S is equal 
to D y ER By . This is a sequence of linear transformations that repre- 

sents the transfer of applied fields into measured voltages. Note that matrix 
multiplication is not in general commutative, XY 4 YX, and the order of the 
sequence is therefore important. 

The rotation from field to sensor axes can be decomposed into a rotation 

from field axes to mirror axes R BM followed by a rotation from mirror axes 
to sensor axes R^, R^ = R^. In terms of these rotationg 

R <V nVl 8 fV RbM Db ’‘® a " d R BM V* S * R MV E '' V 1 - V Q ). 

U K BM U B B B M t ^ le ma g net ic field in components of the orthogonal 
mirror coordinate system. The calibration equation 


b m r mv 


E D y 1 (V - V Q ) is then what allows the in-flight determina- 


tion of orthogonal field components from telemetered sensor output voltages 

and the quantity R My ’ E' 1 D/ 1 = (D y E R^)' 1 must be experimentally 
determined. 

D iscussion - The matrix S is in effect a measured parameter whose nine 
elements may be determined as described later. Assuming that the distortion 
matrices Dg and D y have been independently determined, then ER = D “ * 

S D B = T expresses the unknown scale factor matrix and sensor oSLtation 
in terms of known quantities. Because rotation matrices are orthonormal 
operators, their determinant is equal to one and the sum of the squares of 
the elements of each row or column is also one. Since E is a diagonal matrix 
of V s , R J us t multiplies the i-th row of Rg y by its i-th diagonal element \ 
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Section 6 - Calibration and Alignment 

IDENTIFICATION OF CONTRIBUTIONS TO RESPONSE MATRIX (Continued) 

l 

2 

The sum of the squares of the i-th row of the matrix T is therefore X.. , 
the i-th scale factor, and dividing each row of T by its corresponding \ i 
leaves just Rgy. 

This procedure obtains E and Rgy under the stated assumptions. The 
rotation matrix R My must also be determined for the calibration equation^ 

If R fiM can be independently measured, then R M y = R BM R gy = R BM 

R where T denotes transpose and the equation is a result of the ortho- 
BV , 

normal properties of rotation matrices. Thus the scale factors and orien- 
tation of the sensor are determined as desired when S, D fi , D y and 
are determined. We next consider the determination of the response 

matrix S. 


\ 


) 
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TABLE 6-2. SUMMARY OF CONTRIBUTIONS TO RESPONSE MATRIX 
Instru ment Response Equation 


(V - V ) = s B 

v is the measured 3 component output voltage 
V q is the voltage offset vector 
B is the 3 component applied field vector 
S is a 3 x 3 instrument response matrix 

Contributions to Response Mat-riv S 


Matrix Operator 
E 



Source 

different scale factor X. for each axis 
B axes not orthogonal 
V axes not orthogonal 


r mv 

r bm 


Result: 


sensor axes rotated with respect to 
mirror axes 

field axes rotated with respect to 
mirror axes 


s = D V E R BV V 1 » ° V E *M V R bm V* 

Calibration Equation : 

B M = (D V E * (V - V o ) = G (V - V o ) 

where § M i, magnetic field .pressed in orthogonal mirror 


G = WbmV 1 ' 5 ’ 1 


Scale Factor Determination 


er bv ' D v 


- 1 


sd b = T = [tj 


1J 
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Section 6 - Calibration and Alignment 


DETERMINATION OF THE RESPONSE MATRIX 

The response matrix S and its variation with field amplitude can be deter- 
mined as a least squares error estimate from a large number of sensor 
voltage and magnetic field measurements over the ±6.4 x 10^ V range. The 
precision of the estimate improves with an increasing number of 
measurements. 

The precision with which the response matrix S can be determined is a 
limiting factor in the precision of the scale factors and orientations derived 
from it. A statistical approach is most likely to yield the best estimate of 
the true value as well as allowing an estimate of the precision of the deter- 
mination from the computed rms deviations. The inverse of the response 
matrix, S , is also needed in the calibration equation referred to earlier. 
Since inversion of matrices whose elements are measured quantities is sub- 
ject to potentially large errors, obtaining the inverse directly from the sta- 
tistical analysis is a better procedure. 

We now consider a linear least squares error fit of the measured data 
to the instrument response equation. The following assumes that all sen- 
sor output voltages have been corrected by subtraction of the previously 
determined offset voltages. This is not strictly necessary since they could 
be determined as part of the fitting procedure if enough measurements are 
made, but the independent measurement is likely to be more precise and the 
discussion is simplified under the assumption that the corrections have been 
made. 

Measurement Vectors - For each sensor axis k = 1, 2, 3 the output 
can be measured and expressed in terms of the three applied field compo- 
nents as indicated opposite. The unknown are the elements of the k-th 
row of the response matrix S which is to be determined. If three such 
measurements are made for each sensor axis with different applied fields, 
a system of three equations in three unknowns results for each axis. In 
vector form this can be represented by = BS^ k = 1, 2, 3. The ele- 
ments of the measurement vectors and the rows of the 3x3 matrix B 
now refer to the three different measurements. Assuming the matrix B is 

■+. _ i -► -*■ 

non-singular, the exact solution is = B where is the k-th row 

vector of the matrix S which is then determined. If a fourth measurement 

is made for each axis, its equation will be inconsistent with the other three 

due to errors in all four measurements. In general, for n measurements 

4 

(n > 3) at fields over the full ±6.4 x 10 V range, each will have some error 


6-8 



TIME -ZERO 


and therefore deviate slightly from the best estimate. This can be repre- 

*4 > -4 

sented in vector form as B S k - = A k for each axis where the measurement 

vectors and the deviation vectors have n components and B is a 3 x n 
matrix. 

Least Squares Error - The best estimate of error is the computed devia- 

tion. We wish to find that estimate S. of the true value that minimizes the 

k 

sum of the squares of the deviations. In matrix form, we minimize 


(VM- 


4 nr t -4 “► T T "► "► “► T 

S. B B S, - S. B V. - V. BS, + V. V. 
k kk k kkkk 


where T denotes transpose and the minimization is with respect to S^. 
is we set the partial derivative to zero: 
d 


dS, 


(K T K) 


T -*■ T ** 

= 0 = 2 B BS. - 2 B V. 

k k 


That 


This equation has a unique solution if B B is non-singular 

*4 T . 1 T 

S k = (B 1 B) B 1 V k 


which is the desired least squares estimate for the k-th axis and the best esti- 
mate for the response matrix S has therefore been found in terms of the mea- 
sured fields and voltages. 

Inverse Estimate - The same data used to generate the response matrix S 
can also be used to generate its inverse Q = S This inverse is needed to 

4 _ 1 _ 1 "+■ 

define the calibration equation B^ = ) S (V - V q ), previously 

derived, that allows the in flight determination of orthogonal field components 
from sensor output voltages. Re-expressing the instrument response equation 

-4 - J *4 4 

as B = S V = QV, we see that the data from the n measurements made can 
be rearranged in the following way: for each applied field axis k = 1, 2, 3 
the k-th field component can be expressed in terms of the three measured 
voltage components as indicated below. The unknown Q k are the elements of 
the k-th row of the inverse response matrix Q. For n measurements we have 

4 -4 4 4 

n such equations which in vector form are B k = V Q k where B k , Q k are n 
component vectors and V is a 3 x n matrix. 

This form is identical to that used in the determination of S, and by direct 

analogy the least squares error solution for Q is determined by 

-4 t - 1 T "► 

Q k = (V 1 V) v 1 B k . 

Statistical Weighting - The results above assume uniform statistical 
weights; that is, each measurement is given equal emphasis on the assumption 
that each of them is of equal importance and has equal probability of error. 
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Section 6 - Calibration and Alignment 

DETERMINATION OF THE RESPONSE MATRIX (Continued) 


This may not be the case; for example, perhaps during flight one axis of 

the magnetometer is more likely to be measuring low fields most of the 

time. Given such a situation, low field calibration data for that axis should 

be emphasized. This can be accomplished by minimizing the weighted least 

squares error (A WA) with appropriately chosen weights W. In this case 

-+■ 'p -IT 

the least squares error solution for S becomes S = (B WB) (B W) V. 

Estimate of Precision - The best estimate of the precision of the deter- 
mination of is w A^jy^n - assuming gaussian statistics and 

n measurements. The rms deviation of the experimental data from the 


These same for mu - 
he appropriate change of 


best estimate value (non-linearity) isJ^A^ 
las apply to the determination of the inverse with tl 
variables. 

The least squares error estimates of S and S * can thus be obtained as 
described above, complete with estimates of the variation of these quanti- 
ties over the applied field range. From them the scale factors and orienta- 
tion of sensor axes as well as the calibration equation can be derived as 
described, provided the associated axis distortions and rotations are 
known. We next consider the determination of these parameters. 
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TABLE 6- 3. SUMMARY OF EQUATIONS 

-*■ -*■ 

V = SB 

for each sensor axis k and the i-th measurement: 

ik B H S kl + B i2 S k2 + B i3 S k3 i = *» • • • » n 

- ^ k 

where S^ is the k-th row of response matrix S 
due to measurement error! 


Bs k - v k = 


minimize square error! 

: k -k 
solution is: 


/t T t \ T T -*■ 

K A kj = ° ~ 


s k iBiBs k- s k T bT V V B s k + V 


S. = (B T B)" 1 b t V 


and 


S = 


for inverse Q = S' 1 ; B = QV for each field axis k and the 
i-th measurement: cne 

B ik = V il a kl + V i2 Q k2 + V Q i . 1 


B k = VQ fc 


by analogy solution is: 


Qu = (V T V) -1 v T B. 
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Section 6 - Calibration and Alignment 


HOW ALIGNMENT IS DETERMINED 


Alignment stability can be determined from a matrix equation relating 
changes in sensor axis orientation to the response matrix S and an opti- 
cally measured rotation matrix. An unanalyzed technique for determining 
absolute alignment makes use of a precision gimbal to also relate the 
desired parameters to the response matrix S. 


The instrument response matrix previously dissected into components 
and measured as described relates the output of the sensor axes to the 
magnetic test field axes. To obtain the sensor alignment and its stability 
over temperature and shock, the reference to test field axes must be 
removed since only parameters internal to the sensor assembly are of 
interest. For both absolute and stability measurements this will require 
an optical system to measure changes in orientation of the sensor assem- 
bly with a precision of better than 1 arc-second. 

Rewriting the expression for the response matrix as 


S = E 


( T MV R BM °B 


) ” E * t mv t bm ) 


where T. 


-1 


MV * d v r mv and t bm - r bm d b'*- the " t mv iB 4 ° f the 

direction cosines relating sensor axes to mirror axes. The change in this 
matrix AT^V *- s what we wish to determine for stability, and its value 
describes the absolute alignment of the sensor axes. We next consider 
the properties of infinitesimal matrices that will allow determination of 
AT mV wit hout reference to an absolute external coordinate system. 


Infinitesimal Matrix Properties - An infinitesimal matrix is one whose 
second order differences from the identity matrix I can be neglected. Infin- 
itesimal distortion and rotation matrices commute, XY-« YX» (X + Y - I). 

If the initial alignment of the sensor assembly is kept within 0. 1 degree, 
then the off diagonal elements of Dy, Dg, Rgy are less than 2. 2 x 10 
and they can be considered infinitesimals. 

The equation for S may then be rewritten 


S “ E ( T MV + r bm + d b " 2I) 

where T^ = D y R My - 

Alignment Stability - Changes in the alignment due to heat or shock 
change both TmV = Dy + Ryv - I = Dy RmV an< ^ ^BM> The change in 
Tmv, ATmVi which consists of the change in direction cosines of the 
sensor axes relative to the mirror axes, must be determined. Writing 
the difference equation 


AS = S' - S 


e < at mv + ar bm 


21 ), 


we see that if S and s' are measured, only ARgy is needed to obtain the 
desired result. SinceAR BM is just the change in orientation of the mirrors 
mounted on the sensor assembly, this can easily be determined by two 
theodolites monitoring their orientation along two orthogonal lines of sight. 
This measurement thus does not require any absolute reference since only 
changes are needed. 

Absolute Alignment - The specification for the absolute orientation of 
at least two of the sensor axes requires a well defined external mechanical 
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of ef s"„ e s T^ the orientation 

axis relative to the external reference w' IT ?'T ledg f of the orientation of either 

optica, system that I t'he^r “ ““ 

relate the mirror axes tr, mo „i, . 1IJirrors also either defines or can 

entation with respect to the mirrors* 11 ^! if* er ® nce » t * le absolute sensor ori- 
terms, expressing Rbm d R -1 as Tn W1 ther^T ee ” If terr ^, ined - In matrix 
to the external reference fvstem aSMh^TL TBM ~ TxM Tfix where X refers 
Expressing S as Y a “ d the T 8 are ^trices of direction cosines. 


S = E (T T 1 


then 


MV 


= E _1 S T 


-1 


BM 


E" 1 S(T 


XM T BX ) 


-1 


and the absolute alignment is determined when TVw and T j a 

One possibility for a suitable x XM and ^BX are determined, 

rotational axes ofYprecis!^ X l ? that defined ^ the 

will allow theodolites to obtain the transfn " PPropriately mounted mirrors 
obtained^ follows: R^rrtog to ^ TbX m * y b ' 
around w by a precisely known b f low »_ rotating one sensor axis 

the applied Field B with resnect fn Y et ® rmine th ® direction cosine of 

tion around an orthogonal axis will therTd r & ° rn ? u ^ a S* ven ' Another rota- 
so on. In this way direc io^ cLf r ? n °, ther direction cosine and 
determined and Tg X is defined. nCS ° 6aC ° f the teSt field axes can be 

t h a 1 c n expected*,?” “ h ” ^oT bee' W °. rke , d a » d > h « precision 

procedure may be adequale^f a .JhdX' IcTl 'o",'?*' 1 '”’ 
gimbal mount can be obtained. Y P ec Se L ~ or 3_axis non-magnetic 

TABLE 6-4. SUMMARY OF EQUATIONS 


s = E(d v r bv d b‘ 1 ) = E ( D V R 


'V a MV R BM D B 


-1 


} E(T MV R BM D B -1 


) ~ E(t mv t bm ) 


Stability Equation: AT^ * (E^AS - AR fiM + 21) 


Alignment Equation: T..,. = E -1 SIT T \ _1 

MV ' 1 XM 1 BX ) 


Figure 6-1. Absolute Orientation. 
The direction cosine of an applied 
field B, cos x, can be determined 
relative to a rotation axis w by two 
sensor output measurements: 

Vj • B = \B cos a and Vjj • B = X.B 

cos b separated by a precisely 
known sensor rotation c. 
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Section 6 - Calibration and Alignment 


SUMMARY OF ACCURACY AND EQUIPMENT REQUIREMENTS 

To obtain the required accuracy of measurement the GSFC magnetic test 
facility must be upgraded to 0. IT field accuracy and the orientation of its 
magnetic axes must be specified to 0.4 arc-seconds. 

The accuracy desired in the determination of this magnetometer s 
scale factors and alignment exceeds that of currently available instrumen- 
tation. Almost every parameter must be measured to within a few parts 
per million (ppm). Using the best available equipment and careful statis- 
tical handling of the data obtained from it should allow the requirements to 
be met. It may be necessary, however, to relax the requirement for the 
absolute alignment of sensor axes to within 1 arc-second if the orientation 
of the test facility field axes relative to some optical axes cannot be deter- 
mined with better accuracy. The requirements due to each of the sensor 
parameters are enumerated in the table opposite and discussed below, 
assuming sensor scale factors on the order of 150 pV/Y. 

Offset - Since the magnetic field in the GSFC test volume can be set to 
zero only within 0. 2Y, measurement of the offset to higher precision must 
be accomplished by reversing the magnetometer axes. For a uniform 
field where slight translational differences between "up" and "down" mea- 
surements are negligible, the accuracy of reversal need not be great — 

10 degrees error in angle only produces a 1% change in the measured off- 
set for 0. 2 Y residual field. The factor limiting the precision of the offset 
determination is therefore likely to be the voltmeter used, which must 
measure ~60 pV . 

Noise - To measure sensor noise of 0. 3 Y over the dc to 25-Hz bandwidth 
in a 10-second period will require a strip chart recorder and voltmeter 
capable of measuring 45 pV at up to 25 Hz. Also required is a mu-metal 
flux tank to reduce ambient field variations by a factor of 300 or better. A 
superconducting shield would be even better and could also be used to deter- 
mine offset, but the large sensor diameter (6 inches) may make this 
impractical. In any case, a flux tank is adequate. 

Scale Factor and Linearity - The scale factor must be determined to 
20 ppm, which allows 10 ppm uncertainty in each of the magnetic field and 
voltage measurements. Available voltmeters have accuracies of ~20 ppm 
and current GSFC field accuracies are also ~20 ppm. Neglecting the effects 
of drifts during multiple measurements, the improvement in precision 
from n measurements varies as (n - m)-l/2 where m equals the number of 
parameters to be determined. Since the response matrix has nine elements, 
a factor of two improvement requires a minimum of 13 separate measure- 
ments at each field value. 

Alignment - The alignment stability of 5 arc -seconds and absolute 
alignment of 1 arc-second is the hardest specification to meet. With high 
quality theodolites and mirrors, the changes in orientation of the sensor 
assembly with respect to the optical axes can be monitored to 0. 1 arc- 
seconds, which is adequate. The orientation of the magnetic axes involves 
knowledge of both the magnetic field and the sensor output voltage. Since 
1 arc-second is 5 ppm, both quantities must be determined to better than 
5 ppm. 

A field of 50, 000 Y must therefore be accurate to better than ±0. 25Y, 
and corresponding voltages of 7. 5 V must be measured to better than 
±38 pV absolute. The required field accuracy may be obtainable by taking 
100 measurements with the proton reference magnetometer at GSFC. For 
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a field accuracy of ±0. IT and voltmeter accuracy of ±10 MV, the remaining 
error budget of ±0.4 arc-second from an RSS total of 1.0 arc-second repre- 
sents the total allowable error in determining the optical orientation of the 
magnetometer mirrors relative to the GSFC field axes. 

Accuracies on the order of a few arc-seconds in the determination of the 
field axis orientation with respect to an optical reference may be attainable 
by the previously suggested gimbal mount technique. This is not within the 
allowable 0.4 arc-second budget, however. Further analysis is required to 
determine the accuracies attainable by this or any other suggested procedure 
for the absolute alignment. 

TABLE 6-5. REQUIREMENTS 


• Offset 

0. 2 7 null field 

60 mV voltmeter 

180° ±10° gimbal mount 

• Noise 

Strip chart recorder 
45 mV voltmeter 
Mu-metal fluxtank 

• Scale Factor 

64, 000 ±0. 57 magnetic field 
7.5 V ±100 mV voltmeter 

• Alignment 

6. 4 x 10 ^ ±0. 17 magnetic field 

7. 5 V ±10 pV voltmeter 

0. 1 arc-second theodolites 

1 arc-second 2-axis gimbal mount 

0.4 arc-second GSFC field axis 
orientation 


Note * The analyses of this section have assumed a single range analog elec- 
tronics system and an error budget appropriate to that approach. For a 
multiple range OFG system the budgeted errors can be larger since part of 
the ADC function is included in the sensor electronics. The calibration volt- 
meter requirements are less stringent, but each sensor range must be sep- 
arately calibrated. 
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SIGNAL- TO-NOISE RELATIONSHIPS FOR RING-CORE SENSORS 

Despite contrary suggestions in the fluxgate literature, it is possible to make 
quantitative estimates of the signal-to-noise ratio for ring core ® enSOr h ‘ Droduct 
noise equivalent field strength varies inversely as the square root of the pro 
of frequency a nd core volume* _ — 

The fluxgate literature contains no clear-cut analysis of the signal-to- 
noise performance for fluxgate sensors. This lack is apparently caused by 
the fact that it is customary to cite an early work by William and Nob ( ) 
on Barkhausen noise which does not predict effects of the observed 
tude. In fact, recent explanations of Barkhausen noise do appear to Predict 
the observed fluxgate signal-to-noise ratio, especially when values of t e 
demagnetizing factor appropriate to the cylindrica geome ry are ^ * 

The most authoritative review of Barkhausen noise appears to be that 
of H Bittel, who reviews earlier work on the effect. The principal reason 
why the effect appears to be larger than at first thought is that the elemen- 
tary magnetic dipoles in the material are not reversed one by one but 
rather in clustered groups containing up to 10+1? elementary dipoles. The 
effect of this clustering is to produce fewer but much larger noise pu 
which causes a much higher noise signal than that calculated by G. W. 
McFarlane and reported by Williams and Noble. As a simple approxima- 
tion to a very complex phenomenon, Bittel calculates that the Barkhausen 
noise energy is enhanced over the white noise expected from the ^PP^S 
of single dipoles by a factor porportional to n c vdl 8 , where n c is the nu 
ber of domains interacting in a cluster, v, is the domain volume, and I s 
is the saturation magnetization (expressed in the same units as , 

that The work on which Bittel based his review was performed on gram- 
oriented 3% Si-Fe, a very noisy material unsuitable for fluxgates Cbut 
ideal for the study of Barkhausen noise). Fortunately, T. Higuchi (b) 
performed measurements on a number of core materials used m magnettc 
amplifiers and modulators, including Supermalloy and Mo-Permalloy as 
weU as grain-oriented 3% Si-Fe. The noise for Mo-Permalloy is about 
500 times less than it is for the 3% Si-Fe, which, comparing with Bittel s 
paper, leads to a value of n c v d for Mo-Permalloy of 4 x 10^2, less by a 
factor of 25 x 10^ than for the 3% Si-Fe materials. 

One may now calculate the signal-to-noise ratio expected in a fluxgate 
sensor. To keep things simple, we have assumed that the sensor preampli- 
fier is untuned so that all even harmonics are kept; in practice, this pro- 
cedure would probably give rise to high order harmonic sensitivity and 
phase instability, but the signal-to-noise ratio should not be greatly dif- 
ferent from that of a practical system. 

The result, derived in detail in the subsequent pages, is shown oppo- 
site. The noise-equivalent magnetic field varies directly with the satura- 
tion field of the core and inversely with the demagnetizing factor and the 
square root of the core volume, drive frequency, and integration time. 

This result checks quantitatively as well, predicting noise levels of -0. 04Y 
rms with typical core geometries. (A slight over estimate due to various 

simplifications. ) 
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THE NOISE EQUIVALENT FIELD 
FOR A FLUXGATE MAGNETOMETER 


The noUe-eqrnvalem magnetic field for a ring-core fluxgate is approximated 
by the following formula. The value of n c v d is about 4 x 10-12 cn ^-3 for the 
Permalloy used by Higuchi. 


fn v 
c d 


AB n * fV# B s 


= cluster number 
= domain volume 
V = core volume 
f = drive frequency 

T = integration time (l/(2x noise bandwidth) 
B g = saturation field of core 
K = demagnetizing factor 
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EVALUATION OF THE DEMAGNETIZING FACTOR FOR RING- CORES 


The demagnetizing factor for radially thin ring-core sensors can readily be 
calculated in cylindrical geometry from first principles. The correct result 
differs materially from the long cylindrical rod result often quoted in the 
literature. 


The field inside the core is not as high as it would be in an infinite 
parallel-plane geometry because of the influence of end effects. For long 
rods, the demagnetizing factor is given by Bozorth and Chapin (7). While 
their formula is usually used in the literature to apply to ring cores, it is 
better to use the correct geometry if one wants correct results. 

The ring core behaves as a magnetic shield, blocking the field from 
its interior by guiding it inside the core material (see sketch). Thus, it 
has what amounts to a capture cross section for the incident flux lines. 

For high enough permeability, all of the field is contained in the metal and 
further increases in permeability do not increase the captured flux. Since 
the width of the core tape is large compared to thickness, an infinite- 
cylinder calculation should provide a good estimate of the captured flux. 

Consider an infinite cylindrical magnetic shell with isotropic, homo- 
geneous relative permeability K m = \±/\± Q , placed in a magnetic field AB, 
normal to the axis of the cylinder. Let the inner and outer radius of the 
shell be a and b respectively. Then, the magnetic scalar potential outside 
the cylinder, inside the iron, and inside the cylinder must satisfy the fol- 
lowing equations. 

2 2 

Outside: ft = -( ABr + Ab /r ) cos 0 (1A) 

m 2 2 

In the iron: ft = -(Br + Cb /r ) cos 0 (IB) 

m 

Inside: ft = -Dr cos 0 (1C) 

m 

The constants A, B, C and D are as yet unknown. This form follows 
because, at infinity, the field must be all in the x direction, so that the 
scalar potential must have the form 

ft = -AB x = -AB r cos 0 (2) 

m 

2 

The terms in cos 0/r are the only other solutions of Laplace ! s equation 
allowed which have the same angular dependence. 

The four constants A, B, C, and D can be found by matching the mag- 
netics boundary conditions at r = a and r = b: namely, that the normal 
component of B = -Aft m as well as the tangential component of H = B/jjl 
must both be continuous. The exact result for the field inside the iron at 
right angles to the outside field is 

AB' = 2K m B o [(K m +l)b 2 /a 2 + (K m - l)b 2 /r 2 ]/[(K m +l) 2 b 2 /a 2 - (K m - l) 2 ] 

(3) 

If we write this in terms of the diameter for a radially thin core of high 
permeability by putting 

a = b = d/2, (b - a) = t, and K » 1, the result for the average 
field reduces to 

AB r * K AB/ ( 1 + K t/d) = AB(d/t), K » 1 (4) 

m m m 
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( This result, while still approximate, is more nearly correct than the usual 

long rod approximation in which AB' varies nearly as d^/wt, where w is the 
core width. The "demagnetizing factor" is thus K = t/d. The capture cross 
section is simply the geometric cross section wd. 

We shall now apply this result to calculate the sensitivity of a rine-core 
fluxgate. & 



Cylindrical Geometry used in Calculation. The 
ring core is approximated by a thin-walled cylin 
drical shell. 
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MINIMUM SENSITIVITY OF THE RING CORE SENSOR 


Simple formulas for the minimum sensitivity can now be derived by taking 
advantage of the square -loop approximation and the knowledge of the demag- 
netizing factor. 


To calculate the signal-to-noise ratio, we will need to know the sensi- 
tivity of the fluxgate sensor. 

Estimation of Sensitivity - The sensitivity to ambient field variation 
has been calculated by a number of authors, not always with consistent 
results (1,2, 3,4). We prefer a simplified viewpoint. Instead of calculat- 
ing the second harmonic output, we will calculate the output of the even- 
harmonic synchronous demodulator assuming use of a broadband amplifier 
so that all even harmonics are used. This calculation is easily performed 
in the time domain without the need for Fourier representation (which 
complicates the noise analysis). To simplify matters, we assume use of 
a square-loop core, and a triangular current drive waveform which results 
in a rate of change of flux and concomitant induced voltage along the verti- 
cal sides of the loop, and zero induced voltage elsewhere. The derivation 
is as follows (see the accompanying figure). 

First consider the output of a single half of the core with no field out- 
side. As the field increases along the vertical side of the square loop from 
-B s to B s , the pickup coil produces a voltage output pulse of amplitude 
E = NAB and width 2B S /B,^ where N is the number of turns on the pickup 
coil, A is core area, and B is the rate of change of magnetic field. A 
corresponding negative pulse is produced as the field sweeps through the 
opposite side of the loop. This zero-field output signal contains only odd 
harmonics of the drive field. 

When the core is biased by an external field, AB so that the internal 
field is biased by AB 1 , the pulses tend to separate in one half of the core 
and come together # in the other half, as shown. The pulses move by a time 
interval At = AB’/B. When the pulses from the two half cores are differ- 
enced by appropriate connection of the pickup coil (simply by winding it 
around the core), the difference signal becomes a doubled-frequency sig- 
nal, as shown, with no odd harmonics of the drive frequency remaining. 
When demodulated by an appropriately phased synchronous demodulator, 
and averaged for a time T, the contribution to the average for each cycle 
is 8EAt/T. For a total of fT cycles, we obtain the output voltage 

e = 8fNA AB f = 8fNAd AB/t = 8fNdw AB (5) 

s 

where we have used the demagnetizing factor derived in the last section to 
express AB' in terms of the external ambient field AB. The symbols d, w, 
and t represent core diameter, height, and radial thickness respectively. 

The derived sensitivity is a minimum because we have neglected the 
contribution of partially saturated regions which still have enough per- 
meability to contribute signals but which do not contribute Barkhausen 
noise. 
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Analysis of Fluxgate Sensitivity. Fluxgate Waveforms are shown 
for an idealized square B-H loop and untuned amplifier. 
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BARKHAUSEN NOISE AND THE NOISE-EQUIVALENT FIELD FOR RING- 
CORE FLUXGATES 


Using a simple model for Barkhausen noise, calibrated for typical fluxgate core 
materials by reference to published data, we can calculate the Barkhausen noise 
and noise equivalent field for the ring- core sensor. 


We can now see how the details of the noise calculation are carried out. 

Noise Calculation - Because the simplified model for Barkhausen noise 
is really only a shot noise model (with cannonballs for shot) we can use a 
time domain method of evaluation. 

The effect of the demodulator on the noise is insignificant, since all it 
does is reverse the sign of pulses which are already bipolar. Thus, all 
we need do is calculate the effect of integrating the noise energy. Bittel (5) 
has already done this, in effect; his formula for the noise energy from the 
pickup coil, integrated along one side of the hysteresis loop, is 

E = 8n v,N 2 (A/£)B 2 (V 2 - s 2 ) , (6) 

n c d s 


where n c is the average number of domains per noise cluster, vj is the 
average domain volume, N is the number of turns of the pickup coil, A 
and l = TTd are the coil cross-sectional area and length, and where we have 
replaced the saturation magnetization by the saturation field B s , because 
the permeability is high. 

The noise voltage per cycle generated at the output of an integrator 
can be calculated from CampbelUs theorem: 


t ■ / 


A 2 h 2 (t)Rdt 


(7) 


where R is the rate of arrival of events of mean- square amplitude A, and 
h(t) = 1/T for 0 > t > T and is zero otherwise. BitteUs formula is already 
integrated, so that to apply it we need only to multiply equation (6) by (1) 
a factor fT to allow for the total number of noise bursts; (2) a factor 2 to 
allow for the differential pickup winding; and (3) l/T^ to normalize h^(t). 
The end result for the noise at the output of the integrating filter is 


; 2 = 16 4 n v, N 2 (A/£)B 2 
n Ted s 


( 8 ) 


Noise Equivalent Field Strength - By setting the signal-to-noise ratio 
equal to unity and solving for the ambient field strength, we can find the 
noise equivalent field strength. Formally, 

AB = e /(ere / <rAB) (9) 

n n s 

The result of performing this operation using the sensitivity and demag- 
netization factor K = t/d derived earlier, is 



Here, as before, v is core volume; the factor t/d which represents the 
demagnetizing effect is kept separate to show the fundamental dependencies 
on core volume and demagnetization separately. However, the volume v 
for a ring-core is TTdtw; when this relation is plugged in, the result is 
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AB 

n 


1 / n c V d t 
2\l wfT d 3 


B 

s 


( 11 ) 


It is good to make the core big, thin, and wide. This is a new result, since 
the thin-rod analysis presented in the literature leads to the conclusion that 
the core width w should be small. 

In the table below we present the calculation of the actual noise- 
equivalent magnetic field strength for a 1971 Time-Zero sensor which exhib- 
ited a noise for a 1 Hz single pole bandwidth of 0. 14V peak to peak. The 
calculated rms noise equivalent AB (NE AB) seems only a factor of two high; 
in view of neglected signal contributions and metallurgical differences between 
Time-Zero in California, Bittel in Germany and Higuchi in Japan, the agree- 
ment seems rather good. None of this rationale has been published but it 
will be. 

Conclusion - Fluxgate noise levels are consistent with observed values 
for Barkhausen noise. The noise varies as (t/d^wJl/Z. 


Calculation of NE AB for TZ Sensor 


d 

t 


w 

f 

T 

B 

s 


AB 

n 


4 x 10 ^ cm ^ 

1. 5 inches = 3.8 cm 
0.001 inches - 0. 0025 cm 
0. 28 inches = 0. 71 cm 
15 kHz 


0,5 (1 Hz noise 
0. 8 Tesla 

1 /Vd r~ 

2 V " fT d 3 w 


BW) 


B = 
s 


0. 04 Y rms* 


This value is for ordinary 4-79 Permalloy, used by Higuchi and in the TZ 
sensor* The value of is about two orders of magnitude lower for 

well-annealed 6-81.3 Permalloy. 


A-9 


ESTIMATE OF MAGNETOSTRICTIVE OFFSET 


A simplified estimate of the magnetostricti ve offset shows that the offset, like 
the noise, is reduced by use of thin tape an d a large diameter bobbin. 

It is generally agreed that the principal source of residual offset in a 
properly constructed fluxgate sensor is magnetostricti ve. The research 
conducted at the Naval Ordnance Laboratory in developing the lunar sur- 
face magnetometer for Apollo showed that the best offset performance was 
attained by using as a core material a permalloy containing 6 percent 
molybdenum and 81.3 percent nickel. For this nickel concentration, the 
longitudinal magnetostriction is zero^, Even when such material is used, 
some residual effect may be present, so it is important to understand how 
the offset is generated, Weiner® has presented a theory of the effect, 
assuming sinusoidal waveforms and a particular shape for the magneto- 
striction curve. Weiner also presents data showing that the effect varies 
with drive level, peaking if the core is driven just to saturation, and 
becoming much smaller if the core is overdriven. 

For permalloys whose nickel content is near 80%, the magnetostric- 
tive elongation is small until the material saturates, when the elongation 
increases to the (small) saturation value \ s . It seems reasonable to treat 
the magnetostrictive curve as a step onset at saturation (Figure A). Then, 
the flux and magnetostriction during the cycle will vary somewhat as shown 
in Figure B. 

The voltage induced in the second-harmonic pickup coil is given by: 

e = NA ^ + NB M. ID 

where N is the number of turns on the pickup coil, A is the core cross- 
sectional area, and B is the field intensity. The second term is the 
derivative of the magnetostrictive flux change term in Figure B. The 
magnetostrictive part of the voltage waveform (Figure C) consists of a 
series of spikes at the second harmonic of the drive frequency. 

The area of each spike is roughly given by 



NB AA = NB \ A 
s s 


When divided by the sensitivity calculated earlier, 
magnetic field units is 


B 

o 


I x B 7T 

4 s s d 


( 2 ) 


the offset expressed in 


where d is the core diameter and t is the tape thickness. This formula 
predicts that the offset is independent of tape height, a rough result because 
it is based on the infinite - cylinder calculation of the demagnetizing factor. 
For d = 1 inch, t = 0. 0005 inch, k s = 5 x 10" 6 (79% NiFe) and B s = 0. 8T, 
the offset estimated from this formula is 0. 5Y - not far from what is 
observed. The effect is temperature-dependent which presumably gives 
rise to the observed offset instability of a few tenths of a gamma. 

When the longitudinal magnetostriction is reduced by use of the proper 
alloy, residual volume magnetostriction still exists. This effect is very 
small, about 1 part in 1C)8 for the typical sensor drive conditions, and 
may be neglected. 

Thus, to obtain a sensor with stable offset characteristics, it is 
advisable to use as a core material a Ni-Fe-Mo alloy with 81. 3% nickel. 
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Figure A. Magnetostriction Model 
The length change is assumed to 
occur suddenly at saturation. B s is 
field intensity at the knee of the stiff- 
ness curve. 


Figure B. Variation of Flux and 

Magnetostriction During 
the Sensor Drive Cycle 
The flux variation curve assumes the 
use of a reactive pumping drive circuit. 




Figure C. Magnetostricti ve Error 
Voltage 

The voltage is proportional to the 
derivative of the magnetost rictive 
curve; the spikes occur at twice the 
drive frequency. 
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DETAILS OF STRUCTURAL ANALYSIS 





TIME •ZERO 


APPENDIX C 

ANALYSIS OF STRUCTURAL STABILITY 
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SUPPORTING STRESS AND VIBRATION ANALYSIS 
Sensor Structural Stability 


1. Deflection of spherical coil form 



Note that we are concerned with 
y deflection, not Slope G 


The sphere carl be conservatively modeled as a hollow cylinder with 
0. 2-inch thick walls. 

3 

W| 

then max y = -ygj- 

where E = 12 x 10^ psi for 
g r ap hit e - ep oxy 

1 - ! < R 2 4 - R l“> 



= l[(¥) 4 -(¥) 4 ] 


1 = 0. 17664 


(273 turns No. 36 wire) 
W = wt. of coils and coil form 


= (0.005) 2 (2.5tt) 6 x 273 x 0. 323 +-|ir(1. 3 3 - 0. 9 3 ) 0. 058 

= 0. 46 less 0. 03 for groove and hole 
W = 0.43 lb. 

*2 


y = 


(0. 43)(2. 3)' 

.6 


3 x 12 x 10 x 0. 17664 

,-6 


= 0. 83 x 10 ^ in. 


where 


Ay _ 0.83 x 10" _ ^ 0/ lrt -6 „ ,^-7 . 

"“yr ~ ^ — 3 - 0. 36 x 10 ~4 x 10 in. /in. 

3, 6 x 10 ^ A ni i 

p— »0. 01 arc-second 

5 x 10’ 5 


_ ^ 

5 x 10 in. /in. = 1 arc-second 
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2. Deflection of coil with respect to mirror mount. 

P 2 Pi = 0.449 


Pj = Wt. of feedback coil + Wt. of 
tower stand 



“1 


= 0.43 + -g (0. 6) x 1 x 0.058 = 0. 449 lb. 


F.B. COIL J P 2 _ 2P 1 = °* 898 lb * 


y = 


Wi 

3 El 


where E = 20 x 10 graphite - epoxy with 
fibers oriented axially and 
random 


I = f (0.7) 4 = 0. 189 


y = 


(0. 898)(1. 2) 

3 (20 x 10 6 )(0. 189) 
>-7 


Ay _ 1.4x10 
y 1.2 


= 1.4 x 10" 7 in. 


-7 . 


= 1 x 10 in. /in. 


Dynamic 


Random load factory, = 2. 2 (12. 9 rms gs) = 28 
Sine = 20 £ max. at 36 - 150 Hz 
Use 28 £ 


Find 

deflection of base 
77/////////////S. 

1 


V 

i 

'/////////////// 

h 4 h 


Flat plate clamped all edges dimen- 
sions 4 x 4 x 1/2-inch thick 


From Roark and Young, 5th Edition, Table 26, Case 8a and 8b 


8a uniform load over entire plate. 

.4 (0.0138)(ijP)(4)' 

max y = -a-y = =- 

Et (20 x 10°)(0. 5) 

max y = 0. 15 x 10 in. 


For sensor weight see below. 
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SUPPORTING STRESS AND VIBRATION ANALYSIS (Continued) 


8b load concentrated over small area of radius R q ~0. 2 u 


max y 


= oW bj = (0.06ll)(1.75)(4) 2 = 10>9x 10 -6 in> 
Et 3 20 x 10 6 (0. 5) 3 

.-6 


Actual displacement can be approximated as 1 x 10 because it will be 
much closer to 8a. 


Weight 

Mirror 


0. 079 lb/in 3 


Graphite - Epoxy 0.058 lb/in 
Mirror 

2 x 2 x 1/2 x 0. 079 = 0. 158 


Boom Cable 


20 685 

22 454 


= 1. 37 + 0. 10 (1. 37) = 1. 51 


Sensor 

Windings - No. 36 wire 273 turns 

Feedback Coils 6 t-| (0. 00 5) 2 x 2. 5^ x 273 x 0. 323 ] = 0. 104 

Cancellation Coils 6 [ -r (0. 005) 5^ (64)(0. 363)] = 0.038 

4 Total 0. 142 


Base 


4x4x0.5x0.058 = 0.464 


Tower 


| (0. 6) 2 (1) (0. 058) = 0. 016 


Sensor Support 


x tA = °' 019 

.. . „ . 0.058 202 

Cancellation Sphere - q x 


= 0.410 


Feedback Sphere 
Sensors = 0.015 


0. 058 176 _ , c , 

07063 x 454 " °* 357 

Miscellaneous = 0. 2 
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Total Sensor w/mirror 1. 75 lb. 
Resonant Frequency of Base 
Using curves by Maarten Vet 
1. 75 


P = 


4 x 4 x 0. 5 


= 0. 219 


. 1/2 


f » = ° (?) fcTTTi) 

/ 0. 5 \ / 20 X 10 6 \ 

V 4.2 / \(0. 2 1 9) (0. 99) / 


= 35 


u ~ 0. 13 

1/2 


2 1 9 ) ( 0 . 

= 10, 500 Hz 

Assuming cantilever spring and mass — very conservative 

Using Steinberg comparison method 

, 1/2 


f = f, 

n 1 


0. 0768h 
ab 


(—^j) 

V P (1 - uV 


= 45l (0. 0768)(0. 5) /_ 
4x4 \0. 


20 x 10 


1/2 


219 (0.99) 


= 10, 400 Hz 

Safe to assume f^>4000 Hz, which is twice highest input frequency. 


n 

Deflection and f of Sensor Support 


§- 

I 


a 


1 


W * 0.034 


•1.25 


W = 0.019 + 0.015 = 0. 034 1b. 

See Roark and Young Table 3 Case la 

E = 42 x 10^ Fiber axially oriented 
graphite -epoxy 


max y = 


_ (0. 034)(1, 25P 

3EI 3 (42 x 10 6 ) \ x (0.2) 4 


= 0.42 x 10~ 6 in. 
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SUPPORTING STRESS AND VIBRATION ANALYSIS (Continued) 


-&L = £4| x 10" 6 = 0. 34 x 10" 6 
y 1. 25 

, 1 m 1/2 1 / 386 ' 1/Z 

n 2ir\fj 2 tt 'o.42x 10' 6 / 


f = 4800 Hz 
n 


Maximum Dynamic Stress in Support 

2! 


<r = 28 (0* Q34)( 1. 25)(0. 2) = lg9 pgi 
\ ( 0 . 2 ) 


Well below allowable of 50, 000 psi. 

Dynamic deflection 

max y = 28 (0. 34 x 10 ^) = 9. 52 x 10 ^ 
Maximum Sensor Stress 
Occurs in baseplate 

£ f = 28 K 


<7 = 


B 2 W _ (0. 7542) (28 x 1. 75) 

2 2 
tr (o. 5j 


Load applied over small area 
of center of plate. Maximum 
stress occurs at center of 
edge of plate. See Roark and 
Young, Table 26, Case 8b 


Q = 148 psi 


This is lower than stress in the sensor support. Maximum stress occurs 
in sensor support and is a = 189 psi. 

COIL WIRE STRESS 

y g = 10, 000 psi copper coil wire 
€ _ s 10, 000 = 5. 88 X 10' 4 in. /in. 

" E " 17 x 10 6 
D = 2. 5 in 

AD = 2. 5 in (5. 88 x 10 -4 ) = 1. 47 x 10' 3 in 
TC = 17. 64 x 10- 6 =oc 
AD = aDAT 


AT = 


AD 
a D 


1.47 x 10' 


(17. 64 x 10 ) (2. 5) 


= 33. 3 F 


AT = | (33. 3) = 18. 5°C 


maximum dynamic range coil wires 
can withstand 
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Required dynamic range -50°C to +85°C = [135°C] 

If coils are wound at 27°C with no pre-stress and we assume the graphite- 
epoxy form is significantly stiffer than copper, then 

- 50°C - 27°C = [77°C] 

and copper would yield. 

Copper coils must be wound at: 

-50 + 18. 5 = -31. 5C 

\ 

to avoid yielding when temperature is reduced to -50°C. 

Since coil wire is originally at 2. 5 in. dia. it will try to increase to a 
larger dia. 

AT = -31. 5 - 85 = [116. 5] 

AD = aDAT = (17. 64 x 10' 6 ) (2.5) (116. 5 x |) = 0. 00925 in. 

D 

However, it is not practical to wind coils at -31. 5°C. Wind at 21°C. 

-50 + 21 = (7 1°C) 

85-21 = (64°C) 

AD = (17. 64 x 10~ 6 ) (2. 5) (71 x 9/5) = 0. 0056 in. 


Y t = 32 - 35 ksi 
S = Ee 

S = 17 x 10" 6 (0.0056/2.5) = 38, 080 ksi 
This could cause failure of some coil wires. 


So back coil form groove with elastomer to allow coil 
without inducing stress above yield. 


wires to shrink 


At high temperatures, hold coils in place with flexible adhesives. 
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APPENDIX D 

PACKAGING OF THE ELECTRONICS 



D-l 



ANALYSIS: PC BOARD RESONANT FREQUENCY AND STRESS 
Circuit Card Weight . 

454771b = °- 176 lb 

Resonant Frequency 

Using Steinberg comparison method 


4 in x 4 in x 5 g/in^ x * 



f = f 1 
n 1 


9.0768h 

ab 


(— s-) 

'P(l - U^) f 


1/2 


where 


h = 0.07 a = 4 b = 4 
E = 2 x 10^ lb/in 2 for epoxy glass 
P = 0. 12 
n 0. 176 


4 x 4 x 0. 07 


= 0. 157 lb/in 3 




From Table in Steinberg 


/ / fs . 1 / 2 

f = (0.0768) (0.07) / 2 x 10° \ 

n 4x4 \0. 157(0.9856)/ 


= 211 Hz 


for f of 2 1 1 
n 

Q = tf n ) 1/2 = (2 1 1) 1/2 ~ 14 

PC Board Stress 
Load Factor 

Random ^ = 2.2 (12.9) = 28 
Sine = Q ^ in = 14(5) = 70 

Dynamic Load Intensity 

W£ 


_i . (0.157M70) . 0 . 687 1b/l „2 

a b 4x4 


Bending Moment 

q 


M = 

y 


Dynamic 'Bending Stress 


7 ' - 0.312 in. 


(4 + 4) 

2/l 1 \ 

Va 2 b 2 ! 

MI6- + it) 


lb/ in. 


s = = U50pai 

b h 2 (0.07) 2 

~ 1200 psi 


K,p = Stress concen- 
tration factor 

h = board thickness 
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Overall Mechanical Design 


The general design configuration as presented in the original Time-Zero 
Laborat6ries proposal is sound and will remain basically unchanged. This 
mounting structure consists of a baseplate, sensor mounting tower, a 
feedback coils support sphere, and a field cancellation coils support sphere. 

The following sections describe the rationale for materials selection, 
the influence of this materials selection on the detail design, and an 
associated performance analysis. 




Materials Selection 


The combined requirements of low structural weight and ultra-low coefficient 
of thermal expansion make the selection of a fiber-reinforced, polymer-matrix composite 
a logical choice^ A graphite fiber is recommended because of its low density, high 
strength and stiffness, and in particular because of its highly anisotropic thermal 
expansion properties. The anisotropy of the graphite fiber is caused by prefered 
orientation of its graphitic (lamellar, plate-like) crystalline structure during 
manufacture. For example, mechanical stretching of the precursor fiber during 
the graphitization process will cause this orientation. One result is a high fiber 
axial stiffness and tensile strength relative to the corresponding transverse 
properti6s a Also, the axial coefficient of thermal expansion is very low, 
actually slightly negative in most cases, while the transverse coefficient of 
thermal expansion is always positive and relatively large, comparable in magnitude 
to the coefficient of thermal expansion of common Isotropic metals such as steel, 
aluminum, or magnesium. A wide variety of graphite fibers are commercially available, 
permitting the selection of one having the optimum combination of stiffness, 
strength, and thermal expansion anisotropies for a specific application. 

A number of different polymer matrix materials are also readily available, 
including epoxy, polyester, phenolic, and polyimide, to name the most appropriate 
candidates for the present application. These polymers are all isotropic, and 
differ little in density, strength, or stiffness from one to another. Polyester 
is lowest in cost, but is the least thermally stable. Polyimide is a high temperature 
polymer (300°C or higher service temperature), but is somewhat more difficult to 
fabricate structures from since it is not as of yet as well-developed as the others. 

A comparison of a variety of factors such as these leads to the firm conclusion 
that epoxy is the best choice for the present application. While a large number 
of epoxy resins are commercially marketed, their mechanical and physical properties 




do not differ significantly. Thus, unlike graphite fibers, the selection of a 
specific epoxy matrix material is not a critical factor. 

Graphite/epoxy composites of three different basic configurations were 
considered for the present application: unidirectional, angle-ply, and chopped 

fiber. Each configuration has distinct advantages and disadvantages. Which to use 
depends upon the component design configuration, as discussed in the following 
paragraphs. 

Baseplate - This component is a relatively simple, flat, rectangular plate, 

8.5" long, 5.0" wide, and 0.5" thick. Two alignment mirrors are attached to one 
face, and the sensor mounting tower to the opposite face. The principal requirements 
are a high out-of-plane stiffness (i.e., bending and torsional rigidity) and low 
weight. A uniformly low coefficient of thermal expansion (i.e., low in all in- 
plane directions) is an important secondary property, to minimize out-of-plane 
distortions due to thermal gradients across the surface of the plate. 

solid plate: 

volume = 8.5 x 5.0 x 0.5 * 21.25 in^ 

weight = 21.25 in 3 x 0.058 lb/in 3 - 1.233 lb 
The desire for thermal isotropy and high torsional rigidity suggests a [0, +45, 90] 
laminate, i.e., equal numbers of unidirectionally reinforced plies or laminae 
oriented at 0, +45, -45, and 90 degrees relative to the length axis of the plate, 
the total number of plies being as required to produce the desired plate thickness, 
and the resulting stacking sequence being symmetric about the thickness midplane of 
the plate. 

As a general guideline, the higher the (axial) modulus of the graphite fiber 
selected, the higher the in-plane (quasi-isotropic) stiffness of the composite, 
the lower the coefficient of thermal expansion (always slightly positive), the 
lower the in-plane strain to failure, and the higher the material cost. Commonly 



) 



available commercial graphite fibers range in modulus from about 25 x 10 6 to 
6 2 

80 x 10 lb/ in . Properties of typical graphite/epoxy materials, in isotropic 
(or more correctly, quasi-isotropic) laminated plate form, are given in Table 1. 

One very important property of these composites is their very low coefficient 
of thermal expansion, a, and their low density. For example, invar, a metal 
generally considered when low thermal expansion is required, has a coefficient of 

“6 o 

thermal expansion, ct, of 1.4 x 10 /°C, and density of 0.289 lb/in . That is, 

it has an expansion coefficient comparable to or even higher than that of the 
composites of Table 1, and is about five times heavier. 

One apparent anomoly in Table 1 is the lower value of a for T300/5208 relative 
to HMS/3501-5, this value being much lower than that of AS/3501-5 also. While the 
coefficient of thermal expansion does vary from one epoxy resin to another, this 
alone probably does not account for the total difference in Table 1. A slight 
variation in fiber volume content from one system to another is also a contributing 
factor. A +5 percent variation in fiber volume content is not untypical; slight 
differences in resin viscosity can influence the resin flow during the temperature- 
pressure-time cure cycle. As important, however, is the ability to accurately 
measure these very low values of a. It has only been during the past year or so 
that effort has been concentrated on developing reliable techniques for measuring 
very low values. Thus, much of the data in the literature, from which the properties 
of Table 1 were taken, are of limited experimental accuracy. There is no question, 
however, that the value of a for any isotropic lay-up, graphite/epoxy composite 
is very low. In any current program, this property should be experimentally 
measured for the particular material system being used. 

While considering thermal expansion coefficients, it should also be noted 
that fuzed quartz, the mirror material used in the present system, has a coefficient 
of thermal expansion of about 0.5 x 10~ 6 /°C. Thus, the isotropic laminate, 
graphite/epoxy baseplate matches this property of the mirrors very well. This 
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Properties of Typical Isotropic Lay-up [0,+45,90] g Graphite/Epoxy Composite Materials 
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will permit adhesive bonding of the mirrors to the baseplate, without developing 
subsequent thermal expansion mismatch stresses and distortions in either component. 

For the present baseplate design, any one of the materials of Table 1 would 
be satisfactory. All factors considered, the HMS/3501-5 system has been selected, 
based primarily on a good balance of elastic and shear moduli, an acceptable 
coefficient of thermal expansion, and a moderate cost. The excessive cost and 
therefore the lack of extensive use of the T75S/3501-5 system in prior applications 
has resulted in a limited data base, making the selection of this obviously 
attractive material unadvisable. Another high modulus graphite fiber, GY-70, 
produced by Celanese, Incorporated [4], is in a similar position of having a limited 
data base, although for different reasons. A typical composite form is GY-70/X-505, 
a graphite/epoxy prepreg produced by Fiberite Corporation [5]. While the cost is 
comparable to other, more commonly used systems, the high fiber modulus combined 
with a low to moderate strength results in a low strain to failure. That is, the 
material has developed a reputation for being somewhat brittle in certain 
applications. Efforts by the material supplier to increase the composite strength 
properties, by improving the fiber surface treatment and by modifying the epoxy 
matrix properties, are making this a steadily improving material system, however. 

Should structural weight savings become an even more critical consideration 
than at present, the baseplate configuration should be the first for reevaluation 
during detail design. Two changes could be considered, each of which would result 
in a significant percentage weight savings. The most direct change would be to 
switch to a sandwich panel configuration, incorporating a fiber/epoxy honeycomb 
core and graphite/epoxy face sheets. Shear webs radiating from the region of the 
sensor mounting tower attach point would be added to transmit out-of-plane bending 
loads and to minimize deflections. A sketch of this configuration is shown in 
Figure 1. Fiber/epoxy honeycomb core weighs in the range of 5 to 15 lb/ft 3 



3 

(0.003 to 0.009 lb/in ), depending upon cell size and fiber type. Thus, a total 
weight savings of as much as one-half pound relative to the solid plate design 
could be achieved, at a modest increase in cost. 

The second modification, to achieve an even greater weight savings, would be 
to make the alignment mirrors an integral part of the baseplate. The two quartz 
mirrors are approximately 4" x 4" x V' and 3V' x 3%" x 2" in size, respectively, 

3 

containing a total of about 25 in of fused quartz (after allowing for the non- 
planar contour of the latter mirror). Fuzed quartz weighs 0.080 lb/in 3 . Thus, 
the two mirrors have a combined weight of about 2 lb. Quartz is selected as a 
mirror material because of its good stiffness (E - 10.4 x 10 6 lb/in 2 ) and low 
thermal expansion (a - 0.5 x 10~ 6 /°C). As Indicated in Table 1, the various 
graphite/epoxy isotropic laminate composites have comparable or even superior 
in-plane properties. The mirror side of the baseplate could be molded to the 
required mirror contour shape, polished to the flatness desired, and metal-plated. 
This general technique is currently being evaluated for weapon systems applications 
by industry [6], By using this technique, most of the 2 lb. quartz mirror weight 
could be saved. 

Sensor Mounting Tower - This component (see Fig. 2) serves two principal functions; 
it supports the sensors and provides a mounting base for the two coil spheres. Its 
required properties are compatible thermal expansion coefficients, and a high 
flexural stiffness. 

The 2.5 diameter base ring which supports the field cancellation coil sphere 
(outer sphere) will be adhesively bonded (or possibly co-cured) onto the baseplate. 
Thus, it must have thermal expansion properties which at least approximately 
match those of the baseplate, and it must also be thermally isotropic so that it 
does not induce thermal stresses and distortions into the field cancellation coil 
sphere. Therefore, the 2.5" diameter base ring will be made of the same material 
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as the baseplate, i.e., an HMS/3501-5 isotropic laminate. The thickness of this 
base ring will be such as to match the effective wall thickness of the outer 
sphere it mates with, i.e., slightly greater than 0.10". 

The 1.25" diameter feedback coil base ring will be of the same material, for 
the same reasons. The depth of the sphere mounting step will be slightly greater 
than 0.20", to match the effective wall thickness of the inner sphere. The diameter 
of this step is 1.0". 

The cancellation coil base ring and the feedback coil base ring are to be con- 
structed of an isotropic laminate in order to match the radial thermal expansion of 
the cutouts in the spheres at their junction. The plies of the laminate thus lie 
in a plane perpendicular to the axis of the mounting tower. For the HMS/3501-5 

composite, the elastic modulus and coefficient of thermal expansion in this (tower 

6 2 —6 

axis) direction are 1.1 x 10 lb/in and 25 x 10 /°C, respectively. That is, 

this material has a low stiffness (low resistance to lateral flexural deflection 
of the tower axis) and high thermal expansion (poor thermal stability) in the 
axial direction. 

To offset these undesireable properties, the sensor mounting shaft will be 
made of undirectional GY-70/X-505 high modulus graphite /epoxy (or T75S/3501-5 
graphite/epoxy as a comparable alternate material), with the fibers oriented in 
the axial direction and the shaft running the entire length of the tower, as 
shown in Fig. 2. The properties of this material are given in Table 2. 

The combination of a slightly negative axial coefficient of thermal expansion 
and a very high axial modulus combine to offset the poor thermal stability of the 
mating base rings. The high axial (flexural) elastic modulus provides the required 
flexural rigidity of the tower. 

As given in Table 2, the transverse coefficient of thermal expansion of the 

unidirectional composite is very high, 31.7 x 10~ 6 /°C, and the transverse elastic 

6 2 

modulus is very low, 0.9 x 10 lb/in . 


This high value of a is a disadvantage 
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Table 2 

Properties of Unidirectional GY-70/X-505 
Graphite/Epoxy Composite [5,7] 

Elastic Modulus (10 6 lb/in 2 ) 


Axial 

44 

Transverse 

0.9 

Flexural 

41 

Shear Modulus 0,6 x 10^ 

lb/ in 2 

Coefficient of Thermal Expansion 

Axial 

-1.04 

Transverse 

31.7 


0.061 lb /in 3 


Density 




in the sense that radial interface stresses will be developed between the shaft 
and rings due to any temperature change. However, this will be alleviated signifi- 
cantly because of the low transverse elastic modulus of the shaft. On the other 
hand, the high value of a is an advantage since it will permit shrink- fit ting the 
shaft into the hole in the base rings. That is, the shaft can be cooled to shrink 
its diameter and then slipped into the hole. (There is no significant benefit 
in also heating the base rings since the coefficient of thermal expansion of this 
material is so low.) By controlling the shrink fit interference, a secure connection 
can be assured under all service temperature excursions above and below ambient. 

Feedback Coils Support Sphere - The principal structural requirement of the feedback 
coils support sphere is a uniform and very low coefficient of thermal expansion. 

Each of the three pairs of coil windings must remain parallel within a small 
tolerance, and remain at a prescribed spacing within a small tolerance. 

The requirement of uniform and low thermal expansion would be satisfied 
by any of the isotropic laminate composites listed in Table 1. However, it is 
not practical to fabricate a sphere in isotropic laminate form. It is difficult 
to maintain the fixed fiber-orientation angles around the surface. The most likely 
fabrication technique would be filament-winding; but this would not produce the 
required isotropy. The fiber angles at the winding axis poles would not be the same 
as at the equator. Hand lay-up of unidirectional prepreg material is also impractical 
because of the close tolerances required. 

A viable alternative is to mold the sphere from a chopped fiber/epoxy 
molding compound or structural mat. The latter is a prepreg material available in 
sheet or roll form. It can be cut into random-size pieces and layed over a spherical 
mandrel, then vacuum bag cured to the finish dimensions using spherical segment 
caul plates to control the external surface smoothness and wall thickness. While 
this is not a standard fabrication procedure in the composites industry, it is not 
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expected to present any difficulties. The mat material is prefered over the molding 
compound for the present application because of the better control of fiber orientation 
in the surface plane of the sphere. Less fiber flow will occur during the curing 
process. 

The chopped fiber in the mat prepreg material is randomly oriented in the plane 
of the mat. Thus, the composite exhibits in-plane isotropy. To achieve a very low 
coefficient of thermal expansion, it is necessary to use a very high modulus 
graphite fiber. This also maximizes the in-plane elastic modulus. One such material 
is hy-MAT 7534 structural prepreg [5] produced by Fiberite Corporation. The properties 
of this material are given in Table 3. The fiber is GY-70 graphite, chopped to a 
nominal length of 1 inch. The matrix is Fiberite 934 epoxy. The value of the 
isotropic coefficient of thermal expansion of 1.4 x 10~ 6 /°C matches very closely 
the value of 1.12 x 10 *V°C for the isotropic lay-up HMS/3501-5 laminated composite 
of Table 1. The elastic moduli also are relatively close for these two materials. 

Thus, the molded sphere will mate with the base ring of the mounting tower with 
a minimum of thermal expansion mismatch. Adhesive bonding will be quite satisfactory. 

Annular grooves will be machined in the outer surface of the molded sphere 
to carry the coil windings. Each groove must accommodate 273 turns of #36 copper 
wire (0.005" diameter), the wire having an Insulation coating resulting in a total 
diameter of 0.0055". The shape of the groove is not critical; a rectangular geometry 
is selected as being simple to machine into the sphere. The diameter of each coil 
ring is to be large enough such that the rings laying in the three orthogonal planes 
will have to cross over each other at several points. Thus, in order to maintain 
the circular shape of each coil ring, the grooves into which these coils are 
wound will have to be of three different depths. That is, if one set of parallel 
grooves is of depth b, the next set must be of depth 2b, and the third set of 
depth 3b. Thus, it is desireable to make b as small as practical. It has been 
decided that a practical upper limit to the ratio of coil winding width to depth is 
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Table 3 

Properties of hy-MAT 7534 High Modulus Graphite /Epoxy 
Structural Mat Chopped Fiber Composite [5] 


Elastic Modulus 

(10 6 lb/in 2 ) 

Axial 

12.5 

Flexural 

12.0 

Ultimate Strength 

(10 3 lb/in 2 ) 

Tensile 

35.0 

Compression 

26.0 

Flexural 

50.0 


Coefficient of Thermal Expansion 1.4 x l(f 6 /°C 


Density 

Fiber Volume 
Fiber Length 
Thickness /ply 


0.060 lb/in 3 

55 percent 
1 in. (nominal) 
0.020 in. 


Prepreg Cost 


$80/lb 




about 4. A coil winding of 35 turns width wound to a depth of 8 turns gives a ratio 
of 4*38, a reasonable value. For a coated wire diameter of 0.0055", a groove 
0.193" wide and 0.044" deep is adequate. Thus, the second set of grooves will 
be 0.088" deep, and the third set 0.132" deep. 

The wall thickness of the sphere must be adequate to accommodate these deepest 
grooves while retaining adequate strength and stiffness. A wall thickness of 0.200" 
has been selected, resulting in a minimum section thickness of 0.068" at the deepest 
grooves. The chopped fiber composite has adequate strength (see Table 3) to permit 
this. This material is also essentially insensitive to stress concentrations 
because of its fibrous nature. Thus, the presence of the groove will not Induce 
a local failure at the bottom corner of the groove. Since the highest flexural 
stresses will occur in the region of the attachment of the sphere to the base ring, 
the coil grooves oriented in the planes parallel to this base ring will be designed 
to be the shallowest. It is inanatieral from a strength standpoint as to which of 

the other two sets of grooves is the deepest, the sphere geometry being symmetrical 
in this sense. 

It should be noted that because of the different depths of the grooves, the 
resulting diameters of the coil windings will be slightly different, differing 
from one another by 0.088". Also, because of the necessity of having an opening 
in the sphere of sufficient size to permit the Insertion of the sensor core 
assembly, the maximum diameters of the pairs of coils laying in planes parallel to 
the sensor mounting tower axis is limited. That is, the minimum spacing between 
these pairs of colls is limited by the required size of the opening. Denoting the 
radius of any pair of rings as R, and the spacing between the rings as 2x, the 
minimum permissible value of x will be approximately 0.75R. 

Another important design consideration is the tension with which the copper 
wire is wound into the grooves. Two factors must be considered: the prestressing 



effect on the sphere, and the thermal expansion mismatch between the copper and the 
graphite /epoxy composite during any service temperature change. The properties 
of the copper wire are given in Table 4, the wire being used in annealed form. 

The high coefficient of thermal expansion of the copper wire relative to the graphite/ 
epoxy sphere means that the windings will tend to loosen with any increase in 
temperature, and tighten onto the sphere with any decrease. Excessive expansion 
could cause a local buckling of the wire bundle, and hence a permanent distortion 
of the coil geometry. Excessive contraction could induce a fracture of the sphere 
in the region of the groove. A stress analysis of these effects is presented in 
the analysis section of this report. 

To alleviate this thermal stress problem, the use of selective reinforcement 
of the sphere at the inner surface was considered. A thin band of graphite/epoxy 
composite bonded to the inside surface of the chopped fiber mat, graphite/epoxy 
sphere directly under each coil winding groove would enhance the local strength 
significantly. This can be fabricated relatively easily since the band material 
can be positioned inside the sphere while in an uncured (prepreg) condition, 
and then cured in place. The negative aspect of this approach is that the uniformity 
of the basic molded sphere will be distorted. In particular, the sphere will no 
longer be isotropic in terms of stiffness and thermal expansion, which could lead 
to excessive thermal distortions. While a careful materials selection and design 
would minimize this possibility, the strength reduction problem is not severe enough 
to warrant the additional fabrication complexity and performance uncertainty 
associated with using stiffness rings. Had strength reduction been found to be a 
problem, a better solution would have been to Increase the wall thickness of the 


sphere 




Table 4 

Properties of CDA* No, 102 Oxygen-Free Copper 
in Annealed Form [8] 


Elastic Modulus 
Shear Modulus 
Tensile Strength 

Yield (0.5% offset) 
Ultimate 


17 x 10 6 lb/in 2 
6.4 x 10 6 lb/in 2 
(10 3 lb/in 2 ) 
10-11 
32-35 


Elongation (in 2") 

Coefficient of Thermal Expansion 
Specific Heat 
Thermal Conductivity 
Electrical Resistivity 
Density 


45-55 percent 
17.6 x 10 _6 /°C 
0.166 BTU/lb/°C 
407 BTU/hr/f t 2 /°C/f t 
1.71 mlcrohm-cm 
0.323 lb/in 3 


* Copper Development Association 
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Field Cancellation Coils Support Sphere - The design problems associated with this 


component are essentially identical to those discussed in relation to the feedback 
coils support sphere, and the same graphite /epoxy chopped fiber composite mat will 
be used. The principal difference is that only 68 turns of the same //36 copper 
wire are required per coil. Thus, the depth of the grooves can be less. The same 
assumption of a desirable ratio of coil winding width to depth ratio of 4 or less 
applies. A coil winding of 17 turns width wound to a depth of 4 turns gives a ratio 
of 4.25, a reasonable value. For a coated wire diameter of 0.0055", a groove 0.094" 
wide and 0.022" deep is adequate. Thus, the second set of grooves will be 0.044" 
deep, and the third set 0.066" deep. A wall thickness of the sphere of 0.100" 
will result in a minimum section thickness of 0.034" at the deepest grooves. 

As in the case of the smaller sphere, the shallowest grooves will be those 
oriented in planes parallel to the base ring. Likewise, the coil windings will be 
of different diameters, differing in this case by 0.044". In order for the inner 
sphere to be inserted inside this sphere, the diameter of the opening must be 
slightly larger than 2.5". Thus, the same coil radius to spacing ratio applies for 
this larger sphere as discussed previously for the small sphere. That is, the 
minimum permissible value of x will be approximately 0.75 R. 

The interiors of both spheres will be vented to the exterior, to eliminate any 
induced stresses due to pressure differentials. 
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